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Abstract 


Environmental threats most times occur due to indiscriminate disposal of swine wastes 
into the environment however effective disposal of these waste materials into the soil as an 
amendment for vegetable production remains a vital solution for the challenge. This research 
investigated different application rates of swine wastes on selected soil physical properties 
and cultivation of Amaranthus in Awka. The research was conducted at the Soil Science and 
Land Resources Management Research and Teaching Farm, Faculty of Agriculture, Nnamdi 
Azikiwe University, Awka- Nigeria. The field experiment was laid out in a Randomized 
Complete Block Design. Soil samples were taken from each plot with core samplers 
undisturbed and randomly and was taken to the laboratory for soil physical analysis prior to 
the application of treatment. The treatment used was swine wastes at rates of 0 kg/plot (zero 
application), 4.5 kg/plot, 9.0 kg/plot and 13.5 kg/plot. The test crop used was Amaranthus 
hybridus. The results obtained from the laboratory were subjected to statistical analysis and 
the treatment means were separated using Turkey test at 5% probability level. The results 
showed that the soils of the studied area was sandy loam. There was a significant difference in 
bulk density when the soil treated with 9.0 kg/plot was compared with 0 kg/plot (zero 
application) and it ranged from 1.75-1.97 g/cm’, the total porosity of soil in the studied area 
was significant when the soils treated with 9.0 kg/plot was compared with 0 kg/plot (zero 
application) and this ranged from 23.5 - 29.59 %, the moisture content of soil in the studied 
area was significantly different when soils treated with 9.0 kg/plot was compared with 0 
kg/plot (zero application) and it ranged from 1.95- 3.22 %. The aggregate stability of soils in 
the area studied ranged from 0.37- 0.44 %, hydraulic conductivity of soils in the studied area 
was in the range of 0.01- 0.05 cm/hr. From this research, it was observed that 13.5 kg/plot of 
swine wastes applied improved the yield of Amaranthus when compared to other treatment 
rates. 

Keywords: Amaranthus, application rates, slurry, soil, swine wastes 


Introduction 
The use of organic manures for the production of leafy vegetables has been reported to 
be effective (Iren et al., 2012) and the application of these organic materials as soil nutrients 
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provides growth-regulating substrates and improves the physical, chemical and microbial 
properties of the soil (Iren et al., 2015). Most tropical soils are highly susceptible to 
degradation under continuous cultivation however good quality soil is important for 
sustainable crop production especially under continuous cultivation. Amaranthus, also called 
African spinach, Bush green, Green leaf belongs to the family Amaranthaceae. They are the 
most commonly grown leafy vegetable of the lowland tropics in Asia and Africa. Amaranth 
leaves have wonderful chemical composition with mild spinach-like taste so it comes beneath 
an accurate leafy vegetable (Amicarelli and Camaggio, 2012). It may be eaten as raw in salad 
or by cooked along or with other vegetables especially with potato; due to its chemical 
composition it plays an important role in nutrition (Uusikua et al., 2010). It contains calcium, 
Magnesium, Carotene and Niacin, Vitamin A and C are also present in significant levels. The 
protein found in young plants of Amaranths can be important for people without access of 
meat or other sources of protein thus the need to sustain the cultivation of Amaranthus using 
adequate manure. Manure from swine production could be an asset to a pork producer and to 
agriculture in general if properly managed and utilized in a sustainable food production 
system and because swine manure contains essential plant nutrients, its use as a soil 
amendment for crop production is a practical method to solve the disposal problem. The 
composition and effectiveness of swine waste as a source of nutrients depends on several 
factors including the type of ratio fed, housing system and methods of manure collection, 
storage and handling. Soil is the dynamic link between the biosphere and the lithosphere and 
constitutes a practically non-renewable natural resources, with the key role for the 
environment and for the agriculture (Moraetis et al., 2016). Some factors such as altitude, 
parent rocks, vegetation and anthropogenic activities influences the physical properties of the 
soil like soil texture, bulk density and moisture content. Inappropriate disposal of swine 
wastes into the environment often times leads to environmental pollution which could pose a 
serious health challenge hence the need to channel these wastes into the soil as an amendment 
for vegetable production remains a solution to solving the disposal challenge. 


Materials and Methods 


Description of Experimental Site 

The research was carried out at Soil Science and Land Resources Management 
Teaching and Research farm, Faculty of Agriculture, Nnamdi Azikiwe University, Awka- 
Nigeria. The research farm is located within latitude 6.244175° N and 6.250008° N and 
longitude 7.116615°E and 7.117172° E, with an average annual rainfall of 2590.6 mm as per 
Nigeria meteorological Agency data for the year 2022; the rainfall distribution pattern is 
bimodal with early rain occurring between April and July and late rain between August and 
October. It has an average relative humidity range of 76 - 80 % and average maximum and 
minimum temperatures of 32.1 and 23.5°C respectively. The experimental site has a slight 
slope with a well-drained soil. The soils of the experimental area can be classified as Ultisols 
according to USDA soil classification system. Weeds such as chromolena odorata, Aspillia 
africana, panicum maximum and cyperus spp dominated the vegetation of the experimental 
site. 
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Experimental Design and Layout 

The experiment was conducted in Randomized Completely Block Design (RCBD) with 
four treatments and three replications. The treatment used was swine wastes applied at the 
rates of Okg/plot (Zero Application), 4.5 kg/plot, 9.0 kg/plot and 13.5kg/plot. 


Field Operations and Soil Sample collection (pre-planting) 

The experimental field size of 14.5 m x 15 m (217.5 Sq./m) was mapped out using a 
measuring tape, twine rope and pegs. Prior to land clearing, surface soil samples were 
collected from the experimental plot undisturbed from ten different points using core samplers 
for physical analysis while soil auger was used to collect soil samples at 20 different points at 
the depth of 0-20 cm after which it was bulked and mixed thoroughly from which a 
representative sample was collected. Land clearing was done using a cutlass; the debris were 
packed in between the furrows. Sixteen (16) plots of 3 m x 3 m (9 Sq./m) each was made 
using a hoe, the furrows spacing was 0.5 m x | m respectively while transplanting spacing 
used was 0.25m and 0.5m respectively for intra row and inter row respectively 

Post planting 

After harvest, soil samples for physical analysis were again collected accordingly; both 
samples collected prior to treatment application and after treatment/ harvest were respectively 
labelled and bagged for onward laboratory analysis. 

Pre Nursery and Nursery Operation 

The locally sourced seeds for the experiment were subjected to viability test of which 90 
% germination was attained before planting in the nursery. A nursery bed was made very 
close to the experimental field for ease of transplant and to avoid transplanting shock of the 
Amaranthus seedlings. The swine wastes sourced from the Animal farm of the faculty of 
Agriculture was mixed in water to form slurry after which it was applied on the beds before 
transplanting. The swine wastes was directly applied on the experimental plots one week 
before transplanting with enough water applied daily to undergo a curing process. The 
Amaranthus seedlings transplanted had an average height of 0.32 m. 

Plant Data Collection 
Five middle stands from each plot were tagged and used as sample plants for data 
collection on the following agronomic properties. 

I. Average plant height at 3 weeks after transplanting (3WAT): The plant height was 
measured using a meter rule from soil level to the apex of the leaves at 3 weeks after 
transplanting. 

Il. Leaf area index: The leaf area was measured using a meter rule, it was calculated 
according to the procedure of Samkeliso et al., 2020 
Ill. Fresh shoot weight: Fresh shoot was harvested and weighed using weighing balance 
IV. Dry shoot weight: The fresh shoot was subjected to drying in the oven at a 
temperature of 105°c for 72 hours until a constant weight was achieved it was weighed 
again. 
Vv. Number of nodes after first and second cut: The nodes were counted after first and 
second cut 
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Laboratory analysis 

Particle size Distribution: The hydrometer method as described by Gee and Or 2002 
was used to determine the particle size distribution of the samples while the soil textures were 
determined using the USDA Textural triangle. The bulk density was determined by core 
method as described by Grossman and Reinsch (2002). Saturated Hydraulic conductivity (Kat) 
was determined by method of Young (2001). Soil Aggregate stability was determined using 
wet sieving method of Kemper and Rosenau (1986). Soil Total Porosity was calculated from 
the bulk density as shown in this equation: 


Total porosity (%) = = AQO) ——auideaideovanivancdic (1) 


Where Bd=Bulk density 
Pd=particle density (2.65g/ cm?) 


Moisture content of the soil was determined by oven drying at a temperature of 105°C 
and percentage of moisture in soil calculated mathematically as follows: 


GMC == x 100 Where W;=Weight of the can, W2=Weight of 


wet sample + can, W3=Weight of oven dried sample + can 


Statistical analysis 
Data collected were subjected to Welch’s one way Analysis of variance (ANOVA) using 
SAS Statistical package version 9.4. Separation of treatment means was done using turkey 
test. 
Results and Discussion 
Pre-Treatment Soil Properties 


The results of soil physical properties of the experimental site before the application of 
treatment are shown in table 1. From the results obtained, the study area had a mean sand, silt and 
clay contents of 704 g/kg, 216 g/kg and 80 g/kg respectively, having a textural class of sandy 
loam. The soils of the study area had a mean bulk density of 1.54 g/ cm’, total porosity of 41.77 %, 
soil moisture content of 9.67 %, hydraulic conductivity of 0.18 cm/hr and aggregate stability of 


0.45%. 


Table. 1. Soil Pretreatment properties of the Experimental site. 


Soil Parameters Values 
Sand (g/kg) 704 

Silt (g/kg) 216 

Clay (g/kg) 80 
Textural class Sandy loam 
Bulk density (g/ cm?) 1.54 
Porosity (%) 41.77 
Aggregate stability (%) 0.45 
Moisture content (%) 9.67 
Hydraulic conductivity (cm/hr) 0.18 
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Nutrient Properties of swine wastes 

Prior to the treatment application (Swine wastes), it was subjected to laboratory analysis 
to determine the constituent elements and the results are shown in table 2. 

Table. 2. Nutrient composition of swine waste 


Swine waste properties values 
Organic matter 3.71 % 
Total Nitrogen 3.21 % 
Available Phosphorus 1.68 mg/kg 


Post-Treatment Soil Properties 

From the results obtained as shown in table 3, the plot with zero application had an 
average sand, silt and clay contents of 729 g/kg, 206 g/kg and 65 g/kg respectively, the plot 
treated with 4.5 kg/plot of swine waste had an average sand, silt and clay contents of 719 
g/kg, 201 g/kg and 81 g/kg respectively, the plot treated with 9.0 kg/plot of swine waste had 
an average sand, silt and clay contents of 719 g/kg, 200 g/kg and 81 g/kg _ respectively while 
the plot that was treated with 13.5 kg/plot of swine waste had an average sand, silt and clay 
contents of 749 g/kg, 175 g/kg and 95 g/kg respectively. The results showed that there was no 
treatment effect on the texture of the soil; however the sandiness of the studied area could be 
attributed to the nature of the parent material from which the soils were derived. Soils of the 
experimental area belongs to textural class of sandy loam. Higher sand contents obtained in 
soils of the studied area could suggest low cation exchange capacity, leaching of nutrients and 
can encourage soil erodibility on exposure to high rainfall (Nwosu et al., 2023). The plot with 
zero treatment (control) had a mean soil bulk density of 1.97g/ cm’, the plot with treatment 
applied at the rate of 4.5 kg/plot and 9.0 kg/plot had mean soil bulk densities of 1.86 g/ cm? 
and 1.75 g/ cm? respectively while the plot with the treatment applied at the rate of 13.5 kg/ha 
had a mean soil bulk density of 1.89 g/ cm? ‘There was a significant difference in the soil bulk 
densities when the plot with zero treatment was compared with the plot of 9.0 kg of swine 
wastes. Bulk density reduces with the application of animal manures; this is because the 
organic matter has a lower particle density than the mineral particles and its increase results in 
reduction of the soil particulate density as a whole. Again, increase in soil organic matter 
increases the degree of soil aggregation, increasing the volume of pores. Soil bulk density as 
observed in this study was slightly higher than normal and could be attributed to soil 
compaction owing to machinery compression in the course of road construction around the 
experimental site. According to Celik et al., 2004; Barzegar et al., 2002, Bulluck et al., 2002 
reduction in soil bulk density occurs both with the application of manures alone and in 
association with mineral fertilizers ( Hati et al., 2008; Bandyopadhyay et al., 2010). The plot 
with zero treatment (control) and 4.5 kg/plot of swine wastes had an average mean total 
porosities of 25.72 % and 28.91 % respectively, while the plots treated with 9.0 kg/plot and 
13.5 kg/plot of swine wastes had a mean total porosities of 23.59 % and 29.59 % respectively. 
The plot with 13.5 kg/plot had the highest total porosity of 29.79 % while the plot treated 
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with 9.0kg/plot of swine wastes had the least total porosity of 23.59 %. The results obtained in 
several studies with application of manure presented divergent results regarding their effects 
on the total porosity and the classes of pores. Increase in total porosity as often observed is due 
to the increase of pores with larger diameter ie macropores which is associated with better soil 
structure (Celik et al., 2004; Ros et al., 2013). Marinari et al., (2000) observed that the 
increase of macropores in the soil treated with organic fertilizer was mainly due to an increase 
in the elongated pores which are considered very important for both soil-water plant 
relationships and maintenance of good soil structure. For aggregate stability; the results 
showed that the plot with zero treatment had 0.3 %, the plot with 4.5 kg/plot had a mean 
aggregate stability of 0.44 % while the plot treated at the rates of 9.0 kg/plot and 13.5 kg/plot 
had a mean aggregate stability of 0.40 % and 0.37 % respectively. The plot with 4.5 kg/plot of 
swine wastes had the highest aggregate stability with a mean value of 0.44 % while the 
control and 13.5 kg/plot had 0.37 % mean values respectively being the least. Soil aggregate 
stability is one of the properties most influenced by increase in manure application however 
swine wastes application at the rate of 13.5 kg/plot did not increase the aggregate stability and 
could be as a result of the short period of the experiment. According to Kay et al.,(1999) , a level 
of soil organic carbon of 2.0 to 2.5% is considered necessary to maintain good aggregate stability 
and aggregate stability is considered to deteriorate rapidly when soil organic carbon falls below 
1.2% to 1.5%. Again, the results obtained showed that the plot with zero treatment (control) 
had a mean hydraulic conductivity of 0.02 cm/hr, the plot with swine waste 4.5 kg/plot had a 
mean hydraulic conductivity of 0.01 cm/hr, the plot with 9.0 kg/plot had a mean hydraulic 
conductivity of 0.03 cm/hr while the plot with 13.5 kg/plot had a mean hydraulic conductivity 
of 0.05 cm/hr. The plot with 13.5 kg/plot had the highest mean value of 0.05 cm/hr while the 
plot with 4.5 kg/plot had the least mean value of 0.01 cm/hr. 


Table. 3. Influence of application rates of swine wastes on selected Soil Physical Properties 


TRTs BD (g/em*) TP(%) MC(%) Aggregate Ksat (cm/hr) Sand Silt Clay _Textural class 


Stability (g/kg) (g/kg) (g/kg) 
0 kg/plot 1.97° 25.72" 1.95° 0.37" 0.02" 729 20665 SL 
4.5 kg/plot 1.86°” 28.91 2.26" 0.44" 0.01° 719 201 81 SL 
9.0kg/plot —1.75° 23.59" = 3,22" 0.40* 0.03" 719 200 81 SL 
13.5 kg/plot 1.89" 29.59" 2.97% 0.37" 0.05°° 749 175 76 SL 


Note: TRTs- Treatments BD-Bulk density TP - Total porosity MC-moisture content 
Ksat — saturated hydraulic Conductivity 


There was a significant difference in soil hydraulic conductivity when plot treated with 
13.5kg of swine wastes was compared with the control. Some researchers have reported 
seepage reduction due to sealing of soils by animal manure and by other organic liquids which 
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was observed in this study. In their submission, organic manure is known to improve soil 
physical properties (Andersodu et al., 2013, Aluko and Oyeleke, (2005). 
Agronomic properties 
Some agronomic properties such as plant height, leaf area index and number of leaves, 
Fresh shoot weight, dry shoot weight, number of nodes at first and second cuts were taken 
3weeks after transplanting to determine the influence of swine wastes application rates on the 
Amaranthus plant. 
Plant Height 
From the results obtained, the plots with Okg, 4.5 kg/plot, 9.0 kg/plot, and 13.5 
kg/plot rates of swine waste application had a mean plant height of 37.25 cm, 39.65 cm, 44.90 
cm, and 49.95 cm respectively. Highest height was observed in the plot treated with 13.5 
kg/plot of swine waste while the lowest height was observed in the plot with zero treatment. 
Giacomini and Aita (2008) revealed that animal waste such as deep litter and liquid pig slurry 
applied as a source of N for corn increased the availability of this nutrient in the soil 
throughout the crop cycle compared to the control treatment. 
Leaf area index and number of leaves 
The plots with zero treatment had a mean leaf area index and number of leaves of 
37.60 cm and 43.15 cm respectively, the plot treated with 4.5 kg/plot of swine waste had a 
mean leaf area index of 46.50 cm and 60.90 cm respectively, the plot treated with 9.0 kg/plot 
of swine waste had a mean leaf area index of 43.47 cm and 73.00 cm respectively while the 
plot treated with 13.5 kg/plot of swine waste had a mean leaf area index of 70.80 cm and 
73.90 cm respectively. There was a significant difference at 5% probability level in the leaf 
area index when the plots with Okg/plot,4.5kg/plot, 9.0kg/plot swine waste application were 
compared with 13.5kg/plot swine waste application. Leaves are the main source of food 
synthesized for the plant and thus their absence affects plant growth and development thus 
Leaf area is recognized as a crucial growth index determining the capacity of plants to trap 
solar energy for photosynthesis and has marked effect on growth and yield of plant (Mathowa 
et al., 2014) 


Table. 4. Influence of application rates of swine wastes on some agronomic properties 


TRTs Plant Leaf area Number of FSW DSW NODES NODES 
Height index leaves (kg) (kg) (AST CUT) (2ND CUT) 
Okg/plot 37.25" 37.60" 43,15" 0.51° 0.15° 13.85" 4.80° 
4.5kg/plot 39.65" =. 46.50" 60.90" 0.32? 0.15° 14.75° 5.85" 
9.0kg/plot 44.90" 43.40" 73.00° 0.47" 0.12° 14.50° 4.75° 
13.5kg/plot 49.95" 70.80" 73.90” 0.51° 0.17 13.10° 4.75? 


N/B: FSW =fresh shoot weight, DSW=Dry shoot weight 
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Fresh shoot weight, Dry shoot weight, Number of nodes at first and second cut. 

The plot treated with 13.5 kg/plot of swine waste gave the highest dry shoot weight of 
0.17kg, while the plot treated with 9.0 kg/plot of swine waste had the lowest weight of 0.12kg. 
Nodes were more at the first cut than in the second cut at various application rates however 
there was a significant difference in the nodes at the first cut when plots with 4.5 kg/plot, 9.0 
kg/plot and 13.5 kg/plot were compared with zero application. Generally, it was observed that 
plant height, leaf area index, number of leaves, fresh shoot and dry shoot weight increased 
with an increase in treatment application rates. This could be attributed to high content of 
Total N, P and K in the pig manure. Adebayo et al (2011) found out that organic amendment 
increased both vegetative and dry matter yield of Moringer oliefera. 


Conclusion 

The application of swine wastes at different rates had an appreciable impact on some 
soil physical properties such as the soil bulk density, total porosity and moisture content. The 
research showed that the higher the rate of swine wastes application for cultivation of 
Amaranthus the better the growth and yield which was evidenced in the plant height, leaf area 
index, number of leaves as well as fresh and dry shoot weights measured. 


Recommendation 
Farmers growing annual crops with short growth cycle such as Amaranthus should use 
swine wastes at rate of 9kg/plot and above to obtain relative higher yields. 
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Abstract 


Due to ongoing cultivation, the Teaching and Research Farm of Kwara State 
University in Malete has seen quick and accelerated erosion. Run-off-related soil loss is a 
significant and persistent ecological problem in the study area. Information on soil loss is 
essential for promoting agricultural production and natural resource management. The 
average annual soil loss was calculated and mapped in this study using remote sensing and 
GIS. The soil loss was computed using the Revised Universal Soil Loss (RUSLE) Model. 
Using a topographic map at a scale of 1:50,000, an aster digital elevation model (DEM) with 
a spatial resolution of 20 m, a digital soil map at a scale of 1: 250,000, rainfall data spanning 
39 years (1981-2020), and other data, RUSLE's soil loss variables were calculated. The 
RUSLE parameters were investigated and incorporated using a raster calculator in the geo- 
processing tools in the arc-GIS 10.1 environment to estimate and map the annual soil loss of 
the research region. The results show that the annual soil loss in the study region ranged from 
48.553 to 1,476.606 t ha’! year’, cover around 100 ha of land. Most of the soil erosion 
affected areas are spatially situated in block 2B and 3A part of the farm. These are areas 
where low Ferric Luvisols and high Ferric Luvisols with higher soil erodibility character (21- 
33) values are dominant. Therefore, it was found that the main causes of soil erosion were 
slope gradient and length, followed by soil erodibility parameters. The study therefore 
recommended using sustainable soil and water conservation methods to address the problem 
of soil erosion in the study area. 


Keywords: Soil Loss, Soil erodibility, Slope gradient, Arc-GIS, Ferric Luvisols, 
RUSLE 


Introduction 

During the process of soil erosion, topsoil on the soil surface is lifted off the ground 
by water or wind and transferred to other surfaces. It is regarded as the biggest environmental 
problem the world is currently facing, second only to population growth. According to 
Pimentel et al. (2009), who provided examples of the United States losing soil at a rate of 10- 
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40 times faster than the average replacement rate and China and India losing soil at a rate of 
30-40 times faster, the majority of soil from farmlands is washed away about 10-40 times 
faster than it is being replaced. The trend of soil erosion increased during the 20th century. 
Soil erosion has reduced crop yield by 17%, the majority of which has happened since the 
end of World War II. It has been established in the literature that eroding soils emit harmful 
farm chemicals into streams, rivers, surface waters, and groundwater resources (N20), as well 
as atmospheric pollutants like carbon dioxide (CO2), methane (CH4), and nitrous oxide. 
Thirty percent of the world's arable land has become unusable during the past forty years due 
to erosion. 

It is still difficult to find a model that can estimate soil loss from all sources of erosion 
from a particular site. One of the parametric models utilized to calculate soil erosion is the 
Universal Soil Loss Equation. Studies using the Universal Soil Loss Equation (USLE), 
Modified Universal Soil Loss Equation (MUSLE), Revised Universal Soil Loss Equation 
(RUSLE), Water Erosion Prediction Project (WEPP), Limberg Soil Erosion Model (LISEM), 
and Soil Loss Estimation Model for Southern Africa (SLEMSA) have shown that eroded soils 
carry harmful farm chemicals into streams and rivers and contaminate surface and 
groundwater resources. The USLE (Wischmeier and Smith1978) was previously considered 
to be the best equation for calculating soil loss due to its simplicity and resilience. However, 
it has been demonstrated that when applied to areas of the world other than the United States 
of America, it is insufficient due to variations in the soil types and rainfall patterns found in 
these areas. In recent years, RUSLE has been used in several research across the globe 
(Adediji et al., 2010; Nangia et al., 2010; Sajaul et al., 2012). 

This demonstrates the method's adaptability to locations with a range of climatic 
factors and soil types. The RUSLE erosion model, which was developed to calculate the 
average soil loss over the long term due to runoff from a certain field slope, would compute 
the average soil loss with the proper selection of its factor values. The soil loss predicted by 
RUSLE is brought on by sheet and rill erosion. Sheet erosion, the initial stage of water 
erosion, is the process by which rainfall and runoff remove a layer of soil from the surface of 
the land. Contrarily, rills are often located on agricultural land and are associated with 
agricultural activity. 

There isn't much information available at the KWASU Teaching and Research Farm 
about soil loss or how likely soil erosion is to occur, and neither GIS nor remote sensing 
methods have been used to aid. Therefore, this work aims to estimate and map the 
geographical pattern of annual soil loss rate by water using the Revised Universal Soil Loss 
Equation (RUSLE) simulated using GIS and Remote sensing techniques. 


Materials and Methods 

The Teaching and Research Farm at Kwara State University in Malete (KWASU 
T&R), which is located in the Moro Local Government Area of Kwara State, served as the 
study's experimental location (Figure 1). 
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Figure 1. Location Map of the study area. (Source - Alabi et al., 2017.) 


It is located eight kilometers north of the agricultural buildings and covers about 100 
hectares. The farm is located 360 meters above sea level, between Latitude 08°71'N and 
Longitude 04°44'E, and the terrain is fairly gentle. It is situated in Nigeria's Southern Guinea 
Savannah ecological zone. It consists of a farm center and a space for teaching and research 
by the departments of crop production, agricultural economics, and extension services. It also 
includes a section for animal production, which includes fisheries and aquaculture, small and 
large ruminant, poultry, and rabbitry. It is characterized by distinctive wet and dry seasons 
with a mean annual rainfall of about 1150 mm, with a double maximal pattern between April 
and October. The dry season starts in November and lasts until March, whereas the wet 
season starts in April and concludes towards the end of October. The average yearly 
temperature ranges from 25 to 28.9° © degrees. The land area is a part of the Nigerian 
basement complex's South-Western zone of basement reactivation. The area is predominantly 
used for the cultivation of arable crops such as maize, groundnut and cowpea, perennial trees 
such as cashew and mango. The site also contains a variety of grasses, including spear grass, 
elephant grass, and Guinea Gamba grass, as well as woody species including the baobab 
(Adansonia digitata), neem tree (Azadractha indica), and acacia (acacia species). 
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For this investigation, both primary and secondary data sources were used. Primary 
data were collected using GPS devices through field surveys, ground truth authentication, and 
observation. This gives precise information about what is happening in the study area. The 
six blocks that made up the farm were Crop Museum, Block 2A, Block 2B, Block 3A, Block 
4, and Operation Feed Yourself field. To produce primary data pertaining to ground truth for 
the purpose of training the image utilizing supervised image classification and producing 
thematic land use and land cover maps, intensive field observations using GPS were 
conducted in each block of the farm. For each key land-use/cover category, ground control 
points (GCPs) were also acquired in order to validate accuracy. Additional data includes the 
soil map (1:250,000), the slope length-steepness (LS) factor, a 2019 Thematic Mapper (TM) 
multi-spectral image with a spatial resolution of 30 m (land-use/land-cover map), and 39 
years (1981-2020) of rainfall records from the Lower Niger Basin. These data sources were 
all downloaded from the Global Land-cover Facility (www.landcover.org). The Google Earth 
image was also made and digitized. Each layer of the RUSLE parameters was discretized and 
structured to the DEM's (20x20) cell size in order to obtain a fine grid-based soil loss result 
and to derive a precise slope length-steepness factor value. Digital technology and the Google 
Earth image were used to produce a map of the research area's waterbodies. Published and 
unpublished resources, such as research papers, census reports, and journal articles that were 
gathered from various sources, as well as the aforementioned secondary data, were also used. 

The Revised Universal Soil Loss Equation (RUSLE) was used to calculate the mean 
yearly soil loss that took place in the Teaching and Research Farm using the enhanced LS 
factor estimation approach. The empirical formulation of the Revised Universal Soil Loss 
Equation (RUSLE) is: 


A(M t/ha’ yr')= R* K* LS*C#P oo ceeeeeeeee Equation 1 


where A is the mean annual soil loss (M t/ha-' yr-'); R is the rainfall erosivity factor (MJ mm 
h-' ha-' year-'); K is the soil erodibility factor (M t/ha-’ MJ-' mm-); LS is the slope length— 
steepness factor (dimensionless); C is the cover and management factor (dimensionless); and 
P is the erosion support practice or land management factor (dimensionless and ranges from 
zero to one). The RUSLE model was simulated by GIS and Remote sensing techniques. 
While the raster calculator geo-processing tool, the R-value of each grid cell or block was 
determined from this continuous rainfall data. 


RY 623. DP careaxnaccunameaarceneanapeaxacienes Equation 2 
Where; 


Req= Equivalent R factor for unique climatic condition\ 
Pr= annual precipitation, mm 
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Results and Discussion 

Assessment of the rate of soil loss 

In this work, the Revised Universal Soil Loss Equation (RUSLE) model was used to 
calculate the mean annual soil loss rate (t/ ha-' yr-') cell by cell. The Teaching and Research 
Farm in Malete (Figure 1) also discovered and documented soil erosion danger zones. The 
following are a few of the raster maps that were made and discussed for each RUSLE 
parameter: 


Rainfall erosivity (R) factor 

According to Equation (2), the erosivity factor is estimated to be between 87.391 and 
91.086 MJ mm h- ha-' yr-'. The rainfall erosivity factor, which also considers the amount 
and rate of runoff that is most likely to be associated with precipitation events, measures the 
influence of rainfall. The Lower Niger Basin Meteorological Station in Ilorin, Kwara State 
provided the mean annual rainfall for the research area, which was interpolated in the ArcGIS 
10.1 environment using the inverse distance weighted approach (Table 1) to create unbroken 
rainfall data for each grid cell. 
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Figure 2. Rainfall Erosivity (R) factor 
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Soil erodibility (KX) factor 


The natural resistance of soil particles to the detaching and transporting power of 
precipitation is a sign of soil erodibility (Wischmeier & Smith, 1978). This component 
assesses a soil type's cohesiveness and resistance to dislodging and transport caused by the 
force of overland flow shear pressures and raindrop impacts. For a particular soil type, the K- 
factor is empirically computed to represent the physical and chemical properties of the soil 
that affect its erodibility potential (Figure 3). The Teaching and Research Farm's erodibility 
K-factor ranges from 10 to 33 mt/ha-' MJ’ mnr’ (Table 1). The highest erodibility factor 
values were reported by Ozcan et al. (2008) using USLE and GIS, which is consistent with 
the findings of our work. Koulin et al. reported similar findings in the Chania watershed in 


Greece in 2009. 
Table 1. Estimated K Value for Each Block 


Estimated K _~ value 


Block Soil Name/Class (metric tons ha-’ MJ-' 
mm-') 

Crop museum High ferric Luvisols 21 

2A High to Medium ferric Luvisols 10-21 

2B High to Medium ferric Luvisols 10- 21 

3A Low ferric Luvisols 33 

4 Low ferric Luvisols 33 

Operation feed ; ; : 

Yourself field Medium ferric Luvisols 10 
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Figure 3. Soil Erodibility (K) factor 
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Slope length and steepness (LS) factor 

Using Simms et al., (2003), the topographic component of RUSLE was calculated. 
Slope length was replaced for the upslope contributing area in order to account for the flow 
convergence and divergence in a three-dimensional complicated terrain situation. The (LS) 
factor, under otherwise identical conditions, is the ratio of soil loss per unit area from a field 
slope to that from a 22.13 m length of uniform 9% slope (Wischmeier & Smith, 1978). In this 
case, the slope-length (L sub-factor) refers to the distance between the beginning and end of 
the inter-rill process. The culmination occurs when the flow concentrates into a rill or other 
man-made channel, such as a terrace or diversion, or when the slope begins to drop and the 
subsequent depositional process begins (Wischmeier & Smith, 1978; Renard et al., 1996). 
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Figure 4. Slope Length and Steepness (LS) Factor 


Slope determination utilizing a geographic information system (GIS) has benefits for 
a wide range of environmental models because slope properties are typically needed as input 
for landslides, land planning and construction, and other environmental models (Dunn 
&Hickey, 1998). the length of the slope-drawbacks. When applying the Universal Soil Loss 
Equation, the cumulative uphill length from each cell can be used to overcome the limitations 
of slope-length computation by taking into consideration convergent flow pathways and 
depositional areas (Hickey, 2000). The supposed transport capacity of the flow is measured 
by the LS factor in the RUSLE (Moore & Wilson, 1992). Without accounting for the intricate 
three-dimensional geometry of the landscape. The supposed transport capacity of the flow is 
measured by the LS factor in the RUSLE (Moore & Wilson, 1992). The slope-length 
calculation approach employed in USLE and RUSLE was based on the notion that the longer 
the slope, the greater the soil loss (Robert & Hilborn, 2000). This was done without taking 
into account the complicated three-dimensional geometry of the terrain (Figure 4). The results 
of slope length indicate that the prolonged areas of steep slopes are susceptible to severe 
erosion. In the Turkish Kazan watershed, Erdogan et al. (2007) used USLE and GIS models 
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and reported similar, reliable results. Crop Museum and Operation Feed Yourself Field's 
topography favoured less erosion in general. This suggests that the land is seriously 
endangered, which affects its output. High correlation between steep slopes and poor surface 
cover within the farm was also reported by Kouli, 2009 in the Crete watershed. 


Support practice (P) factor. 

The erosion control practice factor (P-factor) is the ratio of soil loss caused by a 
certain support practice to losses comparable to those caused by upslope and downslope 
agriculture (Wischmeier & Smith, 1978). The generated P-factor map (Figure 5) is not used 
to comprehend the conservation practices being applied in the research region. The 
preventative measures that affect drainage patterns, runoff concentration, and runoff velocity 
and minimize the erosive power of rainfall and runoff are considered in this component. The 
classed land-use/land-cover and slope thematic map format was converted to vector format 
and the corresponding p values were added to each combination of land-use/land-cover and 
slope classes in order to build a raster map of the P factor. 
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Figure 5. Support Practice (P) Factor 


Cover and management (C) factor. 

The cover and management (C) factor is defined as the ratio of soil loss from land 
with a specific vegetation to the corresponding soil loss from a continuous fallow 
(Wischmeier & Smith, 1978; Morgan, 2005). It is the component that can be changed the 
most and is typically taken into consideration when developing a conservation strategy. The 
major method for gathering data on the dominant categories of land use and land cover in the 
study area was unsupervised categorization. On the basis of this data, supervised 
classification supported by GCPs was used to produce themed land-cover maps. The thematic 
land-use and land-cover raster map of Teaching and Research Farm, Malete was then 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 21 


Alabi et al., 2024 


transformed to vector format in order to assign the proper cover and management factor value 
derived from various studies. The result was a C-factor raster map, as shown in Figure 6. 
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Figure 6. Cover and Management(c) factor Map 


The primary land-uses/covers of the Teaching and Research Farm that were identified 
using supervised image classification were water bodies and cultivated land. Maize, cowpeas, 
groundnuts, melon, guinea corn, and other vegetation dominated the area. Crop covers with 
the greatest C-factor values (0.1). This ensures that annual crops like millet and maize do not 
mitigate the direct effects of rainfall on soil resources, unlike forest land. Thus, the final map 
that shows the anticipated annual soil loss of the Teaching and Research farm was produced 
by superimposing the aforementioned five parameters (K, R, LS, C, and P) using Eq. (1), 
raster calculator geo-processing tools, and the ArcGIS10.1 environment (Figure 7). 

The cell sizes for the DEM were resampled (20x20 m), and each layer was organized 
in a grid style. The soil erosion map of the study area was divided into six classes for the 
purpose of management prioritization, as shown in Table | and Figure 7. The level of soil loss 
risk in the research region has also been classified and the potential amount of soil loss has 
been predicted using the statistical approach. As shown in Figure 7, the annual soil loss of 
Teaching and Research farm extends from 0 to 1,476.606 M t/ ha-' yr-’ with a mean annual 
soil loss rate of 47.4 Mt/ ha-' yr-'. The total annual soil loss in the Teaching and Research 
farm is 1476.606 t; which covers about 100 ha of land. This lost is at a rate much greater than 
the tolerable soil loss rate (Table 2). 

Table 2. Annual Soil loss range, Block coverage and severity class 


Annual soil loss Soil erosion Percent of 
: Block coverage ; 
(tone) risk class total soil loss 
0-48.553 Low Crop museum 2.02 
48.553-140.031 Moderate Operation feed yourself field 5.83 
140.03 1- 280.367 High 2A and 2B 11.68 
280.367- 455.398 Very high 4 18.97 
455.398-1476.606 Severe 3A 61.5 
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Figure 7. Soil loss rate Map 


Tolerable soil loss rate (FAO, 1984) is defined as the highest allowable soil loss rate 
that can support economic growth and high levels of production. Since the majority of the 
primary soil types on the Teaching and Research Farms in Malete have soil depths between 0 
and 150 meters, the FAO's (1994) tolerable soil loss rate for deep and very deep soil depths is 
4.2 to 7.2 t/ ha-’ year-'. This suggests that the land is seriously endangered, which affects its 
output. This implies that the land is gravely in jeopardy, which has an impact on its 
productivity. Additionally, the system that keeps life from dying could deteriorate until it 
reaches an irreversible condition. The Teaching and Research farm's yearly soil loss ranges 
from 0 to 1,476.606 Mt/ ha-! yr, with a mean rate of 47.4 metric t/ ha-' year, In the 
Teaching and Research farm, the total yearly soil loss is 1476.606 t, which occupies a space 
of around 100 ha. Compared to the permitted rate of soil loss, this loss is occurring far more 
quickly than the allowable rate of soil loss (Table 2). Tolerable soil loss rate (FAO, 1984) is 
defined as the highest allowable soil loss rate that can support economic growth and high 
levels of production. Since, most of the major soil types of Teaching and Research farm, 
Malete have soil depth ranges from 0 — 150 m deep (Alabi et al., 2017), the soil loss rate in 
the study area is far beyond FAO (1994) tolerable soil loss rate of 4.2—7.2 t/ ha-' yr-' for deep 
and very deep soil depth respectively. Kumar et al., (2022) reported an average annual soil 
loss of 1.60 t. ha’. yr! in the lower Chambal basin in Madhya Pradesh using RUSLE model. 
Adhikary et al., (2014), Suryawanshi et al., (2023) reported similar results by comparing 
USLE and RUSLE models in Central Indian Chambal River Basin. 


Conclusion 

Information about the quantitative and spatial extent of soil loss derived by modeling 
of RUSLE parameters Malete ensures the management of spatially variable data and 
inaccessible location where ground-based observation is challenging by using GIS and 
Remote sensing techniques in Teaching and Research farm. As a result, the method can be 
used in other areas of the community for the evaluation and identification of regions that are 
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prone to erosion, the setting of conservation priorities, and the evaluation of the effectiveness 
of different land management strategies. 

The results of this study take into account the erosion risk map and the spatially 
distributed soil loss rate of Kwara State University's Teaching and Research Farm in Malete. 
The annual soil loss at Malete's Teaching and Research farm ranges from 0 to 1476.606 
metric tons. Slope-length and gradient (LS) factor was shown to be the most important 
RUSLE parameter by the results as well. According to the soil erodibility (K) factor, the 
potential soil loss is therefore frequently greater along the steeper slope and in areas with less 
vegetation cover. 

The areas used for Operation Feed Yourself and Crop Museum Block have been 
discovered to be the least prone to soil erosion. Therefore, in order to reduce soil erosion, the 
research area's conservation and management measures are strongly suggested. Areas where 
significant soil loss was detected should be given careful thought. An appropriate land use 
plan and an Environmental Impact Assessment (EIA) are required for the locations that have 
been made accessible on the University's Teaching and Research farm. For immediate soil 
protection efforts in their agricultural regions, local residents and researchers in the area 
should embrace various soil protection approaches, such as mulching, contour plowing, 
repeated cropping, and other indigenous soil conservation strategies. Additionally, using 
enough vegetative cover and adding organic manure should be considered. 
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Abstract 


Makurdi is in a strategic position in the agricultural map of Nigeria, producing a wide 
range of both annual and perennial crops such as yam, maize, rice, sorghum, soybean, 
cowpea, citrus, mangoes, and a variety of vegetables. One of the factors responsible for this 
wide range of crops is the favorable climate. This study was aimed at evaluating the suitability 
of the soils for the production of groundnut and to have a detailed soil database for effective 
land use planning. Soil requirements for groundnut were obtained from past research works 
and compared with data collected from the survey. The study showed that the soils of the area 
had formed under climatic environment presently characterized by an annual rainfall of about 
1330.20 mm and a mean annual temperature of about 27.80 °C. The soils were well drained to 
poorly drained. The clay content ranged from 7.20 to 29.30 %, increasing with depth. Organic 
carbon was low (0.64 %) in the upland, but relatively high when compared with the low (0.62 
%) value in the low land. The suitability assessment results showed that although, certain 
quantities or characteristics such as mean annual temperature, relative humidity and base 
saturation were optimum for groundnut cultivation, there was however, no highly suitable (S1) 
land for groundnut cultivation in the area. All the soils were classified into moderately 
suitable (S2¢) subclass due to their low nutrient levels. COATRF I and II units of the area were 
moderately suitable (Sx) due to topography and low soil fertility. COATRF III was limited by 
its imperfect drainage to marginally suitable subclass (S3y) for groundnut production. To 
improve on the level of productivity of the soils for optimal groundnut production, 
management techniques such as continuous organic matter incorporation and mineral fertilizer 
application, and efficient use of mineral fertilizers with low levels of chemical inputs with 
adaption of appropriate irrigation techniques would make dry season farming sustainable. 
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Introduction 

The broad principles of land evaluation involve comparing the requirements of land use 
with quality of land thereby assessing the value of each type of land present for each land use 
(Dent, 1981). Decisions on land use are being based on comprehensive analysis of the 
production and potentials of natural resources such as climate, soil and hydrology (Abagyeh, 
2016). Soil evaluation is very important in this direction as it provides information on the 
potentials and constraints of land for a defined land use type in terms of crop performance as 
affected by the physical environment. Soil evaluations are based on knowledge of crop 
requirements, prevailing conditions and applied soil management methods (Ande, 2011). Soil 
evaluation for suitability classification quantifies in broad terms to what extent soil conditions 
match crop requirements under a defined input and management (FAO, 1976). Evaluating the 
soil enables optimum performance and maximum productivity of the crop. In evaluation, the 
specific crop requirements will be calibrated with the terrain and soil parameters (Dent, 1981), 
so that the identified limiting factors could be managed to suit crop requirements and improve 
productivity. Land evaluation thus, enables management guidelines in order to promote a 
more sustainable use of the soil and environmental resources (Maniyunda, 2007). 

Groundnut (Arachis hypogaea L.) is a member of sub-family, papilonaceae of the family 
Leguminosae. It is believed to be the native of Brazil to Peru, Argentina and China. It was 
introduced into India during the first half of the sixteenth century from one of the pacific 
islands of China. Groundnut oil is an edible oil. It finds extensive use as a cooking medium 
both as refined oil and Vanaspati Ghee. “They are rich in protein and vitamins A, B and some 
members of B2 group. Nigeria is the largest groundnut producing country in West Africa, 
accounting for 51 % of production in the region. The country contributes 10 % of total global 
production and 39 % that of Africa. Between 1956 and 1967, groundnut was the country’s 
most valuable single export crop, exemplified by the famous Kano groundnut pyramids. 
Groundnut is a major source of edible oil as well as livelihoods for small-scale farmers in 
Northern Nigeria. Being a labor-intensive crop, it generates employment for the rural poor. It 
is planted on about 34 % of total cultivated area and contributes to 23 % of household cash 
revenue. Groundnut products like oil and cake accounted for a significant percentage of total 
Nigerian export earnings. Before the fossil oil boom, groundnut was one of the major sources 
of revenue and foreign exchange earnings. However, in the post-1967 period, the combined 
effects of drought, increasing prevalence of diseases such as rust, leaf spots and groundnut 
rosette disease (GRD) have caused a decline in groundnut production. The total output of 
groundnut in 1970 was 1.6 m tons, but fell to 0.47 m tons in 1980. Due to insufficient 
groundnut stocks, processors and marketers in Kano State source groundnut from as far as 
Chad Republic. 

The year-round demand for groundnut means farmers can increase production without 
any fear of market glut. Since 1984, production has been increasing at an estimated growth 
rate of 8 %, resulting both from area expansion (6 %), and increase in productivity of 2 % 
(Ndjeunga et al., 2010). Hence, the study was to evaluate the soils and climatic conditions 
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supporting groundnut (Arachis hypogaea L.) production in College of Agronomy Research 
Farm of University of Agriculture Makurdi Benue State 


Materials and Methods 

Study area 

The study was conducted within the teaching and research farm of the College of 
Agronomy University of Agriculture, Makurdi. The University lies North-east of Makurdi 
town which is between Latitude 7.45° N and 7.50° N and Longitude of 8.00° E and 9.00° E. 
Makurdi is in a strategic position in the Agricultural map of Nigeria, producing a wide range 
of both annual and perennial crops such as groundnut, sorghum, soybean, cowpea, rice, maize, 
yam, citrus, mangoes and a variety of vegetables; one of the factors responsible for this wide 
range of crops is the favorable climate (Igomu and Idoga, 2017). Makurdi has tropical climate 
conditions with distinct wet and dry seasons. The dry season starts from November to March 
while rainy season starts from April and end in October. In the past, Benue state experiences a 
bimodal rainfall distribution pattern with one peak in July and another in September. Recent 
meteorological data point to unimodal pattern with the single peak in August. The mean 
annual rainfall for Benue state varies between 1000 and 1600 mm. The mean minimum air 
temperature of Makurdi is 16.20 °C to 17.20 °C during the period of harmattan, which is 
December to February. The mean maximum air temperature is 37.7 °C in March prior to the 
onset of rains. The slope of the area is 0 to 5 % and the mean elevation above mean sea level 
is about 93 m (Ojanuga and Ekwoanya, 1995). 


Field studies 

Field studies involved taking soil samples using a core sampler as well as the collection 
of bulk soil samples for laboratory analysis. Soil samples were collected using grid method 
with a core sampler based on the nature of the slope of the study plot of 600 X 400 meters. 
Five samples were collected to form composite samples from each of the upper crest location, 
the mid slope and foot (lower) slope, respectively to enhance proper study of the soil of the 
areas. The soil samples were described using the pattern outlined in the Soil Survey Manual 
(2010). 


Laboratory analyses 

Standard laboratory procedures were employed in the investigation of the physical and 
chemical characteristic of the soil samples collected from the field. The bulk soil samples 
collected from the field were air dried. The air-dried samples were gently crushed using 
mortar and pestle. The samples were sieved through 2 mm and 0.5 mm sieves and used for 
laboratory analyses. For laboratory analyses like particle size distribution, Boyoucos 
hygrometer method was used. Electrometric method as described by Hesse (1971) was used 
for pH analysis. Walkley-Black method as described by Hesse (1971) was used for organic 
determination. Total nitrogen was determined by the use of macro Kjelidahl method 
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(Bremmer, 1965). Sodium bicarbonates NaHCO; extracting solution Olsen (1965) was used 
for available phosphorous determination, and exchangeable cations (Ca, Mg, K and Na). 
Cation exchange capacity was also determined using neutral solution of ammonium acetate 


method. K and Na in the extracts usually determined using Flame photometer while Ca and 
Mg are determined using AAS. 
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Figure 1. Map of Benue State, Nigeria. 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 31 


Igomu et al., 2024 


Land evaluation procedures 

The conventional parametric method was used for the suitability evaluation of the land 
(FAO 1976), where soils were placed in suitability classes by matching their characteristics 
with the established requirements for groundnut (Table 1). The aggregate suitability class 
indicates the most limiting characteristics of the land. The parameters used for the land quality 
calculation include rainfall, mean annual temperature, slope, wetness, drainage. While those 
of soil characteristics included texture, coarse fragment, soil pH. Other fertility indicators 
were base saturation, organic carbon and available phosphorus. 
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Figure 2. Map of Makurdi Local Government area of Benue State. 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 32 


Igomu et al., 2024 


Table 1. Soil and Climatic Requirement for Suitability of Groundnut (Aracis hypogaea) 
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fimaatne © 
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Annual rainfall 600-1500 | 500-1250 450-850 
Mean Annual Max Temp. ee ee a 


Relative humidity 50-80 | >80 | - J - 
Length of dry season 120-150 110-130 90-130 roe | 


es ee eee eae ee 
Wemess | «dG Ag | 6-30 | 930 
FO Poor 
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Coarse fragment (0-50cm) 


Fertility (1) | 


a 
CEC mo kg) clay | >24 | 16) | <t6) | - 
Base saturation | % | 50 | 20-35 | 20 | 
pH Water | 59-700 | 55-80 | 50-80 | - | 
Organic carbon(OM) | _(0-50)% | >2_— | 08-12 | os | 
fvainble Pome fate [ter 


Key: FO-No Flooding, FI- Seasonal Flooding, FR- Flooding Rare, C-Clay, CL- Clay Loam, LS- Loamy 
Sand, SL- Sandy Loam, LCS- Loamy Clay Sand, CS- Clay Sand, S- Sand, S1- Highly suitable, S2- 
Moderately suitable, S3-Marginally suitable, N 1- Currently not suitable. 


SOURCE: Adesemuyi (2014). 
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Results and Discussions 

The particle size distribution values for the soil mapping units are shown in Table 2. 
The values revealed the sandy nature of all the soil units. Soils of units I and III generally 
have loamy sand surface texture while, soil unit II had sandy loam surface texture. The sand 
fraction of the three soil units ranged from 70.64 to 71.08 %. This high fraction of sand could 
be attributed to the nature of the parent material (sandstone). Sand fraction decreased with 
depth in all soil units. The clay content of the soils was moderate to high ranging from 14.92 
to 16.20 %, while silt fraction ranged from 13.0 to 14.0 %. Because of illuviation, the clay 
content of the soils was higher in the mid crest than the top crest. The weak fine to medium 
crumbs of upper crest of all soil units (Table 2) could be attributed to their low level of 
organic matter and sandy nature of the soils. In the entire field, the three units had moderate to 
strong, medium to coarse sub-angular blocky structures. This may be attributed to the 
corresponding increase in clay content in the sub soils. 

The pH (H20) of the soils ranged between 6.10 and 6.16 (Table 2) indicating strong to 
slightly acid (USDA, 2010). There was a general tendency among these soils to have higher 
pH values at surface, which may decreased with depth. Ojanuga and Ekwoanya (1995) 
attributed the slightly higher pH in the upper most soil horizon to the accumulation of 
exchangeable calcium (Ca) and magnesium (Mg) derived from plant roots and litter 
decomposition, as well as biogenetic cycling of bases. The pH values of the soil units I and III 
were higher than those of Unit II. As the cations are leached, hydrogen ion (H*) from water 
(H20) replaces the leached cations in the exchange complex to balance the charge. Soil acidity 
is known to affect most soil nutrients especially P, N, S, and micronutrients (Agbede, 2009). 

The organic carbon (Table 2) was low according to the ratings of Esu (1991) with 
values ranging from 0.59 to 0.64 %. These values may decreased with depth in all the three 
Soil Units due to the concentration of plant roots and residue on the soils. The low organic 
carbon content of these soil units could be attributed to annual bush burning in the area as 
well as the native vegetation cover and climate. According to Igomu and Idoga, (2017), 
drainage appeared to have strong influence on the organic carbon content of the soil. 
Earthworms and other organic fauna increase organic carbon content of the soil by breaking 
and biochemically altering fresh organic matter through mixing with inorganic materials and 
later excreting them to form a moist homogenous blend of the organic and mineral matters, 
which prevents rapid loss of humic compounds (Esu, 1982). More earthworm casts were seen 
on the surface of soil unit II than units I. Total nitrogen values ranged from 0.13 to 0.18 % for 
all the soil units studied. The distribution of total nitrogen follows the same pattern as organic 
carbon in topo sequence. Osujieke et al. (2018) revealed that savannah soils are generally low 
in total nitrogen. This is due to the low organic carbon content of the soils of the savannah 
agro-ecological zone. It is obvious that climate, vegetation and human activities contribute to 
the low level of nitrogen. The available phosphorus of the soil had values ranging from 2.8 to 
3.2 mg/kg and were generally low when compared with the ratings of Chude et al. (2011). 
The surface values of the amount of available phosphorus, which may decreased with soil 
depth shows that much of the available P is in organic forms associated with organic matter. 
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The values of Ca**, K* and Na* were low while Mg** was moderate when compared 
with the ratings of Landon (1991) in all the soil units. The low values of exchangeable bases 
of these soils maybe due to intensity of weathering and lateral translocation of bases. The 
amount of exchangeable calcium (Ca) ranged between 2.5 and 3.0 cmol (+)/kg of soil. 
Magnesium (Mg) ranged between 2.3 and 2.7 cmol (+)/kg of soil. Potassium (K) ranged 
between 0.24 and 0.27 cmol (+)/kg of soil while Sodium (Na) ranged between 0.21 and 0.24 
cmol (+)/kg of soil. The predominance of Ca over the other cations in these soils may be 
because of all the exchangeable cations, Ca is least easily lost from the soil environment 
(Brady, 1974). The CEC of the soil ranged from 6.37 to 7.31 cmol/kg of soil. The lowest CEC 
values were recorded in unit III owing to their low organic carbon content and perhaps low 
clay content. Donahue et al., (1976) reported that some sandy soils of the tropics with 
dominant clay might have low CEC hence, low fertility. Base saturation by sum of cations 
ranged between 82.42 and 85.51 % which is rated very high (Idoga and Ibrahim, 2012). The 
high base status showed that the soils had high native fertility, which is confirmed by the 
luxuriant growth of the vegetation of the soil. 


Land Suitability Evaluation 

The overall methodology followed involved the selection of attributes relevant to 
groundnut production, generation of data layers, weighting and overlaying of the data layers 
to generate a land suitability map for groundnut production. Elevation dataset was used as 
proxy for climatic variables (rainfall and temperature) and it was generated by reclassifying 
the digital elevation model (DEM). The distance dataset representing accessibility was 
generated using the Euclidean distance tool. Each data layer was reclassified into suitability 
classes based on groundnut requirement. The suitability levels assigned to each dataset were 
defined based on literature, expert knowledge, observation and practical experiences. The 
suitability levels were ranked as highly suitable, moderately suitable, marginally suitable, 
currently not suitable and permanently not suitable based on the suitability classification 
structure (Adesemuyi, 2014). The datasets generated were finally combined using the 
weighted sum overlay to generate the land suitability map for rainfed groundnut production. 
The matching of land characteristics used for suitability ratings of sites for groundnut 
(Arachis hypogaea.) production in Makurdi (Table 3) with land and soil requirements for 
suitability rating of groundnut (Table 1) resulted in the suitability classes of the soil mapping 
units as shown in Fig. 3. The region was optimal or near optimal in mean annual temperature, 
relative humidity, length of dry season, slope and base saturation. 

Textural classes of the COATRE I soils ranged from loamy sand to sandy loam, thus, 
considered moderately suitable (S2s) for the cultivation of groundnut. The soil texture for 
optimum groundnut performance is sandy loam or loam (Sys, 1985). However, topography 
and low soil fertility levels limited the soils into the moderately suitable (S2tf) subclass for 
groundnut production. Based on these assessments, none of the soils was considered highly 
suitable for groundnut cultivation in the area. Soils of COATREF I and II were classified into 
moderately suitable (S2f) subclass as a _ result of their low nutrients _ status 
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Conclusions 

This research offered some valuable insight with outcomes vital to groundnut 
production in Federal University of Agriculture Makurdi, Benue State of Nigeria. The 
suitability assessment results showed that although, certain qualities or characteristics such as 
mean annual temperature, relative humidity and base saturation were optimum for groundnut 
cultivation, there was however, no highly suitable (S1) land for groundnut cultivation in the 
area. Soils of COATRF I and II was classified into marginally suitable (S2f) subclass, due to 
their low nutrients’ levels. COATRF I and II Units of the area were moderately suitable (S2tf) 
due to topography and low soil fertility. It can further be concluded that the generated land 
suitability map with overall classification accuracy which was 320 meters of the farm land can 
rightly be used as a guide in decision making for groundnut production. Therefore, it is 
conclusive that the soils of College of Agronomy Research Farm Federal University of 
Agriculture Makurdi, is moderately suitable for groundnut production. 
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Abstract 


Intermediate crops used for plant cover or green manure affect a positive direction of 
different soil indicators: organic matter content, soil structure, reduced bulk density and soil 
strength, limit the risk of erosion, eventually leading to a reduction in soil and moisture losses. 
The purpose of the present study was to establish the role of year-round cover of arable land 
with crops of cereals and forage crops and agro-technical measures to improve some water- 
physical parameters of Alluvial-meadow soil (Fluvisol) in crop rotations. 

The study was carried out in the experimental field of Tsalapitsa, Plovdiv at the 
experimental base of “N. Pushkarov” Institute of Soil Science, Agrotechnology and Plant 
Protection, using the standard block method. In the next crop rotation - maize (Zea mays, L.) — 
barley (Hordeum vulgare, L) — maize (Zea mays, L.), were included intermediate crops 
(winter and spring) grown for green mass before the main spring crop maize (Zea mays, L.). 
From the obtained results, it was found that when in crop rotations were included winter 
precultures the risk of negative influence of climatic factors on the physical characteristics of 
the Eutric Fluvisol is reduced. Despite leaving the stubble from the pre-sowing crop, after 
spring mixture crops, there is an increase in the values of bulk density and soil strength. These 
results are probably due to the poor condition of the soil that was plowed in the spring and 
insufficient time to carry out the sowing. The study confirms the role of soil cultivation in 
improving the studied physical parameters - moisture content, bulk density, and soil strength. 


Key words: intermediate crops, crop rotation, soil water-physical parameters 


Introduction 

The sustainable management of soils in agricultural systems places great demands on 
the research related to the mechanisms for maintaining all quality parameters of the soil 
(physical, chemical and biological) in optimal limit values (Dumanski, 1991; Krastanov, 
2006). The main task is to achieve long-term maintenance of soil fertility through the use of 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 42 


Gerasimova et al., 2024 


natural remedies and organic materials - organic fertilization, plant residues, green 
fertilization, applied within a well-structured crop rotation (Mitova et al., 1999; Estrade et al., 
2010, Dimitrov et al., 2013; Karlen et al., 2013, Grifith et al., 2013). These agricultural 
practices ensure the return to the soil of a significant amount of biodegradable biomass and 
are the basis for the construction of highly productive ecological systems, together with their 
positive significance for limiting ecological risks. In this aspect, the intermediate crops used 
for plant cover, which occupy the fallow period before the spring crops and are sown to be 
incorporated into the soil (green manure) or left as mulch, positively change the various soil 
indicators - the organic matter, the soil structure and the stability of the soil aggregates, bulk 
density, soil strength, reduce the risk of erosion, which leads to reduced losses of soil and 
moisture (Marinari et al., 2000; Mitova et al., 2015). According to research by Tsvetkova et 
al. (1995) green manure in conventional crop rotation improves basic soil physical parameters 
and reduces erosion losses. In addition to the above-ground biomass, a significant amount of 
root mass remains in the soil, which is subject to mineralization. According to (Talgre et al., 
2010), the root mass of leguminous crops for green manure constitutes 37-54% of the total 
biomass. The results of long-term field experiments have been used to increase the effect of 
their management, for the structural formation, stabilization of the composition and 
improvement of the water-physical parameters of the soil (Ulrich et al., 2006; Dimitrov et al., 
2011; Abdollahi et al., 2014; Kinoshita, 2017). The provision of year-round plant cover of the 
arable land makes it possible to protect the soil surface from the risk of erosion processes and 
preserve soil fertility (Milev & Iliev 2014; Mitova et al., 2011; Laird et al., 2013). Apart from 
being a soil protection agent, plant residues from the cultivated main and intermediate crops in 
field crop rotations are a source for increasing the content of organic matter in the soil. 
(Pommeresche et al., 2006; Gerasimova, 2018). 

The present study aims to establish the role of year-round plant cover of the arable 
land with crops of cereal and forage crops included, crop rotation with cover crops as well as 
the agrotechnical measures for their cultivation to improve some water-physical parameters of 
Alluvial-meadow soil (Fluvisol). 


Material and Methods 

The researched area was located at the experimental field of ISSAPP "N.Pushkarov" in 
the village of Tsalapitsa, Plovdiv region, South Bulgaria for the period 2016-2018. 

The experiment is based on the block method — on long plots on an area of 0,3 ha, in 
four replications with a size of the experimental plots of 40 m?. The study included year-round 
coverage of the soil surface with crops of cereals and grass forage crops and a control variant 
without used coverage crops. During the study period crop rotation included cover crops for 
vegetation before the main spring trench crop — maize (Zea mays). 

The following 4 variants (Aj; A3; As and Ag) were included in the study with 
alternating preculture - maize (Zea mays) — barley (Hordeum vulgare) — preculture - maize 
(Zea mays): 
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Table 1. Crop rotation 2016-2018 years 


Variants 2016 2017 2018 
Preculture (winter peas Preculture (winter 
V Aq | +barley) - maize Barley peas + barley) - maize 


Preculture (spring peas 


V A3_ | + spring barley) — maize Barley Maize 
Control - maize 

VA5 without pre-crop Barley Maize 
Perennial wheat-bean Perennial wheat- Perennial wheat-bean 


mixture (sainfoin (60%) + | bean mixture (sainfoin | mixture (sainfoin (60%) 
VAg_ | hedgehog head (20%) + | (60%) +hedgehog | + hedgehog head (20%) 
timothy (20%) head (20%) + timothy | + timothy (20%) 

(20%) 


Variant A, - preculture (winter peas (Pisum sativum) + barley (Hordeum vulgare)) - 
maize (Zea mays) for grain — barley (Hordeum vulgare) — preculture (winter peas (Pisum 
sativum) + barley (Hordeum vulgare)) - maize (Zea mays) for grain. 

Variant A3 - preculture (spring peas (Pisum sativum) + spring barley (Hordeum 
vulgare)) — maize (Zea mays) for grain - barley (Hordeum vulgare) - maize (Zea mays) for 
grain. 

Variant As - control — maize (Zea mays) for grain without pre-crop - barley (Hordeum 
vulgare) — maize (Zea mays) for grain. 

Variant Ag - perennial wheat-bean mixture (sainfoin (Onobrychis sativa) (60%) + 
hedgehog head (Dactylis glomerata) (20%) + timothy (Phleum pratense) (20%). 

The variants with winter cereal-bean mixture and variants with spring cereal-bean 
mixture are included in crop rotation. Variant As is used for control. A perennial grass field 
(Ag) was also included in the study, which was also used as a control variant. 

According to the stock of soil with nutrients after agrochemical analysis in both 
studied crops for barley, was used fertilizer norm Nj29P i099 per ha, and maize — Nj40P joo. 

For winter cereals phosphorus and nitrogen fertilizers are imported in the autumn with 
the last pre-sowing tillage; for barley nitrogen fertilizer is introduced at 1/3 rate pre-sowing 
and 2/3 of the norm as feeding at the end of the ,,tillering“ phase; for maize — phosphorus is 
imported for preculture (except for the control variant — pre-sowing), and nitrogen - with pre- 
sowing tillage. 

The tillage system includes plowing + spring tillage (shallow plowing or chiseling) for 
the main crop and pre-sowing tillage for cereals with a heavy disc harrow. The main tillage 
for maize compacted with winter mixtures was carried out in autumn, and the one compacted 
with spring mixtures in spring, in order to leave unploughed stubble from the predecessor - 
barley. 
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The following water-physical parameters have been determined: 

- moisture content of the soil (%) - by the weight method, as a percentage of the 
weight of the soil, through 10 cm, at a layer depth of up to 60 cm; 

- bulk density of the soil in g/cm? - according to Kaczynski's weight method, with 
rings of 100 cm’ per layer depth up to 40 cm, through 10 cm layer thickness. 

- soil strength — with penetrometer with falling weight type DORNII, with a cone 
angle of 30°, and then 5 cm. in depth to 40 cm; 

The climate in the area is hot with dry summers and mild winters and irregular 
rainfall distribution (Levicharska, 1991). The period accepted as a norm 1961-1990 has an 
amount of precipitation of 492 mm. The average annual air temperature is 12.5-12.8 °C, in 
July — 23.4 °C, and in January it is around 0 °C. In the study area, the soil is Alluvial-meadow 
(Eutric Fluvisol), (USS Working Group WRB, 2015). This type of soil is used for intensive 
agriculture. The soil profile has a coarse texture, low water holding capacity and fairly high 
water permeability, which creates conditions for active migration of chemical elements along 
the profile. The arable horizon has the following mean characteristics: light texture, low clay 
content (18.6%), bulk density between 1.54 and 1.66 g.cm®, slightly acidic reaction pH H»O 
(6.0), total nitrogen (0.052%), humus content - 1.23%, and cation exchange capacity — 22.4 
cmol.100 g”' (Stoichev, 1997). 

The weighted average value of the capillary field capacity for the layer 0-100 cm is 
19.6%, the wilting moisture is 11.1%, and the productive water is 8.4%. With these values of 
the indicators, the one-meter soil layer retains at field capacity a total of 306 m°/da of water, 
176 m° of non-productive water, and only 131 m? of plant-usable water reserve. According to 
the indicator of usable water reserves, this soil type can be referred to as poorly stocked. 

The mathematical-statistical analysis of the experimental data was performed with the 
statistical program SPSS. 


Results and Discussion 


During the experimental period, significant dynamics were established in the values of 
the studied physical parameters of the soil type Alluvial-meadow (Eutric Fluvisol). 
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Figure 1. Monthly amount of precipitation in 2016-2018 and for a long-term period 
(1935-1985), (l/m?) 
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Soil moisture content. 

The content of soil moisture is a basic physical parameter on which the dynamics of 
the processes taking place in the soil depend, as well as the emergence, growth, and 
development of cultivated crops. During the study period, lower precipitation values were 
found for the study years compared to the long-term precipitation database (Figure 1). 

After the experiment was set in the autumn of 2015, low content of soil moisture in the 
arable layer was reported - around wilting moisture (4-6%). No significant differences in 
study variants were found. In the spring of 2016, soil moisture was measured before the 
harvest of the winter and after the emergence of the spring mixture crops and the beginning of 
the vegetation in the perennial mixture. The data shows that as a result of the rainwater, the 
values for the parameter are significantly higher than those reported in the autumn. The 
content of soil moisture has the highest values in the variant with preculture - winter mixture 
peas + barley (Aj) - from 12.34 to 15.98% (Figure 2). Compared to the variant (As) — without 
preculture, the differences in the 20-60 cm layer are from 1.40% to 2.14%, and in the 10-20 
cm layer the trend is the opposite. The obtained lowest soil moisture was in the variant with a 
perennial mixture crop, which can be explained by a higher consumption of moisture than the 
plants that have started their vegetative development. 


OA1 
BA3 
B45 
BA6 


o 5 10 15 20 
Sail moisture, % 


Figure 2. Content of soil moisture by weight (%) - during the spring period of 2016 


The moisture content in the soil marks significant changes after the sowing of maize in 
the variants with spring mixtures and in the crop stage "3-5 leaf", in the variants with winter 
mixtures, and in the without preculture variant. The humidity was increased by 5.5% to 
7.25%, except for the surface 0-10 cm layer, where with drying the values increased to 4.0%. 
The highest values for the parameter were found in the variant with winter mixtures, as in the 
layer 30-40 cm the differences compared to the other variants were statistically proved with a 
probability of error p <1%. According to research by Mitova et al., (2011), the presence of 
plant residues in the variants with an intermediate crop for green manure helps to limit the 
loss of moisture in the high summer temperatures (Mitova et al., 2011). 

The soil moisture has the lowest values in the variant without preculture - 17.21% to 
21.60% (Figure 3b). 

During the maize tasseling phase (June, July) the rainfall was a little 5 to 15 mm, 
during which drought occurred (Figure 1). As a result of the drought that occurred after the 
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tasseling phase (for variants A3 and As - phase 9-10 leaves) the moisture content in the soil 
marks a steady downward trend (Figure 3a). The largest differences in the values compared to 
the previous reading were found in the layers 0-20 cm and 40-60 cm of the Eutric Fluvisol. In 
the plowed layer, the low values for the parameter are the result of increased evaporation from 
the increased atmospheric and soil temperature, and the deepest layer of the studied profile is 
probably due to the movement of moisture by capillary. 


Depth.cm 


0 5 10 15 20 10 20 30 
Sail moisture, % Sail moisture, % 


oO 


Figure 3. Content of soil moisture by weight - maize 2016 a) phase “tasselling”’; 
b) phase “3-5 leaf” 


In the layers 10-20 cm and 30-40 cm, the differences in soil moisture between the 
variant compacted with winter pre-crops (A,) and the other variants are significant with a 
probability of error p <5%. Most significant in this phase is the influence of crop rotation on 
the percentage of moisture content in the soil — 92,49% (Table 2). 

Up to the wax maturity phase, the maximum values are 11.74% and 11.83% reported 
in the 30-40 cm layer of the compacted variant (As). The lowest moisture content is in the 
deeper layers of the studied profile in the variant with the perennial mixture (A¢), which is 
probably due to the small amount of rainwater and the higher consumption of the root system 
of plants. 


Table 2. Dispersion analysis of moisture content — maize ’tasselling” phase 2016 
Sum of | Sum of | Value of| Medium | F-rate |Sign.Lev 
variation Squares | Sq. % the ch. Sq. el(p) 


Depth (D) 4 
Tillage (O) 
Rotation (S) _| 2936.501 008" 
ores 

O12# 

0,25 010 ** 


333 - 


517 | 137 [| 40 | 1083 | | 
Total Sum___| 3174885] | 59 | 


SPD 5%=1,722 
SPD 1% =2,303 
SPD 0,1%=3,025 
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In the next crop in crop rotation — barley (2017), the reported moisture content in the 
“tillering” stage is 60-70% of maximum soil moisture, except for the surface layer 0-10 cm, 
where it is below 50% of maximum soil moisture. The low values for soil moisture are due to 
the lack of precipitation in the spring (Figure 1). The comparison between the variants in the 
amount of moisture shows higher values in the variant without preculture - from 5.95% in the 
surface 0-10 cm layer to 12.16% in the layer 40-60 cm (As). There is evidence of differences 
in the layers 30-40 cm and 40-60 cm compared to variant A3, and only in the deepest profile 


layer in variant A Z (Figure 4a). 
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Figure 4. Soil moisture content in weight % barley 2017 a) phase “‘tillering’ 
b) phase “head emergence of a class” 
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Table 3. Dispersion analysis of moisture content —barley “tillering” phase 2017 


| Sum of | Sum of | Value of | Middle F-rate |Sign.Lev 
variation __| Squares Sg. Yo ch. Sq. el(p) 

Depth (D) 4 

Tillage (O) 

Rotation (S) 


| 400 | 40 | 0668 | | 
Total Sum___| 662.847 | | 59 OE 


SPD5% =1,344 
SPD1% =1,798 
SPD 0,1%=2,361 


When taking into account the soil moisture in the “head emergence of a class” phase of 
barley, higher values were reported in the layer 0-10 cm, as the highest is in the variant barley 
after spring mixtures (A3)- 13.5% (Figure. 4b). Lower values of moisture were reported in the 
layer 30-40 cm, which is due to poor water permeability as a result of soil compaction. In the 
variant with perennial grasses, uniformity in the moisture content in the layer of 10-60 cm is 
established, which is due to the lack of moisture inflow from the deeper (below the measured 
profile) soil layers. 

In the next 2018, in stage "3-5 leaf" of maize, a deficit of soil moisture is established. 
This led to opposite trends compared to previous years in the measured values by survey 
variants. The parameter in variant A3 (after spring mixtures) has larger dimensions, probably 
due to the earlier phase of maize development ("second leaf" stage) and lower consumption of 
productive moisture by the plants. Differences in the measured values at depth over 40 cm in 
this variant compared to the other variants were proved at p <1%. This trend is maintained in 
the phase of "tasseling" of plants from the variants with spring mixture and control. This fact 
confirms the assumption that in variant A3 the higher values for soil moisture content are due 
to lower consumption by less developed plants. The lowest values are reported in the variant 
with perennial grass mixture - between 4% and 6%, which is due to drying as a result of 
elevated temperatures (Figure. 5a). Another reason is the stronger compaction of the soil in 
the whole studied layer due to the lack of tillage. Most significant in this crop stage is the 
influence of crop rotation on the percentage of moisture content in the soil (Table 4; Figure 
5a). 
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Figure 5. Content of soil moisture by weight % - maize 2018 
a) phase “3-5 leaf”; b) phase “tasselling” 


Table 4. Dispersion analysis of moisture content — maize” tasselling” phase 2018 
variation __| 
Depth (D) 


Tillage (O) 
Rotation (S) 


Total Sum 


SPD 5% =1,166 
SPD 1% =1,560 
SPD 0,1%=2,049 
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The content of soil moisture in the arable layer in the phase of "tasselling" of maize 
has values below the wilting moisture - 7.40%, as a result of the higher consumption of the 
enhanced vegetation of the crop and the loss due to evaporation after inter-row cultivation. 
The lowest values were reported in the layer 0-10 cm (A¢— 3.91%) and 10-20 cm (Ag — 5.7%), 
(Figure 5b). The lower humidity values are due to the lack of precipitation during the 
vegetation period and from the vegetation treatments (Figure 1). 

Maize yield was reduced by water stress, in the different stages of development, but 
the culture is very sensitive during the phases of tasseling and detection of the cob. 


Bulk density of the soil. 
Bulk density is one of the most frequently measured parameters, primarily as an 
indicator of the presence of soil compaction. 
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Figure 6. Bulk density of soil in g/cm’ - maize 2016 a) phase “3-5 leaf”; b) phase 
“tasselling”’ 


The bulk density of the soil measured in the phase "3-5 leaf" of maize in the first year 
of the experiment was differentiated by layer. The surface layer has a value between 1.3 g/cm? 
- 14 g/em’, then slightly increases in the layer 10-20 cm, increases its values in 20-30 cm 
(1.51-1.59 g/cm?) and reaches the highest value in the layer 30-40 cm - 1.66 g/cm? in A3 
(Figure 6a). The bulk density of the variants with used cover crops was lower compared to the 
variant without compaction 

In the "tasseling" stage of maize, the measured values for soil bulk density remain 
almost unchanged and show a steady upward trend from the surface to the deepest layer of the 
soil profile. Although statistically unproven, there are differences between the values of the 
parameter in the variants with preculture (A; and A3) and the one without preculture (As) in 
the 0-30 cm layer (Figure 6b). This fact shows that compaction of the area with preculture 
helps to maintain the values of bulk density in a lower range. 

The bulk density determined during the barley vegetation is high. In the “tillering” 
phase, they show an upward trend from the surface to the deepest layers. 
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Figure 7. Bulk density of soil in g/cm? —barley 2017 a) phase “tillering” b) phase 
“end of heading: inflorescence fully emerged” 


Depth, cm 


The values are lower in the variants A; and Az (barley after maize with preculture 
winter and spring mixtures) in the layer 0-10 and 10-20 cm bulk density is 1.41 - 1.43 g/em> 
and the layer 20-30 respectively to 1.49 - 1.52 g/cm® (Figure 7a). The differences compared to 
variant As (control) in the layer of 20-30 cm were proved at p <1% and p <5%. Most 
significant in this crop stage is the influence of crop rotation on the percentage of moisture 
content in the soil — 47.74% followed by soil depth — 28.60% (Table 5). In stage ,,end of 
heading: inflorescence fully emerged “for all studied variants, the obtained data for bulk 
density are higher than those reported in the ,,tillering “stage”. In all studied variants of this 
vegetative stage, the highest obtained values were recorded in 10-20 cm depth and the lowest 
in the 0-10 cm layer. In the 0-10 cm layer, the lowest value was obtained in variant As and the 
highest in A6. At a depth of 30-40 cm, the lowest value was recorded in variant A3 and the 
highest in Ag (Figure 7b). 


Table 5. Dispersion analysis of bulk density — barley “tillering” phase 2017 


Source of Sum of |Sum_ of|Value of| Middle | F-rate |Sign.Lev 
variation Squares |\Sq.% ch. Sq. el(p) 
-0OO*** 


Depth (D) 3 
Tillage (O) 
Rotation (S) 
“O21** 
| 7.62 | 32 | oo01 | | | 
Total Sum a a 


SPD 5% =0,052 
SPD1% =0,072 
SPD 0,1%=0,093 


In the barley “tillering” stage, this trend is also found in variant A;, while in A3 in the 
arable layer, the bulk density increases to 1.71 g/cm’, which is due to the formation of a thick 
soil crust due to drought. The highest values were reported for the variant with perennial grass 
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mixture from 1.68 g/cm? to 1.74 g/cm’. As a general trend showed, lower values are reported 
in the variants with the growing of the crops in three-field crop rotation. 

In the third year of the study, the bulk density determined during maize vegetation was 
high. The values show an upward trend from the surface to the deepest layers (Figure 8a). 
The lowest values for the bulk density in the vegetative stage "3-5 leaf" were obtained again 
in the variants with winter and spring mixtures, as in the layer 10-20 cm the differences 
compared to the control were proved with a probability of error p <1%. 
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Figure 8. Bulk density of soil in g/cm — maize 2018 a) phase “3-5 leaf”; b) phase 
“tasselling”’. 


Statistically proven (p <0.1%) is the difference between the measured bulk density in 
the layer 0-30 cm between the variant with perennial crops and those with winter mixtures, 
and the differences reach 0.14 g/cm’. 

In the ,,tasseling” crop stage of maize as a result of the deficit of soil moisture, the bulk 
density in the deeper layers increases significantly. The lowest value for all studied variants 
was reported in the A; (variant with compaction with winter mixture) in 0-10 cm layer. The 
highest value is the bulk density of the area with the perennial mixture from 1.59 to 1.71 
g/cm*, as compared to the control (Figure. 8b). There is a proven difference only in the layer 
0-10 cm, which is 0.07 g/cm’, compared to variant A,. There are significant differences in the 
layer 0-30 cm (Table 6). 

From the dispersion analysis of the data, it was established that all the studied factors 
have an influence, the most significant being the influence of the depth of cultivation, the 
alternation of the crops, and their interaction. During the study period, the layer 0-20 cm with 
the largest share of active root mass and the presence of microorganisms are the variants with 
compaction of crop rotation. As a result of properly conducted agro-technical measures, the 
bulk density is maintained in values up to the critical limits, despite the lower soil moisture 
during the growing season. Compaction of the area with pre-crops helps to maintain the bulk 
density with lower values. 
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Table 6. Dispersion analysis of bulk density — maize ’tasseling” phase 2018 


Source of [Sum  of| Sum. of |Value of| Middle | F-rate |Sign.Lev 
variation |Squares | Sq. % ch. Sq. el (p) 


Total Sum 
SPD 5% =0,074 
SPD1%=0,100 
SPDO,1%=0,132 


Soil strength 

Soil strength is a parameter that is closely correlated with soil moisture and bulk 
density. In the Alluvial-meadow (Eutric Fluvisol) in Tsalapitsa the humus horizon is shallow 
and with low organic content, and the fraction of coarse sand predominates in the soil 
structure. Therefore, at greater depths, the soil strength is often higher. During the study 
period the values in the treated layer up to 30 cm are within the optimal range, whereas there 
is no risk of impeding the penetration and function of the roots of cultivated plants. 

When growing maize for the first year in the crop rotation, as a general tendency it was 
found that the soil resistance was maintained within the limits of the optimal values for the 
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Figure 9. Soil strength in kg/cm2 - maize 2016 a) phase “3-5 leaf”; b) phase 
“tasseling”’. 


In the "3-5 leaf, phase, the soil strength from 27.0 to 39.15 kg/cm? in the surface layer 
(0-10 cm) increases in depth and reaches values 2-2.5 times higher than indicated (Figure 9a). 
In the variant with winter preculture in the 15-25 cm layer, the values are lower. 

As the moisture content decreases in the “tasseling” phase, the values increase by 
about 10.30 kg/cm? compared to those reported in the spring. The soil in the layer of 10-40 
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cm of the variant with perennial grass mixture has the greatest strength. This fact proves the 
role of tillage as a regulator of its physical condition (Figure 9b). The compaction with 
precultures in variants A; and A; has a positive effect on the values of the parameter in the 
layer 10-20 cm. The soil strength reported in the barley “tillering” phase (second year of the 
study) ranged from 13.78 kg/em? to 104.61 kg/cm’. The hardness of the soil in the layer of 15- 
25 cm increased most significantly in all tested variants, probably due to strong soil 
compaction. At low moisture content, the strength of the soil aggregates and the adhesion 
forces between them increase significantly. At the end of the vegetation in barley after maize 
(spring mixtures) the increase of the resistance forces is the most significant — over 120 
kg/cm’. In the variant barley, after maize compacted with winter mixtures, the reported values 
of soil strength are the lowest compared to the control and other variants, probably due to the 
full coverage of the soil surface with plants throughout the year. 

During the spring period of the third year of the study in the "3-5 leaf" stage of maize, 
the reported values for soil strength ranged from 31.78 kg/cm? to 94.61 kg/cm? (Figure 10a). 
In the treated layer up to 20 cm, the values are within the optimal limits, where there is no risk 
of impeding the penetration and function of the roots of cultivated plants. 

Due to strong soil compaction in the depth 5-10 and 20-25 cm layers, the hardness 
increases significantly - from 72.40 to 95.08 kg/cm? (Ag). In maize the variant after winter 
mixtures (A;), the reported values of soil strength are the lowest compared to the other 
variants (on average by 20.96 kg/cm’), probably due to the fuller coverage of the soil surface 
with plants during the period of severe drought. Statistically proven (p <1%) is the difference 
in the layer 20-30 cm. 
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Figure 10. Soil strength in kg/cm’- maize 2018 a) phase 3-5 leaf; b) phase 
, tasseling “. 

In the ,,tasseling“ phase of maize, there is a differentiation of values depending on the 
factors tested. The tendency of lower values of the parameter in the variant with winter 
mixture shows that the coverage of the area during the autumn-winter period has a better 
effect on maintaining the soil strength with lower values than the stubble (variant As with 
spring mixture). In the control - As without compacted in a layer of 15-25 cm - 108.40 


5 


Depth, cm 


kg/cm’. The strength of Eutric Fluvisol remains with the highest values in the variant with the 
perennial mixture (A6), and in a layer of 20-30 cm, it is over 120 kg/cm’, i.e. in the layer with 
the strongest drying (Figure 10b). 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 55 


Gerasimova et al., 2024 


Conclusions 

From the study was established that the investigated factors affect the general physical 
state of the soil. The year-round plant cover of the soil surface with the sowing of crops 
provides the most stable maintenance of the values of the studied physical parameters in the 
optimal range. The compacting with winter pre-crops the risk of negative influence of climatic 
factors on the physical characteristic of Eutric Fluvisol is reduced. After the spring mixture, 
there was an increase in bulk density and soil infiltration resistance, possibly due to the poorer 
soil condition resulting from late tillage for the following spring crop. As a result of properly 
conducted agro-technical measures, the bulk density is maintained in values up to the critical 
limits, despite the lower soil moisture during the growing season. Compaction of the area with 
intermediate crops helps to maintain the bulk density with lower values. During the research 
period, the soil strength values in the arable layer up to 30 cm. are within the optimal limits, 
and there is no risk of hindering the penetration and functioning of the roots of the cultivated 
plants. The study confirms the role of soil cultivation in improving the studied physical 
parameters - moisture content, bulk density, and soil strength. 
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Abstract 


Soil moisture plays a great role in physical, chemical and biological activities in the 
soils. A research was conducted to determine the soil moisture contents across three land use 
types (grassland (GL), forested area (FR) and cultivated land (CL)) over a period of three 
months during rainy season in Federal University Gashua, North East of Nigeria. Soil samples 
was taken from each land use type using soil auger at 0-15 and 15-30cm depths for physico- 
chemical properties analysis, while core samples were taken fortnightly and after every rainy 
event for gravimetric soil moisture content (SMC) and bulk density determination. Results 
showed that the soils across the land uses are sandy loam in texture, high bulk density at GL 
and CL (>1.80g cm™) and moderate in FR (1.45g cm”), very low organic matter content and 
neutral in pH value. Significantly high SMC was recorded in FR and lower in GL with a high 
coefficient of variation (43%). Temporally, the SMC varied moderately (35% coefficient of 
variation) with significant increased after every rainfall event. The high bulk density and very 
low water retention capacity of the soils of the study area can be enhanced by conservation 
farming practices and constant addition of organic matter. 


Keywords: Gashua, moisture, north east, savanna, soil. 


Introduction 

Physical, chemical, and biological processes of the soil system as usually influenced 
by soil moisture content (Brevik et al., 2015). Water contained in the soil is called soil 
moisture, its presence and movement are key variables for different applications in the fields 
of agriculture and vegetation growth monitoring, climate system, hydrology and stream flow 
prediction (Seneviratne et al., 2010; Kumar et al., 2014). Understanding soil water variations 
under different land use types is critical for agricultural water management and land erosion 
control (Gao and Shao, 2012; Ziadat and Taimeh, 2013). 

Soil moisture is one of the key factors limiting the growth of plants most especially in 
the arid and semi-arid environment where the crop water requirements in most cases are 
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beyond the water supply, but even in the humid environment where poor rainfall distribution 
and water management results in occasional water stresses (Musa and Adeoye, 2010). 

Land use change was found also to contribute to the variability of soil moisture, 
leading to soil deterioration, agricultural productivity decline and land degradation (Biro et al., 
2013; Leh et al., 2013; Parras-Alcantara and Lozano-Garcia, 2014). Evaluation of the 
temporal and spatial variation of soil moisture in numerous land uses is important for 
vegetation restoration and watershed management (Li et al., 2018). The soil internal factors 
affecting the spatial and temporal changes of soil moisture were mainly roots, soil texture, and 
soil organic carbon, and reported that soil clay content was the major factor affecting the 
spatial distribution of soil water on a semi-arid land (Gao and Shao, 2012). The temporal 
variation in soil moisture was also driven by meteorological properties like seasonal changes 
in temperature and evapotranspiration (Jia et al., 2017). In arid and semi-arid regions, the 
most important limiting factor to sustainable crop production is soil moisture. This valuable 
resource for sustainable agriculture must be monitored, conserved and managed well. 
Knowledge of the spatial and temporal distribution of soil moisture and the impact of land use 
types on its availability is indispensable for correct management strategies. In some tropical 
regions the availability of reliable data for water retention in respect to soil type, texture and 
organic matter content is low, it therefore became desirable to analyze the soil so as to 
estimate its water retention capacity through moisture content determination within the 
selected land uses for sustainable water management. The study was conducted to determine 
the soil moisture variation on the selected land use types in Gashua, Yobe State, North East, 
Nigeria. 


Materials and Methods 

The study was conducted at Federal University Gashua, Bade Local Government Area 
of Yobe State, Nigeria. Bade LGA is located between longitudes 12° 30° and 12° 65‘N and 
latitudes 10° 35‘and 11° 15‘E (Figure 1). The area is situated in the Sahel savanna agro- 
ecological zone of Nigeria. It has an area of 809.661 km? with a population of 139,804. The 
most important vegetation type is that of the Sahel savannah, which consists mostly of open 
thorny savannah with short trees and grasses (Alhassan et al., 2018). 
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Figure 1. Location map of the study area 
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Geology and soils of the study area 

The area is located within the Chad basin, lithological description consist of 
continental depositions, which was deposited unconformable over the Precambrian crystalline 
basement rock (Agada et al., 2020). The soils of Bade LGA showed a dominance of sand 
particles and most of them belong to the textural class of sandy loam with very low level of 
most of the essential nutrients (Alhassan et al., 2018). 


Climate of the study area 

The climate is characterized by high atmospheric temperature and low seasonal 
rainfall. The mean minimum temperature ranges between 10-12°C usually in 
December/January, while the mean maximum is about 35-42°C mostly in March to May. The 
annual average rainfall is about 300-450 mm per annum and is unimodal and lasts from July 
to October while the dry season starts from late October to May (Alhassan et al., 2020). 


Soil sampling techniques 

The sampled land use type’s locations and coordinates were cultivated land 
(11°900°90"N, 12°52’57"E), forested area (11°00’95"N, 12°54’56"E) and grass land 
(11°00’75"N, 12°53’11"E). Three soil samples collected using soil auger from each of the 
selected land use types at 0-15 and 15-30 cm depths at random. A total of 18 soil samples 
were then collected from the three land use types. Composite samples were made for each 
sampling depth thus, arriving at six composite samples which were then package separately in 
polythene bags according to depth and site for processing and laboratory analysis. 


Soil analysis 

The samples were air-dried separately, ground using a pestle and mortar, and sieved 
through a 2 mm sieve for particle size analyses, pH, EC, organic matter and other chemical 
properties. Particle size distribution was determined by hydrometer method. pH (in water) was 
determined potentiometrically using pH meter in a 1:2.5 soil: water ratio. Electrical 
conductivity (EC) was determined from soluble salts in an extracted solution of 1:2.5 soil: 
water ratio using electric conductivity meter. Organic carbon (OC) was determined by 
Walkley-Black dichromate wet oxidation method, percent of OM was obtained by multiplying 
percent soil OC by a factor of 1.724. Total nitrogen content of the soil was analyzed by micro- 
Kjeldahl method. Exchangeable bases (Ca, Mg, K and Na) were extracted with 1.0M 
NH,OAC solution buffered at pH 7.0. Calcium and Magnesium in the filtrate were 
determined with an Atomic Absorption Spectrophotometer, while K and Na with Flame 
Photometer. Cation exchange capacity (CEC) was determined by neutral (pH 7.0) NHsOAc 
saturation method. Available phosphorus (AP) was determined using after the extraction of 
the soil samples with 0.5M sodium bicarbonate (NaHCO3) following the Olsen extraction 
method (Estefan et al., 2013). 


Determination of soil moisture, soil bulk density and soil porosity 

Soil samples for soil moisture determination were taken also at 0-15cm and 15-30cm 
depths using core sampler. Samples were taken fortnightly and two days after every rain to 
determine the temporal trend and effect of rainfall respectively on soil moisture content under 
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the selected land use types. Samples were taken from ha August to 31“ October, 2021. The 
samples were then weighed for moist weight (ww) and placed in an Oven to dry at a 
temperature of 105°C for 24 hours. The dried samples were then weighed to obtain the dry 
weight (wd). It is then obtained using the expression in Hillel (2014): 


SMC (%) = ww ta) wil (2), 100 


wd (g) 

(Eq.1) 

Where: SMC (%) = soil moisture content in percentage, ww = weight of wet soil (g), 
wd = weight of dry soil (g). 

Bulk density was determined on individual samples by dividing the oven-dry weight of 
the sample by the volume of the soil probe at the respective depth (Reeuwijk, 1992). 

Bulk density (pb) == 

(Eq.2) 

Where: Ms = is the soil oven- dried mass, Vt = volume of the soil core sampler. 

Porosity was determined using equation as (one minus the solid volume fraction of a 
sample (the bulk and particle densities). The ratio of the bulk density (Db), to the particle 


density (Pp), describes the fraction of the total volume occupied by solids (Hao et al., 2008). 
ae eee 
TP(%) = [1 (2) x 100 
(Eq.3) 
Where; TP = total porosity, pb = bulk density, ps = particle density (average value of 
2.65 g cm” is used as particle density). 


Data analysis 

Statistical analyses performed to test the influence of soil moisture on different land 
use using ANOVA (Analysis of Variance), and comparisons of the difference between mean 
water content in different land uses were separated by least significant difference (LSD) 
method using STAR statistical package (STAR, 2013). 


Results and Discussion 

Effect of land use types on soil physical properties 

The soil particle size distribution (PSD) showed that the sand (p=0.08) and silt 
(p=0.21) fractions did not significantly differ among land use types (Table 1) and clay 
(p=0.04) was significant different with the land use. The highest overall clay content was 
recorded under forest land (16.00%) and the lowest in the grass land (13.00%). High clay 
content in the forested area may be due to the nature of the erosional sediments deposited. The 
soils of the study area is sandy loam in texture. Similarly, Shehu et al. (2015) reported that 
most of the soils in savanna zone of northern Nigeria had sandy loam texture. 

Soil bulk density (BD) (p=0.0001) found to be significantly different with the land 
uses. The soil bulk density was found to be very high in GL with the mean value of 1.85g cm 
> which is at par with CL (1.83g cm”) and lower bulk density was recorded in FR with the 
mean value (1.35¢ cm”®). The lower bulk density found in FR could be due to high organic 
matter content as the area is covered with trees and grasses, while very high BD in GL and CL 
it may be due to low organic matter content. Result in GL and CL is very high and could 
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restrict root penetration as reported in Alhassan et al. (2018) that normal range of bulk density 
for clay is 0.90 to 1.40g cm” >and a normal range for sand is 1.40 to 1.90 g cm” with potential 
root restriction occurring at > 1.40 g cm’? for clay and > 1.60 g cm” for sand, he also reported 
mean bulk density value (1.63 g cm”) for soils of Bade LGA in Yobe State, Nigeria. The 
mean bulk densities recorded are relatively high for crop growth, a desirable bulk density for 
the soil should be < 1.5 g cm® and bulk densities greater than 1.6 g cm™ usually restrict root 
growth (Brady and Weil, 2013) 

Porosity (p=0.0001) significantly differs within the land use types, the lowest porosity 
was recorded in GL 30.19% and CL 30.82% and the highest porosity was found in FR 
49.06%. Higher porosity found in FR may be due to the high organic matter in the land use, 
low bulk density and less disturbance. The low porosity recorded in GL and CL it may be due 
to high bulk density and intensive farming and grazing activities on these lands. This 
corroborates the work of Abdullah et al. (2018). 


Table 1. Effect of Land Use on Soil Physical Properties 


Land use Sand Silt (%) Clay Textural BD Porosity 
(%) (%) Class (g/em*) (%) 

GL 61.00 26.00 13.00b SL 1.85a 30.19b 
CL 60.00 25.50 14.50ab SL 1.83a 30.82b 
FR 55.50 28.50 16.00a SL 1.45b 49.06a 
SE 2.27 1.67 0.97 0.0186 0.7038 
Pr(>F) 0.08 0.21 0.04 0.0000 0.0000 
CV% 5.46 8.84 9.48 6.38 11.02 


SL = Sandy Loam/ BD=Bulk density/GL=Grass Land/ CL= Cultivated Land/ FR=Forested Area/ 
SE=Standard Error/ CV=Coefficient of Variance 


Effect of land use on soil chemical properties 

The results showed that pH (p=0.19), EC (p=0.42), OC (p = 0.92), Na* (p=0.87), 
available phosphorus (AP) (p=0.18) where not significantly different with the land use types, 
while OM (p= 0.01), CEC (0.00), N (p= 0.02), Mg** (p= 0.01), Ca** (p= 0.00) and K’, 
(p=0.02) concentrations where significantly different within the land uses (Table 2). 

The soil pH showed little variability across the land uses with a mean value of 6.98 in 
GL, 6.76 in CL and 7.10 in FR, the soils are then categorized as neutral as the pH range 
between 6.51 and 7.5 (QldGov, 2013). The pH values of 6.5 to 7.5 is considered ideal for crop 
production, therefore, the soil pH values of the study area is considered good for crop 
production and will not require any amendment (Estefan et al., 2013). 

Organic carbon (OC) was not significantly different with the land use types (p=0.92). 
The mean OC recorded are 0.90 g/kg in grass land, 0.93 g/kg in cultivated land and 0.87 g/kg 
in forested area. The OC content within the soils were considered low according to Esu (1991) 
rating (<10 g/kg) and is in agreement with Shehu et al. (2018) who reported less than 0.8% 
OC within the top 30 cm in Southern Guinea Savanna of North Central Nigeria. The soil OC 
with (CV=17.11%) is moderately varied within the land uses. Bationo and Buerkert (2001) 
reported 0.4% as a mean concentration of organic carbon within the topsoil of an area in 
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Sudan zone and 0.2% for the Sahel zone and further opined that the soils in the Sudano- 
Sahelian zone are inherently low in OC. 

Soil organic matter OM (p= 0.01), were significantly different within the land uses, 
0.53g/kg was recorded in grass land, 0.55 g/kg in cultivated land and 0.41 g/kg in forested 
land. All were found to be very low within the land use types which could be attributed to low 
vegetative cover as a result of low rainfall that characterizes Sahel savanna areas like Gashua. 
The study area is characterized by low precipitation and high temperature, this may also 
contribute to low level of OM as reported to have effect on low of organic matter 
accumulation (Hamza and Anderson, 2010). Similarly the low organic matter content of the 
soils in all areas could be attributed to factors such as continuous cultivation, frequent burning 
of farm residues commonly carried out by farmers in the area which tends to destroy much of 
the organic materials that could have been added to the soil (Sharu et al., 2013). 

Available phosphorus (p=0.18) was found to be not significant within the land uses, 
8.75mg kg! was recorded in grass land, 8.31mg kg" in cultivated land and 9.34 mg kg" in 
forested area also generally rated low (<10 mg/kg) according to Esu (1991). This result falls 
within the range presented by Alhassan et al. (2020) that available P in the soils of Bade 
ranged from 7.31- 11.37 mg kg’. Total nitrogen (p= 0.02) is significantly different within the 
land use types. The mean value of 0.04 g/kg was recorded in grass land, 0.08 g/kg in 
cultivated land and 0.07 g/kg in forested area, indicating overall low level of N according Esu 
(1991) ratings. The low TN observed in the study area might be due to continues cultivation 
by farmers over the years and crop removal by uptake. 

The CEC (p=0.00) was significantly different within the land uses. The mean value of 
7.32 cmol(+)/kg was recorded in grass land, 8.11 cmol(+)/kg in cultivated land and 7.74 
cmol(+)/kg recorded under forested area. The CEC content is relatively medium as shown that 
the CEC content between 6 - 12 cmol(+)/kg is median reported by Esu (1991). As the soils 
were generally low in organic carbon contents, it is probable that the medium CEC obtained 
could be as a result of the presence of high activity clays and high content of basic cations 
from parent materials. 

The Ca** content (p=0.01) was significantly different within the land uses. The mean 
value recorded in grass land was 6.83 cmol(+)/kg, in cultivated land 6.24 cmol(+)/kg and in 
forested area 6.11 cmol(+)/kg. The calcium concentration was found to be very high within 
the land uses. The Mg™ contents (p=0.01) were found to be significantly different within the 
land use 1.11 cmol(+)/kg was recorded in grass land, 1.68 cmol(+)/kg in cultivated land and 
forested area 1.44 cmol(+)/kg. The high level of Mg** might be due to the dominance of Mg 
bearing minerals in the area. The K* concentration (p=0.00) was significantly different within 
the land use types. Moderately K* content was recorded in grass land (0.24 cmol(+)/kg) and 
forested area (0.22 cmol(+)/kg) and significantly higher concentration of K* was recorded in 
cultivated land with the mean value 1.13cmol(+)/kg. 

The Na’* content of the soil (p=0.87) was not significantly different within the land use. 
Na* recorded within grass land has a mean value 0.15 cmol(+)/kg), while in the cultivated 
land and forested area were 0.14 cmol(+)/kg each. High soil Na* (>0.3) was also reported by 
Amalu (1997). The electric conductivity (EC) (p=0.42) is not significant different within the 
land uses. The mean value recorded in grass land was 0.15 (ds/m), 0.10 (ds/m) in cultivated 
land and 0.15 (ds/m) in forested area. 
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Table 2. Effect of Land use on Soil Chemical Properties 


Land use pH EC OC OM CEC TN Mg" Ca** K* Na™ AP 
(ds/m) g/kg g/kg cmol(+)/kg g/kg emol(+)/kg cmol(+)/kg cmol(+)/kg cmol(+)/kg mg/kg 

GL 6.98 0.15 0.9 0.53a 7.32c 0.04b 1.11c 6.83a 0.24b 0.15 8.75 
CL 6.76 0.1 0.89 0.41b 8.1la 0.08a 1.44b 6.11b 1.13a 0.14 8.31 
FR 7.1 0.14 0.93 0.55a 7.74b 0.07a 1.68a 6.24b 0.22b 0.14 9.34 

SE 0.17 0.02 0.11 0.03 0.05 0.01 0.09 0.18 0.03 0.02 0.51 
Pr(>f) 0.19 0.42 0.92 0.01 0 0.02 0.01 0.01 0 0.87 0.18 


CV% 3.55 32 1711 9.79 0.99 26.69 __‘9.87 4.08 10.61___2.05 8.17 
GL=Grass Land/ CL= Cultivated Land/ FR=Forested Area/ SE=Standard Error/ CV=Coefficient of Variance 
Effect of land uses on soil moisture content 
The mean soil moisture content in three land uses are shown in Table 3. The soil 
moisture content within the land uses were highly significantly different (p = 0.002). In 
general, higher mean moisture content (13.53%) was recorded in forested area (FR) which is 
at par with 13.01% recorded under cultivated land (CL), while the minimum (10.39%) was 
found in grassland (GL). The high moisture content recorded under forested area may be due 
to high organic matter content through shade provided by the trees canopies which may 
reduce evaporation and decayed leaves and grasses which improve the soil organic matter 
content and help in soil moisture retention. Also there was high coefficient of variability (CV 
= 43%) within the land uses as per soil moisture content. This agrees with Yao et al. (2012) 
who reported that land use was one of the main factors that affected the spatial variations in 
soil moisture. 
Table 3. Mean soil moisture content under Land use types 


Land use Soil Moisture Content (%) 
Grassland (GL) 10.39b 

Cultivated land (CL) 13.0la 

Forested area (FR) 13.53a 

SE 0.922 

Pr(>F) 0.002 

CV (%) 43.00 


Effect of sampling depth on soil moisture content (%) 

The soil moisture content within the sampled soil depths showed no significant 
difference (p=0.495) between the moisture content (Table 4). This is consistent with Gao and 
Shao (2012) who observed insignificant variations of soil water at 0-20 cm among various 
land uses they studied (cropland, shrub land, woodland, and grassland). 


Table 4. Effect of sampling depth on soil moisture content 


Sampling depth(cm) Mean 
0-15 12.05 
15-30 12.58 
SE 0.775 
Pr(>f) 0.495 
CV% 44.28 
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Effect of temporal variation on soil moisture content 

The effect of sampling time (dates) on soil moisture content was presented in Table 5. 
The result shows that high significant differences exist between the soil moisture content at 
the different sampling periods. The result showed high moisture content (16.77%) was 
recorded at the beginning (08/08/2021) and then continue to decline up to the end of August, 
2021 (9.52%). The moisture content then increased to 14.55% on 05/09/2021 and decline 
also to the end of experiment on 31/10/2021 with 5.16%. The sudden increase in soil moisture 
content after its decline was due to the heavy rainfall experienced on 03/09/2021 and decline 
to the end of the experiment as no much rainfall was recorded. The increase and decline in 
soil moisture shows that rainfall has greater influence in soil moisture content, which brings 
the need for moisture conservation in days without rain. 


Table 5. Effect of sampling date on soil moisture content 


Sampling date Mean SMC (%) 
08-08-21 16.77a 
14-08-21 16.42ab 
18-08-21 15.23ab 
22-08-21 12.49abed 
27-08-21 11.84bed 
31-08-21 9.52cde 
05-09-21 14.55ab 
19-09-21 13.06abc 
03-10-21 12.39abcd 
17-10-21 7.99de 
31-10-21 5.16e 
SE 1.45 
Pr(>f) 0.0001 
CV% 35.23 


Effect rainfall on soil moisture content 
There was high significant difference in soil moisture (p=0.004) measured after a 
rainfall event and when there was no rain (Table 6). 


Table 6. Effect of rainfall on soil moisture content 


Rainfall effect Mean SMC (%) 
No rainfall 11.31b 

After rainfall 13 .5la 

SD+ 5.4 

Pr(>f) 0.004 

CV,% 43.42 


Highest soil moisture content was observed (13.51%) after rainfall and (11.31%) when 
there was no rainfall. This was in agreement with the report that, in arid land, the soil moisture 
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content often reached its highest value after a heavier rain, and when the rain had stopped for 
a couple of days, the soil moisture content decreased sharply due to evapotranspiration, 
absorption by roots or runoff etc. (Fu et al., 2003) further opined that precipitation was the 
main factor that influenced the temporal patterns in soil moisture (Yao et al., 2003). 


Conclusions 

This research work was conducted to evaluate soil moisture content of different land 
use selected in Federal University Gashua, Bade Local Government Area of Yobe State. The 
land use types are: grass land, cultivated land and forested area. The study showed that the 
soils of the study area of Federal University Gashua Bade LGA area are Sandy loam, high in 
bulk density, low porosity, very low in organic matter content. pH was neutral, EC values, 
total Nitrogen, organic carbon were very low in the soils, while sodium, available P, Ca, Mg, 
K and CEC ranged between medium to high concentration. The Land uses have shown greater 
influence on soil moisture retention, which also vary when measured periodically and 
consistently also been influenced by rainfall. 
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Abstract 

Suitability evaluation of soils of Hyuku in Wukari Local Government Area of Taraba 
State, Nigeria was carried out with the view to characterize and evaluate the soils for its 
suitability for maize production. Two soil units were identified based on their physical, 
chemical and morphological features using grid method of soil survey. Two profile pits were 
sunk in each unit, making a total of four pits in all. The soils were deep (80 to 115 cm), well 
drained with sandy loam surface textures and sandy loam to clay loam subsurface textured. 
The soil reaction was moderately acidic (5.65 to 6.75) with high base saturation (95.07 to 
96.40 %), low organic carbon (0.13 to 0.34 %), low nitrogen (0.17 to 0.67 %), and moderate 
to low available phosphorus (2.77 to 5.52 ppm). The soil characteristics were matched with 
the requirement of maize and the overall suitability rating of the soils were obtained using 
limitation method. Units I and II soils were rated as moderately suitable (S2) for maize 
production. 


Keywords: land suitability, evaluation, maize production and soil 


Introduction 

Maize is one of the main staple crops grown in Nigeria. Improving yield and 
efficiency in resource use is vital to ensuring adequate food security and environmental 
sustainability. The slight increase in maize production in Nigeria has been attributed to 
increase in hectares devoted to its cultivation and not optimum yield of each farmland. Recent 
report of FAOSTAT (2012) has also shown that the cultivation of areas used for maize 
production (Kaduna state) have declined by about 20 % from the periods of 1991 to 2000 and 
2001 to 2010. In addition, the potential yield of maize crop using the improved seedling in 
Nigeria is not enough compared to that of the developed countries (Oyekale, 2008). Also, 
very key is the improper land-use resulting to land degradation and decline in agricultural 
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productivity generally with maize inclusive. One of the best ways of soil conservation is the 
allocation of the soil to its most suitable use. Allocation of the soil/land to its best uses is 
referred to as Land use planning in which the characteristics of the land are compared or 
matched with the requirements of particular land use. Land use planning also known as land 
suitability ensured optimum or even maximum output and at the same time guarantees 
minimum damage to the land resources. 

Land suitability is the fitness of a given type of land for a defined use. Land suitability 
assessment for agriculture is meant to evaluate the ability of a piece of land to provide optimal 
ecological requirement of a certain crop variety. That is, assessing the capability of land is 
enabling optimum crop development and maximum productivity. Thus, evaluation needs a 
specification of the respective crop requirement and calibrating them with the terrain and soil 
parameters (Dent and Young, 1981). The identified limiting factors could be managed to suit 
various crop requirement and improve crop productivity. This is a pre-requisite to 
productivity maximization in agricultural sector (Pareta and Jain, 2010). Land evaluation 
provides avenues for sustainable land use since land will be used according to its capability. 
This therefore, makes it mandatory in this study to carry out land suitability in order to ensure 
that the selected site is suitable and capable of sustaining long term production of maize. 


Materials and Methods 

Study Area 

The study area is located in Hyuku, Wukari Local Government Area of Taraba State, 
North East of Nigeria. It lies between latitude 7°51' to 7° 85'N and longitude 9° 46' to 9° 78' E, 
of the Greenwich meridian. The entire area is a gently undulating plain, with a mean altitude 
of 200m above the Sea level. The drainage systems drain northward and serve as tributaries 
network to the river Benue while Eastward discharge of rivulets and other smaller tributaries 
from Wukari town drains towards the Donga River, which is a major tributary of River 
Benue. The mean annual rainfall value ranges from 1000-1500mm. The onset of the raining 
season is usually around April, while the offset period is October. The mean maximum 
temperature is being experienced around April at about 40°C while the mean minimum 
temperature occurs between the period of December and February at about 20°C (NIMET, 
2015). The study area is situated over Cretaceous sandstone popularly known as Bima-sand- 
stone. The vegetation consists of grassland interspersed with trees and shrubs. The soils are 
moderately deep to deep, well-drained to poorly drained, with loamy sand to sandy loam 
surface over sandy clay loam to sandy clay subsurface (Peng et al.,2006). Land use includes 
the cultivation of Maize (Zea mays L), Cowpea (Vigna unguiculata L. Walp), Groundnut 
(Arachis hypogea), Rice (Oryza sativa), Guinea corn (Sorghum bicolor), Beniseed (Sesamum 
indicum), Bambara nut (Voadenzia subterranea), Yam (Dioscorea spp), Okra (Albemuschus 
esculentus), as well as Cashew (Anacardium occidentale) and Mango (Mangifera indica) as 
well as Oil palm (Elaeis guineensis) plantations. Grassland areas support animal grazing 
activities. 
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FIGURE 1: MAP OF WUKARI LGA SHOWING LOCATION OF STUDY AREA, 


Field and Laboratory Studies 

The area was soil surveyed using the grid method. Auger point investigations were 
carried out at 100 m intervals along traverses cut at 100 m apart on the baseline. Based on 
these investigations, two soil units were identified and two profile pits were sunk in each. The 
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pits were described according to the guideline for soil profile description (SSS, 2014) and the 
samples were collected and taken to laboratory for physical and chemical analysis. The air 
dried, crushed and sieved (d<2mm) samples were analyzed for particle size distribution 
(PSD), pH, organic carbon, CEC, exchangeable bases (Ca, Mg, K and Na), total N, available 
P, and base saturation. 

PSD was determined by Bouyoucos hydrometer method (Day, 1965). Soil pH was 
determined by electrometric methods as described by IITA (2015). Walkley and Black 
method as described by Nelson and Sommers (1982) was employed for organic matter content 
determination. Total N was determined using the modified macro-kjeldahl method as 
described by HTA (2015). Bray No.1 method was used for extractable P. For exchangeable 
bases, Ca and Mg were determined by AAS while K and Na were done by flame photometer. 
CEC was determined by IITA (2015) procedures while base saturation was calculated by 
dividing TEB by ECEC and multiplying by 100%. 


Land Suitability Evaluation (LSE) 

The FAO framework for soil suitability classification was used for the study, 
according to Sys (1985). The suitability of the soils was evaluated by a direct method where 
the established soil land characteristics (Sys, 1985) were cross-matched with the soil land 
characteristics obtained from the study. The outcome was cross-matched with the suitability 
class rates and agricultural uses, according to Ezeaku (2011). 


Results and Discussion 

Soil Morphological Characteristics 

The soil morphological characteristics examined in the study area include; the soil 
colour, texture, structure, depth and consistence (Table 1). The soils were predominantly dark 
reddish brown (2.5 YR 2.5/6), reddish brown (2.5 YR 5/4) and light reddish brown (2.5 
YR6/3) colours (moist) in their A and Ap horizons of all the profiles. This could be attributed 
to the presence of relatively high organic matter which is the main colouring agent on surface 
soils (Ufot, 2012). In the subsurface horizons of the profiles, the presence of light reddish 
brown (2.5 Y 6/4 moist) and reddish colour is an indication of complete organic matter 
decomposition (oxygen rich soil) and seasonal fluctuation of ground water table resulting in 
the reddish yellow (SYR 7/6) mottles observed. Such mottles were mainly confined to root 
channels because of the oxidizing effect of oxygen, which enters the channels as soon as the 
water table receded. The soils were deep ranging from 80 to 110 cm. They are generally well 
drained with sandy loam surface textures while the sub-surfaces varied from sandy loam to 
loam to silty clay loam to loam. The textures of these soils reflected the parent rocks from 
which they are formed, the rate and the nature of some weathering processes (Ahukaemere et 
al., 2016; Idoga and Azagaku, 2005). The soil structures were well developed being moderate 
medium crumb to moderate medium subangular blocky at the surface while the subsurface 
was moderate medium subangular blocky to strong medium crumb. The structures are 
functions of organic matter content of the soil (Akinyemi and Vivian, 2001). 
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Physical and Chemical Properties of the Study Area 

The physical properties of the soils of the study site are shown in Tables 2. For all the 
profiles, sand fraction ranged from 31.2 to 75.4% in surface and subsurface soils respectively. 
Similarly, silt fraction varied from 20.4 to 38.8% respectively; and a range of 4.0 to 38.0% for 
clay fraction at the surface and subsurface respectively were observed. The texture of the soils 
studied was sandy loam at the surface while subsurface horizons varied from sandy loam to 
clay loam. 

The chemical properties of the soils are shown in Table 2. The overall results show 
that the soils were moderately acidic with pH ranging from 5.65 to 6.75. Organic carbon 
varied from 0.04 to 0.34 while that of total nitrogen varied from 0.17 to 0.75%. Available 
phosphorus content ranged from 2.65 to 5.52 mg kg’'. The organic carbon together with total 
nitrogen and the phosphorus were rated low. Generally, these may be attributed to release 
from plant tissues, gaseous loss, surface runoffs, leaching, climatic factors, vegetation, human 
activities and initial soil pH. Loss of N through denitrification and volatilization may also 
contribute to the low level of N in the wetland (Brady and Ray, 2014). The exchangeable 
bases were low as a result of the nature of the underlying parent materials, intensity of 
weathering, leaching, low activity clay, low O.M, erosion and lateral translocation of bases. 
Ca was the most dominant cations with values ranging between 4.42 and 4.6 cmolkg” in the 
exchange complex. This may be linked to the occurrence of exchange sites which have 
specific affinity to Ca (Idoga, 1985) or may be due to the fact that Ca is least easily lost from 
exchange site or has high displacement ability over other cation exchange reaction. The Mg 
values varied between 0.66 and 0.93 cmolKg"! while that of K and Na varied from 0.24 to 
0.31cmolKg” 0.76 to 0.94cmolkg"' respectively. These values confirmed the predominance of 
Ca over Mg, K and Na as observed by Idoga, 1985, and Ogunkunle, 1989. CEC content 
varied between 6.62 and 6.80 cmolkg” and were rated low. The CEC of the soils were low to 
medium with values ranging between 3.75 and 8.34cmolkg'. The low CEC values may be 
attributed to the nature of clay minerals (kaolinite) and low organic carbon level of the soils. 
The percentage base saturation values were very high (95.07% to 96.69%). The high base 
saturation could probably have associated with the presence of weathered minerals which 
release nutrients into the soil and their alluvial nature (Table 2). 


Suitability Rating for Maize Production 

Land suitability assessment for agriculture is meant to evaluate the ability of a piece of 
land to provide optimal ecological requirement of a certain crop variety. That is, assessing the 
capability of land is enabling optimum crop development and maximum productivity. Thus, 
evaluation needs a specification of the respective crop requirement and calibrating them with 
the terrain and soil parameters (Dent and Young, 1981). The identified limiting factors could 
be managed to suit various crop requirement and improve crop productivity. This is a pre- 
requisite to productivity maximization in agricultural sector (Asadu and Ezike, 2017). Land 
evaluation provides avenues for sustainable land use since land will be used according to its 
capability. This therefore, makes it mandatory to carry out land suitability in order to ensure 
that the selected site is suitable and capable of sustaining long term production of crops with 
maize inclusive. 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 73 


Gani et al., 2024 


sp 
Sp 


Arepunog 


S]OO1 OUT MO] 
$]OO1I OUT MO] 
S]OOI OUT MO] 

$]OO1 OUTJ UOWWOD 

$}001 wuntIpaw Aue 


$]00Y UL] MOT 
$]00Y JUL] Mo] 
S]OOI OUT MO] 
$]OO1 SULT MOT 

s]OOY WOWUOD UINIpaf 


$]001 WINTpoul M9,J 
$]OO1 UINIpet UOUTWIOD, 

$]001 asIv09 AURIAL 

$]OOI asIv09 AUR] 


WDQ8 IB YIOI JUSIBYOO 
piey/s}ooy sul Moy 
$]00Y JUL] Mo7 
$]00Y JUL] MO] 
s}OOY WINTPI|A] May 


uoIsn[ouy 


ASS 
ASS 
ASS 
ASS 
ASS 


ASS 
ASS 
ASS 
ASS 
ASS 


MSA 
ASS 
ASS 
ASS 


MSA 
ASS 
ASS 
ASS 


goua}sIsuod, sIMONNG sNIXI], 


MaSWC 
MASI 
MaSWC 
MaSWC 
MaSWC 


MaSWC 
MaSWNC 
MaSWC 
MaSWC 
MASc 


Woe 
MASc 
MaSWC 
MASc 


MAS 

MASc 

MASc 
WW 


IS 
IS 
IS 
IS 
IS 


IS 
IS 
IS 
IS 
IS 


TD 
IS 
IS 
IS 


TD 
IS 
IS 
IS 


dcW 


Te1p 
TOW 


9/L AAS 


Soop 


poy = ¥ ‘UMorg 1481] 
= AT ‘UMOIg YsIppoy SIT = {YT ‘UMoIg YSIMOT[E A Ye ={AC ‘UMOIg YsIppoy = Cy ‘UMoIg YsIppoy eC = CAC “Mo[[2A YsIppoy =AY “UMorg ysikeip =o ‘umorg yeq=_qq :nojoy 
YWOOUSs [enpeis = $3 ‘YJOOUIS dsNJIp = sp :AIBpunog 
Jom onseyd-uoN = MdNy 10M AyONS-UON = MSN ‘10M AOS = MS OM Ayons ATOA = Md A ‘10M Axons AIOA = MSA ‘30M AXONS ANYSITS = ASS :2dUa}SISUOD 
Ayo] Jepnsuvgns wumnipour swesopoyl = YASINTAUINI winipoul seIOpoy[ = YOWT ‘Quinso asreos SuoNS = YOOE soanjon.ys 

weol Ae[D =JO ‘weol Apues =7J¢g :34n4xa 7, 
juoutuoid asie0o AurW=deW :Spleyeqd SulpWO; 


dav/SaAs'c SI1-06 
aD7/SAAS'T 06-09 
CAQ9/ra AO! 09-8¢ 
AaAI/IAAS 8e-0C 
TIV/9AAG L 0c-0 
°C MHO1d 7 Huy) 

aw l8/9 YAS’ OI T-68 
ad 19/L YAS 68-SS 
aw Ic/9 YAS SS-LE 
av c/o YAS LE-CT 
aqadde/e AAS cC-0 
T 901d 7 Huy) 
adD7/SAAS'T 06-0” 
TIV/MAAS L OV-Ve 
AAI/IAAS ve-Ic 
CAQ9I/raAAOl 17-0 
°C OJOId T uly) 

U7, MAST O8-SV 
auT"/, YAST Sv-Of 
du /. WAST Oc-81 
aud’/..-aAS 7 81-0 
*T afJOrld T yu) 

(ui9) 

({SIOUW) IOTOD widoq 


by q 
44 
4q 

av 
dy 


by q 
44 
4q 

av 
dy 
44 


av 


OZLIOH 


nynky fo sjiog ays fo saijsadosg [va1sojoydiopy’{T IQqvy, 


74 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 


Gani et al., 2024 


69°96 
CcS°S6 
81°96 
€8°S6 
S9°S6 


096 
99°S6 
C€8°S6 
Iv'S6 
e1'S6 


Iv'v6 
ve c6 
6r' 76 
LS'S6 


ST'v6 
VC C6 
L0°S6 
6S°S6 


(%) 
Sd 


g9'9 
OL'9 
08°9 
cL'9 
L9°9 


99°9 
89°9 
OL'9 
9L'9 
8L'9 


c9'9 
9L'9 
IL9 
LL‘9 


L9°9 
08°9 
OL'9 
08°9 


cc'0 
ve0 
08°0 
STO 
810 


670 
Sc0 
LY'0 
vI0 
910 


ce0 
eC 
LY0 
910 


ce'0 
cc 
LV'0 
S10 


ev9 
ovo 
vs9 
vv9 
829 


cv9 
669 
cv'9 
sV9 
cV9 


Sc9 
Te9 
ve9 
LvV9 


8c9 
ve9 
LEO 
os'9 


6L'0 
08°0 
08°0 
6L'0 
8L°0 


9L'0 
08°0 
6L'0 
130 
c8'0 


vL0 
9L'0 
c8°0 
160 


LL‘0 
8L°0 
30 
£60 


wumntsousep = S| ‘WINnToTeD = BD SWINIpOS = BN ‘uINIssejog = y ‘snioydsoyd s[quireAy = dV ‘UdsONIN = N 


vco0 
970 
L70 
6c'0 
Sc0 


90 
8c'0 
870 
8c0 
80 


9¢'0 
Sc0 
L70 
00 


970 
90 
8c0 
Te0 


88°0 
8L°0 
c8°0 
cL'‘0 
£60 


6L'0 
LL‘0 
130 
08°0 
180 


L9°0 
OL'0 
69°0 
L9°0 


L9°0 
69°0 
89'0 
99°0 


CoV 
9S'V 
S9o'T 
vov 
cy vy 


19'y 
LOY 
vov 
9ST 
VST 


897 
097 
9ST 
6ST 


897 
19°V 
LOV 
09°77 


> (,.8yjowo) <—_______ 
va dal 


JHOA 


BN 


M 


3 


78) 


O8'V 
88'e 
OL'c 
vrs 
067 


8L'V 
LOE 
18°C 
ors 
COV 


09'7 
Or 
cyy 
S9'T 


09°€ 
cls 
css 
LLT@ 


(,-34sur) 
dV 


870 
L70 
ce0 
890 
Sv A0) 


60 
LY'0 
be0 
9¢°0 
SVv'0 


€S'0 
SL‘0 
cS°0 
cS'0 


Sv'0 
L9'0 
Sv'0 
vv0 


(%) 
NI 


e100 
vo0 
LT0 
veo 
ce0 


00 
veo 
STO 
90 
IT0 


610 
cT'0 
970 
ve 0 


80°0 
ITO 
S10 
€c0 


(%) 
JO 


nynky fo spiog fo saijsadosg JDIMAYD puv JVIISKYd °Z Iqvy, 


co's 
vL’S 
08'S 
8S 
pos 


OLS 
go's 
OLS 
SL’S 
S8'¢ 


eL’s 
99'S 
Lv9 
9¢°9 


c6'S 
o8'¢ 
99°9 
SLO 


(0°H) 
Hd 


‘IS 
IS 
‘IS 
IS 
IS 


IS 
IS 
‘IS 
‘IS 
‘IS 


YO 
IS 
IS 
IS 


ste) 
‘IS 
IS 
IS 


vIt TOE 8r S11-06 
OLI Sle SIs 06-09 
Sel Ore Tze; 09-8€ 
TOI «OEE == «89S 8E-07 
rs tre  L09 07-0 
27 PYoig 7 Wg 
O6l +v0e~ 909 OII-68 
Ost s9t 9¢¢ 68-SS 
Ort 88E TLr CS-LE 
OOl sve TLS LE-C 
09 Sre 76S 77-0 
27 org 7 Huy, 
Ore SET TTS 06-0r 
OTL Tw gs9 Orr 
T6020 OO ST—C«S'S'Qd VE-I7Z 
gE g0t SL 17-0 
:7 aod [Wu 
ose S80E TIE 08-Sr 
OOl gst TZI9 Sr-0€ 
O08 88% TED O€-81 
Ov St TEL 81-0 
27 WYoig T WuQ 
—> (%) +«— (WD) 
WIS pues yydaq 


ainyxay, Avg 


1q 
Qq 
4q 


uO0ZLIOH 


75 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 


Gani et al., 2024 


The suitability rating of Hyuku soils were carried out by comparing the qualities of 
the soil with the requirements of maize. The chemical characteristics of the soils such as pH, 
OC, CEC, exchangeable bases, available P and total N (Table 2) were found to be either 
conducive to maize production or can be amended by individual farmers and therefore cannot 
be considered to be permanent limitations. 


Table 3. Land Use Requirement for Maize 


Land Quality and Characteristics 100-95 94-85 84-40 39-20 19-0 
Sl S2 S3 Nl N2 
Climate (c): 
Annual rainfall (mm) 850-1250 850-750 750-600 600-500 - 
1250-1600 1600- >1800 
Length of dry season (days) Mean annual 150-220 130-150 1800 90-110 <90 
max temp. (°C) 22-26 22-18 110-130 36-30 $36 
18-16 
Relative humidity (%) 42-50 50-80 >32 3 = 
>80 
Topography (t): 
Slope (%) 0-2 2-4 4-8 8-16 
0-4 4-8 8-16 16-30 >30 
Wetness (w): 
Flooding FO MO Fl Aeric Poor 
Drainage Good Moderate Poor Poor Drainable 


Soil physical Characteristics (s): 


Texture/Structure CL, L SL, LS LCS CS,S S 
Coarse fragments (%) 0-10 cm <3 3-15 15-35 35-55 

Fertility (f): 

CEC (cmolkg | clay) >24 16-24 <16(-) 16(+) - 
Base saturation (%) >50 35-50 20-35 <20 - 
pH 5.5-7.0 5.5-7.0 5.0-8.0 5.0-8.0 

OC (%) 0-15cm >2 1.2-2 0.8-1.2 <0.8 - 
Av. P (mgkg 1) >22 13-22 7-13 3-7 <3 
Total N (%) >0.15 0.10-0.15 0.08-0.10 0.04-0.08 <0.08 
Extr. K (cmolgk 1) >0.5 0.3-0.5 0.2-0.3 0.1-0.2 <0.1 


Key: FO=No Flooding; Fl=Seasonal Flooding; CL=Clay Loam; SL=Sandy Loam; LS=Loamy Sand; 
SCL=Sandy Clay Loam; S=Sand. 
Source: Modified from Sys (1985) 


Soil depth, drainage, slope, porosity, texture, structure and soil fertility status are 
important physical and chemical characteristics that influence maize growth, development 
and yield. The climatic condition of the area (temperature and rainfall) were highly suitable 
(S1) and favorable for maize production as values were higher than 640 mm and little above 
28 °C (FAO, 2007). The textural class was found to be highly suitable (S1) for maize 
cultivation across the soil units. The soils were well drained, free from flooding and were 
therefore highly suitable (S1) for maize production. The organic matter content and CEC 
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were low and so rated marginally suitable (S3) for maize cultivation. Total nitrogen was low 
to moderately high and therefore rated moderately suitable (S2) for maize. Generally, the 
aggregate suitability class of the soils suggested that they were moderately suitable (S2) for 
maize production (Table 3 & 4). 


Table 4. Suitability Class Scores of the Profiles in the Study Area for Maize 


Land Quality and Characteristics Profile Profile Profile Profile 
I II Ill IV 

Climate (c): 
Annual rainfall (mm) $195) S$1(95) S$1(95) $195) 
Length of dry season (days) $1(100) 
Mean annual max temp. CC} $1(100) $1(100) $3(40) $1(100) 
Relative humidity (%) pa) 83(40) $3(40) 
Topography (t): 
Slope (%) $1(100) S$1Q5) $2(85) $2(85) 
Wetness (w): 
Flooding 
Drainage S$1(100) S1Q5) $2(85) S3(40) 
Soil physical Characteristics (s): 
Texture/Structure 

$1(95) $1(95) $2(85) $2(85) 
Fertility (f): 
ECEC (cmolkg | clay) $3(40) S3(40) $3(40) $3(39) 
Base saturation (%) $1(100) $1(100) $1(100) $1(100) 
pH $2(85) $2(85) $2(85) $2(85) 
OC (%) 0-15cm N1(20) N1(20) N1(20) N1(20) 
Av. P (mgkg 1) N1(20) N1(20) N1(20) N1(20) 
Total N (%) S$1Q95) S$1Q5) N1(20) N1(20) 
Extr. K (cmolgk 1) S$3(84) S3(40) S$3(84) S3(84) 
Mean value 74 63 66 63 
Aggregate suitability class S2 S2 82 S2 
Limiting Characteristics d,n k,n n dn 


Aggregate suitability class scores 
100-75=S1, 74-50=S2, 49-25=S3, 240=NI1 
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Figure 2. Interaction of Particle Size Distribution (PSD) of the Soils in Unit 1 
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Figure 3. Interaction of Particle Size Distribution (PSD) of the Soils in Unit 2 
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Figure 4. Interaction of pH, OC, TN and Av. P of the Soils in Unit 1 
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Figure 5. Interaction of pH, OC, TN and Av. P of the Soils in Unit 2 
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Figure 6. Interaction of Exchangeable Cations of the Soils in Unit 1 
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Figure 7. Interaction of Exchangeable Cations of the Soils in Unit 2 
Conclusion 


Selected soils of Hyuku were studied with the aim to characterize and evaluate the 
soils for its suitability for maize production. Two major soil units were identified in the area. 
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The soils were deep, well drained and had sandy loam to clay loam texture. They were 
generally moderately acidic in reaction with pH values ranged between 5.65 and 6.75. This 
pH range is safe for any sustainable maize production. They also have high base saturation, 
low value of available P, organic carbon, total nitrogen, and exchangeable cations meaning 
that they have high potentials to respond positively to any addition/application of fertilizer be 
it organic or inorganic. Based on the physical and chemical properties of the soils, they were 
rated as moderately suitable (S2) for maize production. 


References 

Ahukaemere, C. M., Ndukwu, B. N., Ihem, E. F. and Onwuso, G. I. (2016). Profile 
characteristics of selected rice producing soils of Abia State in South-Eastern Nigeria. 
Nigerian Journal of Soil Science, 26: 146-153. 

Akinyemi, O and Vivian A. (2001). Influence of combined application of Cowdung 
and Inorganic Nitrogen on Microbide Respiration and nitrogen transformation in an Affisols. 
Nigeria journal of soil science, 2:15-20. 

Brady, N.C and Ray W. R. (2014). The Nature and Properties of soil. 14™ Edition 
Prentice Hall, New Education Inc. India 621 pp 

Day. P.R. (1965). Particle fraction and particle size analysis. In: C.A. Black(ed), 
Methods of Soil Analysis, ASA, Madison, Wisconsin, USA, Agron.9: 547-577 

Dent, D. and Young, A. (1981). Soil survey and Land Evaluation. George Allen and 
Unwin Pub. Ltd, London, 278pp. 

Ezeaku, P.I. (2011). Methodologies for Agricultural Land-use Planning - sustainable 
soil management and productivity. Great AP Express Publishers Ltd, Nsukka, Nigeria. 69p. 

FAO, 2007. Land evaluation, towards a revised framework. Land and Water 
Discussion Paper 6. 

FAOSTAT, Food and Agriculture Organization/Statistics. (2012). Trade data base, 
Production data base http://faostat.fao.org. Accessed 01/02/2013. 

Idoga, S. (1985). The Relationship of Topography to Soils and Vegetation in the 
Northern Guinea Ecosystem, Nigeria. MSc thesis Submitted to the Department of Soil 
Science, Ahmadu Bello University, Zaria, Nigeria. 182p. 

Idoga, S. and Azagaku, E. D. (2005). Characterization and classification of soils of 
Janta area of Plateau State, Nigeria Journal. of Soil Science. 15:116 — 122. 

IITA, (2015). Selected Methods for Soil and Plant Analysis. Manual Series NO.13, 
IITA Ibadan, Nigeria 

Nelson, D.W. and Sommers, L.E. (1982). Total carbon, organic carbon and organic 
matter. In: Methods of Soil Analysis part 2. (Eds A. I., Miler R.H. and Kenney D.R.) 
Madison, W.I., Amer. Soc. Agron. 42:539 — 379. 

Ogunkunle, A.O. (2004). Soil survey and sustainable land management. Proceedings 
of the 29th Annual Conference of the Soil Science Society of Nigeria, 2004, Dec. 6-10. 
University of Agriculture, Abeokuta, Nigeria, 12-16. 

Ogunkunle, A.O. (1989). Variation within a taxonomic unit; The Ibadan series 
example in South-Western Nigeria. Nigerian Journal of Soil Science, 9, 9-20. 

Oyekale, A.S. (2008). Land degradation, soil conservation practices and poverty 
incidences in southwestern Nigeria”. Agricultural Journal, 3(6),482-487. 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 81 


Gani et al., 2024 


Pareta K and Jain C. K. (2010). Land Suitability Analysis for Agricultural Crops 
using Multi-Criteria Decision Making and GIS Approach (2010) 
http://agropedia.iitk.ac.in/sites/efault/files/Land%20Suitability %20Analysis%20for%20A gric 
ultural%20Crops.pdf 

Nigeria Meteorological Agency (NIMET), (2015). Climate Weather and Water 
Information for Sustainable Development and Safety. 

Peng, K. Li, Luo, C. Shen, Z. (2006). Vegetation composition and heavy metal 
uptake by wild plants at three contaminated sites in Xiangxi area, China, Journal of 
Environmental Science and Health, Part A, 40, 65-76. 

Soil Survey Staff /USDA (2014). Soil Taxonomy. A basic system of soil classification 
for making and interpreting soil surveys (14"ed). Natural Resources conservation, service, 
USDA, Washington DC. USA, 869p. 

Sys C. (1985). Land Evaluation. State University of Ghent, Belgium. 

Ufot, U. O. (2012). Soils and the Environment for Colleges and Universities, Shadow 
Publishers Ltd, Owerri, Nigeria 1634p. 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 82 


Primary Nutrients Status of Three 
Vermicomposts Produced in a Nigerian Sahel 
Region 


Mahmud Sani, Alhassan Ibraheem'™, 
Sufiyanu Sani”, Aliyu Abdulkadir2, 
Ibrahim Hassan*, Zara Musa“ ,» Musab Ibrahim*” 


I - Department of Agronomy, Federal University Gashua, 
Yobe State, Nigeria 
2 - Department of Soil Science, Federal University Dutsin- 
Ma, Katsina State, Nigeria 
3 - School of Vocational and Technical Education, Sa’adatu 
Rimi College of Education, Kumbotso, Kano state 
4 - Department of Forestry and Wildlife Management, 
Federal University Dutsin-Ma, 
5 - Department of Agricultural Extension and Rural 
Development, Federal University Dutsin-Ma, Katsina State, 
Nigeria 


Corresponding Author: Sufiyanu Sani e-mail:ssani2 @ fudutsinma.edu.ng 
Received: 3 November 2023 Accepted: 10 April 2024 


Abstract 

A pot experiment was conducted at the soil and water laboratory of Federal University 
Gashua Nigeria, to assess the reaction and primary nutrients — NPK, and Organic Carbon 
(OC) status of three white grub-produced Vermicompost using different ingredients 
(Municipal Waste — Mw, Typha Grass — Tp, and Municipal Waste + Typha Grass — MyT,). 
Laid on a Completely Randomized Design (CRD), the number of live worms, worms' weight, 
Vermicompost weight, and worms' mortality percentage as production indices were recorded 
and statistically analyzed upon Vermicompost maturity. Matured Vermicompost Reaction (pH 
and EC) and Primary nutrients (NPK), and OC were determined in another laboratory upon 
taking the respective compost samples as Soil health Indices. The result of the analysis 
showed that Mw has the highest Nitrogen (%) = 1.60, the highest available Phosphorus = 0.47 
mgKg", and the highest Exchangeable Potassium = 720 mgKg", it also presented a mildly 
acidic reaction (pH=6.8), with the highest percentage of OC (2.80%) and the lowest basic 
salts EC = 0.35 dS/m. In terms of earlier determined live worm percentage, live worm weight, 
Vermicompost weight, and worm mortality rate, the highest findings are Mw = 64.5 %, Mw = 
40 grams, Mw = 301 grams, and Tp = 66.9 % respectively. In conclusion, Mw Vermicompost 
proves to be the best, taking into consideration its better indices. As the premier of its kind in 
the study area, among other recommendations, we encourage setting this study as a pilot, and 
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the conduct of more studies on the efficiency of Vermicompost on soil health and crop 
production in the study area. 


Keywords: vermicompost, vermicompost production, primary nutrients, sahel 
savannah, soil health 


Introduction 

Vermicomposting is the decomposition of organic residues in an aerobic environment 
by exploiting the optimum biological activity of earthworms and microorganisms (Garg and 
Gupta, 2009). Studies on the production of important vegetable crops like tomatoes 
(Lycopersicon esculentum), and eggplant (Solanum melongena) using Vermicompost have 
yielded very good results in Asia (Adhikary, 2012; Ashiwar & Hussain, 2015). Similarly, the 
overall productivity of many classes of crops was significantly boosted by Vermicompost 
applied at reasonable amounts per unit area (Adhikary, 2012). Positive findings after the use 
of Vermicompost as a source of organic manure in supplementing and/or substituting 
chemical fertilizers lead to it becoming popular among farmers around the world (Joshi et al., 
2015). Application of Vermicompost was reported in several literature to enhance most 
inherent biogeochemical aspects of crop production (Sinha et al., 2009; Ashiwar & Hussain, 
2015; Dhanalakshmi et al., 2014; Joshi et al., 2015), as well as its economy (Lim et al., 2015). 
A Study also suggested that Vermicompost enhances the formation of humic acid and plant 
growth-promoting bacteria, which is a plus for soil health and sustainable agriculture 
especially in the tropics (Laufer and Tomlinson, 2013). 

Vermicast, a cocoon compound produced during Vermicomposting, produced a higher 
percentage of the protein and carbohydrate contents as compared to chemical fertilizer in the 
study of Adhikary (2012).Vermicast was established to be the active constituent, and its 
action was believed to produced higher garden pea pods per plant, higher green grain weight 
per plant, and higher green pod yield per hectare as compared to chemical fertilizer, and the 
overall vegetable production (Adhikary, 2012; Ahirwar & Hussain,. 2015). 

Recently, Organic farming inputs have been adopted by farmers to boost soil fertility, 
OC content and the general health of arable soils (Lim et al., 2015). However, the knowledge 
of deliberate Vermicomposting is among the limiting techniques in agricultural production to 
the farmers, especially in the Sahel Region of Nigeria. Also, evident is the fact that, 
chemically synthesized inorganic fertilizers pose economic, health, and degradative risks in 
the long run (Adugna & Abegaz, 2015) . 

Located in the Sahel Savannah Agro-ecological zone, in the northern parts of Yobe 
State of North Eastern Nigeria, Gashua is a predominantly farming and fishing community. 
Despite its unfavorable climatic condition and low annual rainfall (500 to 1000mm) (FDALR, 
1990), Gashua soil is fertile and suitable for cultivation of cereal, tubers, fruits, and vegetable 
(Alhassan et al., 2018). The soil, however, was reported to have low OC content, a crucial site 
for ion exchange and buffering (Sani et al., 2019 ). 

In the Nigerian Sahel, with inorganic fertilizers being costly, sometimes scarce, and 
elsewhere regarded as a security threat, subsistent farmers find it economically impossible to 
apply enough for good growth (Alhassan et al., 2018). These Farmers depend largely on 
locally sourced organic fertilizers on one hand (Makinde et al., 2010; Sani et al., 2023). On 
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the other hand, a huge amount of waste in the form of Municipal Waste (Mw), Farm yard 
Manure (FYM), and numerous indigenous and exotic grasses that colonize farmlands and 
waterways such as Typha species, are generated and heaped on dump sites, posing potential 
environmental pollution and hazard (Agboola and Omueti, 1982). Incorporating these waste 
materials into the soil for crop production would serve as a source of organic farming inputs 
to build up an organic matter layer that is needed for a steady supply of nutrients in tropical 
soils (Agboola and Omueti, 1982; Noma and Sani, 2008), as well, a measure for pollution and 
bio-degradable waste control. 

The potential of Vermicomposting to turn municipal, waterways, and on-farm waste 
materials into an important farming input makes it an attractive proposition. Vermicomposting 
offers benefits such as enhanced soil fertility and soil health that engender increased 
agricultural productivity, improved soil biodiversity, reduced ecological risks, and a better 
environment (Adhikary, 2012; Sinha et al., 2009; Ashiwar & Hussain, 2015; Barik et al., 
2011; Dhanalakshmi et al., 2014; Dominquez & Edwars, 2004; Garg & Gupta, 2009; Joshi et 
al., 2015; Lim et al., 2015; Mistry, 2015; Sinha et al., 2010). However, many farmers, and 
especially those in developing countries and this study area, find themselves at a disadvantage 
as they fail to make the best use of organic recycling opportunities using earthworms. Thus, 
Vermicomposting could be one of the valuable options for Sahel region farmers to restore or 
enhance their agricultural soil's health. 

The major objective of this study is to know the reaction, primary nutrient and OC 
status of three Vermicomposts produced with different raw materials. Specific aims of the 
study include: Producing Vermicompost and comparing the production indices of different 
Vermicomposts 


Materials and methods 

This research was carried out at Federal University Gashua (FUGA). Gashua is a 
community located on Latitude: 12°52'26” N Longitude: 11°02’26" E Elevation above sea 
level: 339 m = 1112 ft in Yobe State, northeastern Nigeria. 


Experimental design 

The experiment was laid out on a Completely Randomized Design (CRD) where three 
vermicomposts are produced using three different ingredients (Fiber) in three times replica 
containers. The main raw materials are that make the three Vermicompost treatments are: 


1. Municipal Waste (Mw) 

2 Dried Typha Grass (Tp) 

a 50% Municipal Waste + 50 % Dried Typha Grass 
(MwTo) 


Vermicompost production 
Organic waste in the forms; Farmyard Manure (FYM), Municipal Refuse, and the 
White Grub worms were obtained from the Animal Science Teaching and Research Farm, 
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Faculty of Agriculture, Federal University, Gashua. Typha grass was obtained along Garbi 
River, a delta on the Hadejia-Nguru Wetlands, all in Yobe State, Nigeria. 

During preparation for the compost production, all organic wastes were sorted and air 
dried for a week. 11 white grub worms were introduced to perforated pots each for the three 
compost ingredients mentioned in 3.3 above, upon which an equal weight of the raw materials 
and farm yard manures was added. 100ml of distilled water was then sprinkled on the mix and 
the containers were covered. The same procedure was replicated three times for each 
treatment. Typha grass was shredded to increase its surface area and increase decomposition 
during rate. This setup was monitored and turned daily for three weeks and data was taken on 
the 21* day after set up. 

Sampling and determination of Vermicompost nutrient status 

A sample was collected from each of the matured Vermicomposts and used to 
determine Electrical Conductivity (EC), pH, Total Nitrogen (N), Available Phosphorus (P), 
Exchangeable Potassium (K) and Organic Carbon (OC) using standard laboratory procedures 
at the soil science laboratory, Bayero University Kano. 

EC and pH were determined with its respective meters using the method as described 
in (Mitchell and Soga, 2005). Total Nitrogen (N)was determined using the Kjeldhal digestion 
method in (Bouajila and Sanaa, 2011). Phosphorus (P) in Cmol/Kg was determined using 
Extractable Phosphate Bray and Kurtz —Method (Department of Sustainable Natural 
Resources, 1995), and converted to mg/Kg using the conversion factor Icmol/Kg 
=10,000mg/Kg. Exchangeable Potassium (K) was determined using the Ammonium 
saturation method described by (Burt, 2014). Organic Carbon (OC) was determined using wet 
oxidation using the method of Walkley and Black (1934). 

Method of data collection 

Data on mortality rate, number of live worms, weight of worms and compost weight 
were collected and recorded from Sample containers at compost maturity. 

White grub worms species Mortality rate: 

The mortality rate of white grub worms which were placed in the different prepared 
Vermicompost medium for the experiment was recorded by turning the sample with a spatula 
and watching with eye visibility count at compost maturity. 

Percentage of live White grub worms: 

Live worms in the respective pots after compost maturity were counted and the 
percentage was determined from the initial number at the start of the experiment. 

Weight of White grub worms: 

The weights of the worms in 3.5.2 above were determined using a digital weighing 
balance. 

Vermicompost weight: 

Matured Vermicompost of all pots was weighed after the removal of live and dead 
worms using a weighing balance. 

Data analysis 

The Vermicompost production indices data recorded from all pots was subjected to 
Analysis of Variance (ANOVA) using STAR 2.2 edition. 

L.S.D was employed for mean separation where significance existed at 0.5% 
probability. 
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Results and Discussion 

Table 1 and figure 1 above show the production indices of the three vermicomposts 
produced in this study. Vermicompost weight shows a statistically significant variation among 
the composts at maturity, Mw with 301g is higher in biomass and statistically different from 
Tp with 204 g which was also significantly different from MwTo which weighs 249g. The 
percentage of live worms is highest in Mw as a substrate at 64.5, followed by 24.7% for 
MwTp. The MwTp seemed to be the most favorable substrate to the worms, as the average live 
worm weighed 59.1g showcasing statistical supremacy over My = 40g, and Tp =20.8g. In 
terms of worm mortality, however, the highest value indicates a worse substrate for the 
worm's survival and in that, single Typha Tp proved to be the worst. At the same time, MwTo 
was the best attesting to the magnitude of the worms’ average weight in grams. 


Table 1. Vermicompost production indices 


Vermicompost Vermicompost Live Earthworm Mortality 
type weight (g) worms (% ) weight (g) rate % 
Mw 301° 64.5° 40.0° 30.0° 
Tp 204° 55° 20.8° 66.9° 
MwTo 249° 24.7° 59.1° 28.0° 
Mean 251.3 31.6 39.9 48.3 
LSD 15.79 5,02 4.08 18.8 


Mw= Municipal Waste Vermicompost; Tp = Typha Grass Vermicompost; 
MyTo = 50% Municipal Waste + 50 % Dried Typha Grass 


Production indices of the vermicomposts 
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Figure 1. Vermicomposts indices during production 
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Looking at the production indices, the two vermicomposts that contain municipal 
waste (My) performed better in this study. This could be probably due to its constitution of 
many species of biodegradable wastes sourced largely from table drops and household wastes, 
thereby supplying a variety of substrates for the worms to thrive. (Garg and Gupta, 2009; Lim 
et.al., 2015) mentioned a similar possibility as a reason for production and nutrient level 
disparity in agro-industrial waste. Also, Type grass as a substrate may be rich in Carbon and 
Hydrogen but less in peptides, which may diminish its performance as against Mw containing 
vermicompost (Ahirwar & Hussain, 2015) 


Table 2. Vermicomposts reaction, primary nutrients (NPK), and Organic Carbon status 


Parameter Mw Tp MwTo 
pH in H20 6.85 7.26 7.37 
EC (dS/m) 0.35 0.36 1.10 
Total Nitrogen —N (%) 1.95 0.80 0.80 
Phosphorous — P (mgKg’') 0.47 0.35 0.37 
Potassium —K (mgkg”') 0.72 0.49 0.4 
Organic Carbon (%) 2.80 1.20 1.80 


Mw= Municipal Waste Vermicompost; Tp = Typha Grass Vermicompost 
MyTo = 50% Municipal Waste + 50 % Dried Typha Grass 


Soil Health Indices of the vermicomposts 
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Figure 2. Soil health Indices of the produced Vermicomposts 
Table 2 and figure 2 shows the reaction and primary nutrients (NPK) status of the 


vermicomposts produced in Table | above. Both Tp and MwTo vermicomposts are neutral in 
terms of reactivity (pH), which attests to a good quality soil amendment and provides the 
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favorable conditions needed for plant and Rhizosphere conditions for good production. 
Percentage Organic Carbon (OC) in the vermicomposts produced is at least 1% higher in Mw 
than either of Tpand MwTo, 2.80%, 1.20%, and 1.80% was the respective percentages for Mw, 
Tp, and MwTo, The recorded variations in the selected soil health indices could be as a result of 
the nutritional constitution of the major ingredients of the vermicomposts as reported in 
(Barik et al., 2011). Also, the OC content of all the produced vermicomposts, is way above 
the average OC content of the soils of the study area reported in many literature (Alhassan et 
al., 2018; James & A, 2019; Sani et al., 2019; Wikipedia, 2020). This is an indicator of OC 
supplementing power of all Vermicompost produced (despite different levels) when and if 
incorporated to the soil. 


Conclusion 

In conclusion, the use of organic fertilizer especially Vermicompost has a great impact 
on improving the physical, chemical, and biological properties of the soil. Most importantly, 
Vermicomposting is the easiest way of conserving the soil and reducing environmental waste. 
Findings in this study establish the superiority of organic wastes from Municipal dump site. 
The product from such fibers and organic materials produced the best set of Vermicompost 
production indices and Vermicompost-soil health indicators. 


Recommendation 

Taking this research as a pilot, we recommend subsequent studies that will research 
the substrate type and rate, for a more informed compost recommendation and the undertaking 
of organic farming in the Sahel region. 
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Abstract 

Pesticide was considered complex substances intended for preventing, destroying, or 
controlling any pest including vectors of human or animal diseases. They are specific 
chemicals, which contain both “active” or energetic and “inactive” or inert ingredients. The 
term pesticide includes herbicide, insecticide, nematicide, sanitizer, bactericide, insect 
repellent, animal repellent, antimicrobial, fungicide, disinfectant, molluiscide, and avicide 
among others. They are classified based on the mode of entry, mode of function and action, 
and their chemical composition. These pesticides played a key role in controlling and 
mitigating variety of agricultural pest and diseases in farm, garden. They contrary, affect soil 
functional service and great biological biodiversity including human and his environment. 
Many factors and processes are involved in the transfer, transport and degradation or 
biochemical conversation of pesticides in soil. The soil factors are the soil organic matter, soil 
pH, temperature, nutrients, moisture, infiltration rate, aerobic and anaerobic, soil texture and 
structure. The combined processes are physical, biological, and chemical with great 
conversion on weathering, hydrolysis, oxidation-reduction, photolysis humification, 
cheluviation, organic sorting, translocation, eluviations, illuviation, calcification, 
podsolisation, gleying, and leaching. Seven important subjects have been treated in a synopsis 
manner — major groups of pesticides, pesticide toxicity (acute and chronic effects), processes 
affecting pesticides in soil, pesticides and its ecological effects (soil, crop, biota and human), 
step by step procedure and soil sampling, techniques involving analysis of pesticides residues 
in soil, and bioremediation and biodegradation approaches. 


Keywords: pesticides: groups, processes, effects, techniques and bioremediation 


Introduction 

Sustainable agriculture means ensuring food security and better livelihood in rural and 
urban areas, and this depends on the broad diversity of unnaturally produced chemicals, 
including pesticides, inorganic fertilizers, and bio-chemicals (Cooper and Dobson, 2007; 
FAO-WHO, 2014; European Union, 2021). The word pesticide comes from two terms namely 
‘pest’ and ‘cide’. Pest is a name used to define any plant or animal or soil organism which is 
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detrimental or harmful to soil productivity, agronomic plant growth, pastures, forest shrubs 
and fishes (Usman, 2020). On the other hand, the term ‘cide’ is a Latin word that means ‘to 
kill?’ or ‘to destroy’ (Usman, 2020). Hence, the term ‘pesticide’ can be considered as a 
complex name that includes all chemicals used to kill or control pests either in agricultural 
field or in other environments such as store rooms, human houses and gardens. According to 
FAO (2002) pesticide can be defined as any substance or mixture of substances intended for 
preventing, destroying, or controlling any pest including vectors of human or animal diseases, 
unwanted species of plants or animals causing harm during, or otherwise interfering with, the 
production, processing, storage, or marketing of food, agricultural commodities, wood and 
wood products, or animal feedstuffs, or which may be administered to animals for the control 
of insects, arachnids or other pests in or on their bodies. It may also be regarded as a chemical 
or biological agent that deters, incapacitates, kills, or otherwise discourages pests in 
agriculture; and the target pests may include insects, plants pathogens, weeds, mollusks, birds, 
mammals, fish, nematodes (round worms) and microbes that destroy property, cause nuisance 
or spread diseases, or disease vectors (FAO, 2002). Kortekamp (2011) also viewed pesticides 
as any substance or mixture of substances intended for controlling pest and diseases in crop 
production. 

Twenty-seven years ago (1995 to date), the global pesticide consumption is believed to 
have reached 2.6 million metric tons of active or energetic ingredients; however, this 
consumption of pesticide has now (2023) reached 3.53 million metric tons (BRC, 2023). The 
global commercial market of the active pesticide chemicals, also reached $90.5 billion in 
2022 and attained an increased 8.8% ($98.42 billion) in 2023 (BRC, 2023). Majority of these 
pesticide chemicals, are marketed illegally. According to the Transnational Alliance to 
Combat Illicit Trade in agrichemicals and pesticides (TRACIT, 2019), illegal pesticides 
commercially marketing globally, can be considered as outdated and unauthorized pesticides, 
pesticides without appropriate regulations, pesticides without tests and licenses, fake pesticide 
products, un-labelled and fake-labelled pesticide products, and refilled pesticide containers 
and packages found anywhere in the world. 

The major criminal enterprise regarding illicit pesticides marketing in Africa is serious 
and alarming (Medina-Pastor and Triacchini, 2020; United Nations, 2020). Typically, a global 
operation called Silver Axe was reported to have seized 1346 tons of illicit pesticides, which 
were estimated to worth up to €94 million (equal to NGN56400000000), during just the first 
four months of 2020 (Europol, 2020). 

Existing knowledge indicates that the agricultural use of pesticide chemicals 
contributes greatly to global ecosystem contamination (FAO-WHO, 2018, 2019; Sharma et 
al., 2019; Bonvoisin et al., 2020; Brescia, 2020; Lin et al., 2022; Geissen et al., 2022). With 
an increased pesticide use, questions are rising on potential effects regarding public health, 
soil and water contamination and environmental ecosystem (United Nations, 2018; FAO- 
WHO, 2020; Braga et al., 2020). This has led to an extensive review of the global pesticide 
legislation and the scale of challenge in reaching the global harmonization of food safety 
standards, which highlighted the presence and future consequences of pesticides to global 
ecosystem and human resources (Mohring et al., 2020). Thus, selection of pesticides to reduce 
human and environmental health risks was highly recommended in agriculture (Jepson et al., 
2020). And, pesticide management approach towards protecting the safety and health of 
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farmers is binding (Mohammad et al., 2018). Similarly, assessing the farmers’ safety and 
protective behaviour to use pesticides in the farms, was also looked to be necessary 
(Moradhaseli et al., 2017). 

The risk assessment of pesticide molecules is surface soil, water components and plant 
biomass including animal feeding on the plant biomass has been also demanded in many calls 
(Ngolo et al., 2019). Mapping the impact of pesticides trading is equally important (Tsimbiri 
et al., 2015; FAO-WHO, 2014; European Union, 2017; United Nations, 2020). This paper 
focused on background concept of pesticides in soil, and addressed the following subject 
areas: major groups of pesticides, processes and factors affecting pesticides in soil and 
microbial degradation, effects on soil organisms, crop damage and human health, techniques 
involving analysis of pesticides residues in soil, and bioremediation approaches. This 
overview, demonstrated and addressed the key issues regarding pesticide and its diverse 
health issues to soil, plant, human and environment. It is believed that the information has 
relevance to National Agency for Food and Drug Administration and Control (NAFDAC), 
and National Environmental Standard and Regulations Enforcement Agency (NESREA), all 
based in Nigeria. Likewise, it would undoubtedly serve as relevance materials to both 
academic and public health centers, globally. 


Major groups of pesticides 

Pesticides are specific chemicals, which contain both “active” or energetic and 
“inactive” or inert ingredients. The former is referred to ingredients, which are sparkling and 
capable of preventing, destroying, repelling or mitigating a pest of plant and human vectors, 
present in all pesticide products (US-EFA, 2017). They are considered as plant regulators, 
defoliants, desiccants, or nitrogen stabilizers, and behave very active in soil causing severe 
consequences to biological organisms (Usman, 2020). These active ingredients are classified 
as conventional (pure chemicals), antimicrobial (mixtures of substances), and biopesticides 
(biological pesticide) (US-EFA, 2017). On the other hand, the inactive ingredients are 
stationary and immobile in soil, useful for product performance and usability, and contain at 
least one active ingredient and other calculatedly added inactive ingredients (US-EFA, 2017). 

According to Drum (1980), pesticides can be grouped based on three classes, namely, 
function and action, mode of entry and chemical composition. Fishel and Ferrell (2013) noted 
that pesticides can be classified based on their function on target pest into the following 
classes, as outlined by Usman (2020). 
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Table 1. Pesticides class, description and examples (after Usman, 2020) 


No. Pesticide class Description Example 
These are pesticides formulated to kill damage or 
alleviate one or more species of insect found to be a 
1 Insecticides problem to crops. Functions include the disruption of the Aldicarb 
pest nervous system, and damaging their exoskeletons. 
Kill insects and arthropods. 
This class of pesticides produced mainly to kill fungi and 
2 Fungicides other related organisms such as blights, mildews, molds Azoxystrobin 
and rusts. 
These are produced specifically to kill bacteria and acts 
3 Bactericides —_ against any bacteria that are parasite to plant or biological Copper complexes 
organisms useful to soil productivity. 
rn Herbicides These groups of pesticides ass produced to destroy weeds Mae 
and other unwanted plants in farms and gardens. 
5 Kea es ae are pests that kill se and other related Wiad 
organisms feeding on plants and animals. 
6 iedentaiixe Ths group of pesticides is peated to control rodents, Warfarin 
mice and other related organisms. 
7 Algaecides These are used to control or kill growth of algae. Copper sulfate 
8 Larvicides These are used to inhibits the growth of larvae. Methoprene 
9 Desiccants These act as reagent on plants by drying their tissues. Boric acid 
10 Repellents This ches of pesticides is used to repel pests by its taste Methiocarb 
or smell in farm and garden. 
u Ovicides This class of pesticides s designed to hinders the growth Henan 
of eggs of insects and mites. 
Dv Virucides These pesticides are used to control viruses in soil and Gaytown 
plant. 
This class of pesticides is used to inhibit or kill mollusc ’s 
13 Molluscicides _—_i.e. a snail that usually disturbing growth of plants or crops Metaldehyde 
in farm and garden. 
4 Nematicides These ake mainly ae to kill soil nematodes that act Aldicarb 
as parasites of plants in farms and garden. 
15 Avicides These pestenes are used to kill birds, which attack rice, per 
pearl millet, wheat and sorghum. 
Lo These are used to target larvae of lampreys, which are _—_‘Trifluromethyl 
16 Lampricides ig as : . : 
jawless fish-like vertebrates in the rivers and oceans. nitrophenol 
7 Piscicides These pesticides are formulated to act against parasite Wecnone 
fishes. 
18 Silvicides "These are used to acts against woody vegetation. Tebuthiuron 
19 Termiticides _ These are used to kills termites. Fipronil 


According to Yadav and Devi (2017), classification of pesticides based on mode of 
entry, can be understood as systemic (non-contact absorbed by plant), non-systemic (contact), 
stomach poisoning and stomach toxicants, fumigants and repellents. The non-contact is 
absorbed by the plant and soil, the contact enters the body of pests through their epidermis 
and causes death, the stomach poisoning enters via mouth and causes death, the fumigants 
enter the body of pests through their tracheal or respiratory system and caused death, and the 
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repellents served as disturbance agent to pest, and are capable of deterring the pests from the 
farm (Usman, 2020). 

By and large, classification based on chemical composition of a given pesticide is 
generally the most important way of understanding pesticide behavior in soil and plant 
(Handford et al., 2015; Usman et al., 2017). Normally, this classification, gives the knowledge 
about the physical and chemical properties of pesticides, as well as the active ingredients 
included (Yadav and Devi, 2017). It also, provides the efficacy and extent of pesticides 
toxicity and the information that is useful in determining the mode of pesticide application, 
safety measures, and immediate clarification in case of any accident (Usman, 2020). Four 
groups of pesticides, were reported namely, organochlorine, organophosphates, carbamates, 
and synthetic pyrethroids pesticides. Thanks to the works of Yadav and Devi (2017) and 
Usman (2020) who have provided a comprehensive overview of these three groups of 
pesticides as presented below. 

1. Organochlorine pesticides: These pesticides are used in 1940s and 1960s for the 
control of mosquito and agricultural pests. They are chlorinated hydrocarbons, which refers to 
organic compounds attached with five or more chlorine atoms (chlorinated aromatic 
molecules). Examples are DDT, methoxychlor, dieldrin, chlordane, lindane, endosulfan, 
benzene hexachloride, kepone, mirex, toxaphene and aldrin. These chemicals worked 
effectively by disrupting the entire nervous system of the insects leading to convulsions, 
paralysis and death. 

2. Organophosphates pesticides: These pesticides are widely used in farms and 
garden, homes and offices, grazing lands and sport areas. They are the derivatives of 
phosphoric acid characterised as biodegradable compounds with less soil contamination 
effects. They worked as stomach, contact and fumigant poisons leading to nerve damage, 
paralyse and death of an insect. Examples are diazinon, chlorphyrifos, carbophenothion, 
endothion, parathion, glyphosate and malathion. These pesticides are highly toxic to insects 
but are relatively minus to human and domestic animals such as cows, goats, sheep and 
camels. 

3. Carbamates pesticides: These pesticides are derived from carbamic acid primarily 
to kill insects. The mode of action is shorter compare to that of organophosphate but interfere 
in the same way in the nervous system of insects and as stomach, contact and fumigant 
poisons. Examples are carbaryl, carbofuran, propoxur and aminocarb. 

4. Synthetic pyrethroids pesticides: These are pesticides with special active 
ingredients, which can be synthesized by duplicating the structure of natural pyrethrins 
“‘pyrethrum-like’. They are similar to natural pyrethrins that are produced by the flowers of 
pyrethrums (Chrysanthemum cinerariaefolium and C. coccineum), but with longer residual 
effects than natural pyrethrins. They are highly toxic to insects and fish but slightly toxic to 
mammals and birds, widely considered to be amongst the safest insecticides for use in foods 
and on grains. Examples are Cypermethrin and Permethrin. 

Pesticide toxicity: acute and chronic effects 

Toxicity of pesticide, length and magnitude of exposure, has been considered as the 
key determinants of the degree of harmful impact on soil organisms and human health 
(Lorenz, 2009). This can be understood in two forms namely acute and chronic toxicity. 
According to Yadav and Devi (2017) acute toxicity is the ability of a chemical pesticide to 
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cause harmful effects in human or biological organism within a very short period after 
inhalation or contact, whereas chronic toxicity is the ability of a chemical pesticide to cause 
adverse health effects which follows from long-term exposure to a chemical; and this is 
expressed in terms of lethal dose 50% (LDs0o) or lethal concentration 50% (LCs0). This LDso is 
regarded as the single exposure dose of the poison per unit weight of the organism required to 
kill 50% of the test population, where the population is genetically homogeneous expressed in 
milligram per kilogram (mm/kg) body weight (Usman, 2020). According to Singh and 
Mandal (2013) acute effects occurs as a results of accident and misused, which follows due to 
untaught concept. Pesticide Action Network (PAN, 2012) reported that the major ssymptoms 
that occurs from the chronic effect include birth defects, fetus deformity, genetic changes, 
blood and nerve disorders, endocrine disruption, and reproduction abnormalities. Continuous 
exposure to chemical pesticide particularly during application without protective clothes may 
cause chronic illness. This has been the fact that the joint FAO/WHO Meeting on Pesticide 
Residues (JMPR) and the Joint FAO/WHO Meeting on Pesticide Specifications (JMPS) has 
led to the compilation of broad lists of pesticides which can serve as a guide to the concept of 
pesticides in the field of agriculture and environment in general. The information may be 
useful to both academic environment, National Agency for Food and Drug Administration 
and Control (NAFDAC), and National Environmental Standards and Regulations 
Enforcement Agency (NESREA) operating in Nigeria. ~ 

Processes and factors affecting the transfer and transport of pesticides in soil 

The processes that involved in the pesticides communication with soil and bio- 
chemical conversion, are complex and dynamics. Usman et al. (2017) considered these as 
sorption, adsorption, physical and chemical processes, and microbial degradation. According 
to Usman (2020) as noted from the detailed study of soil formation and soil profile 
development (Jenney, 2009), numerous processes are involved in the transfer and transport of 
pesticide molecules in soil and plant. These include the physical, chemical, biological and 
ecological processes, which are set of mechanisms, considered complex in soil formation 
(Usman, 2013). It is evident that when pesticide get contact to soil, the following processes 
play a key role in the potential degradation of pesticides (Ngolo et al., 2019). These soil 
processes can be understood from the following viewpoints (Waugh, 2000; Usman, 2020): 

1. Weathering: This is a breakdown of soil particles by physical, chemical and 
biological factors. During this process, pesticide molecules can be further broken, and affect 
the soil materials. 

2. Humification: This is a process of decomposition of organic matter into humus. 
Pesticide molecules attached to the organic materials, can be included into this decomposition 
and affect the soil particles. 

3. Cheluviation: This is the transport of chelates or nutrients and organic acid 
downwards through the soil profile. Molecules of pesticides are said to be transported into the 
various parts of soil profile during this process. 

4. Organic sorting: This is the reorganization of minerals and organic matter into the 
horizons. The arrangement may contribute to the great distribution of pesticide molecules in 
soil voids. 

5. Translocation of soil materials: This is the movement of soil components in a form 
of solution, suspension, biological activities or direction (downward and upward). 
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Translocation is one of the most serious soil process that moves pesticide molecules easily 
from one area to another by means of solution, suspension and/ biological agent. 

6. Eluviations: This is a process of washing out of materials — the removal of organic 
and mineral matter from the A-horizon. Pesticide molecules that decompose with organic 
material during humification, can be washed out to A-horizon. 

7. Illuviation: This is a process of in-washing — the re-deposition of organic and 
mineral matter in B-horizon. Similarly, pesticide molecules that are decomposed along with 
organic material during humification, can be washed-into B-horizon. 

8. Calcification: This is a process of calcium accumulation in B-horizon — occurs 
where precipitation is equal or slightly high than evapotranspiration (tropical areas). During 
this accumulation of calcium, molecules of pesticide can be carried away by rainfall to other 
part of the surface soil. 

9. Podsolization: This is the removal of iron and aluminium oxides together with 
humus — occurs where precipitation is greatly more than evapotranspiration (cool or humid 
areas). During this removal process, molecules of pesticide can be removed by the rainfall to 
other parts of the surface soil. 

10. Gleying: This occurs when the output of water from the soil system is restricted, 
giving anaerobic or waterlogged soil conditions. Particles of pesticide can be deposited along 
with gley materials. 

11. Leaching: This is a removal of soluble material in soil solution — occurs where 
precipitation exceeds evapotranspiration while soil drainage is good. Pesticide molecules can 
be leached into the soil profile zones, and contaminate the soil. 

Soil organic matter, soil pH, temperature, nutrients, moisture, infiltration rate, aerobic 
and anaerobic, soil texture, and structure, are considered the most common factors affecting 
and controlling the effectiveness of soil-applied pesticides as well as the bio-chemical 
degradation in soil (Katagi, 2004). For example, Gold et al. (1996) noted that in the bio- 
chemical conversion of pesticides, soil pH and clay content are found to affect the persistence 
of bifenthrin, chlorpyriphos, cypermethrine, fenvalerate, permethrin, and isofenphos. This 
conversion was noted to be faster in the silt loam soils compared to sandy-loam and sandy soil 
for soil affected by imidacloprid and diazinon (Jones and Ananyeva, 2001). The persistence 
and biodegradation of these two pesticides in soil together with pencycuron degraded faster in 
coastal saline soil than in alluvial soil and soil amended with decomposed cow manure (Hafez 
and Thiemann, 2003). Racke et al. (1994) noted that herbicides such as atrazine and 
trifluraline disappeared more rapidly under anaerobic conditions than under aerobic 
conditions. However, this might be very limited due to the absorption process that prevents 
the release of pesticides into the soil solution and the loss by causing dislocation during 
leaching and runoff (Pal et al., 2005). 

Pesticide degradation in soil 

The term ‘degradation’ as applied to the concept of pesticides in soil means 
breakdown of pesticide molecules in soil solution. Katagi (2004) considered this breakdown 
as a process of pesticides dissolving in soils, which results in the conversion of its molecules 
into simpler compounds such as H20, C02, and NH3, Buchel (1983) noted that this breakdown 
might further be increased as a result of chemical reactions such as solubility, volatilization, 
adsorption and desorption, hydrolysis, photolysis, microbiological catabolism, and 
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metabolism; considered the major route of mineralization (Usman, 2020). These sets of 
processes are controlled by three modes of degradation in soil namely: biological (breakdown 
by micro-organisms), chemical (breakdown by chemical reactions, such as hydrolysis and 
redox reactions), and photochemical (breakdown by ultraviolet or visible light). Microbes in 
soil used pesticides as a food source under biotic and abiotic transformation processes 
(Gavrilescu, 2005). Usman (2020) noted that the biotic conversion is mediated by 
microorganisms, while abiotic conversion embroils the processes such as chemical and 
photochemical reactions. The key factors involved in these two conversions mostly in soil- 
water conditions, are hydrolysis, oxidation-reduction and photolysis (Burns, 1975). However, 
hydrolysis and redox-reactions are the main regulators in sediment conditions, with redox 
reactions mediating in aquatic environments by direct or indirect photolysis (Burns, 1975). 
Similarly, hydrolysis and redox reactions are considered the major processes controlling the 
breakdown of pesticides in soil liquid and solid forms (Usman et al., 2017). Soil moisture, 
temperature, aeration, pH, and organic matter are considered the key components of microbial 
degradation in soil (Drum, 1980). 

According to Usman (2020), surface and sub-surface soil medium are the positions 
where chemical degradation of pesticides occurs although very high at the surface medium. It 
has been described that soil temperature, moisture, pH, aeration, light and the ability of 
pesticide molecules to adhere and adsorb to soil particles, are the major factors enhancing the 
chemical degradation of pesticides in soil (Drum, 1980; Usman, 2020). Photochemical 
degradation, on the other hand, is an abiotic process that describes the dissipation of 
pesticides where molecular excitation by absorption of light energy results in various organic 
reactions (Katagi, 2004). This form of pesticides degradation in soil is influenced by the 
intensity and spectrum of sunlight, length of exposure, and the properties of the pesticide 
molecules (Katagi, 2004). 


Pesticides and its ecological effects: soil, crop, biota and human concerned 

Pesticides as chemical agent is used for preventing and destroying pest to achieve high 
crop yield and sustainable food security (Handford et al., 2015). However, these chemicals 
have more negative impact the key components of the ecosystem: soils, crops, water, forest, 
biota, and human (Zikankuba et al., 2019; Ndungu et al., 2019; Bonvoisin et al., 2020; 
Brescia, 2020). The national and international concern about their various effects remains an 
area of interest (Yao, 2020). Soils, crops, water and the human environment are four 
components, that must be put into consideration with regard to pesticides and ecological 
effects. According to Joko et al. (2020) illicit pesticides, are highly toxic and might poison 
almost all components of human and biological lives including water bodies, soil voids, 
aquatic organisms and the overall ecosystem. Usman (2018) emphasized that soil quality, 
plant health, water quality and human health, must be protected from the chronic effects of 
pesticides. These four subject areas have been regarded as the major target of pesticides when 
released into the environment for a particular purpose (Lai, 2017; United Nations, 2018; 
FAO-WHO, 2019). The pesticides molecules released during spraying in the field, 
immediately come into contact with soil and plant residues (Bollag et al., 1992; Braga, 2020; 
Lin et al., 2022; Geissen et al., 2022). The soil and its biological components (parent 
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materials, biota and biodiversity, water), the plant and biomass, and also, the atmosphere and 
air circulating, are at high risk of becoming contaminated by the synthetic molecules of 
pesticides (Yao, 2020). 

Soil biota play a vital role in soil quality and soil fertility development, and plant 
growth as well (Usman et al., 2016). The interaction between these biological organisms and 
pesticides has effects on their biodiversity, and soil-crop health relationship (Usman et al., 
2017). Sebiomo et al. (2011) noted pesticides decrease the microbial population and cause 
significant damage to soil and plants. Pesticides such as butachlor were reported to suppress 
the microbial population and deter their enzyme activities in soil (Min, et al., 2007; Xia et al., 
2012). Other pesticides were observed to have severe effects on crops, and are capable of 
destroying great number of non-target organisms (Constenla et al., 1990). These severe effects 
reduce the functional services in soil system, eliminate the biological diversity, and damage 
plant systems at a very early stage of growth (Sharma et al., 2019). According to Kalia and 
Gupta (2004), pesticides disrupt the food web and put the entire ecological system at very risk 
of altering the biosynthetic mechanism, affecting the protein synthesis and cellular membrane, 
and disorderly plant growth. Study by Mayeetreyee et al. (2013) noted a significant reduction 
of heterotrophic aerobic bacterial count and fungal population in soil treated with paraquat 
and glyphosate pesticides. Similarly, a group of herbicides such as butachlor, pyrazosulfuron, 
paraquat and glyphosate have been observed to cause damage to plant biomass and delaye a 
substantial microbial function on the proper decomposition of organic matter (Bolter et al., 
2006). Other studies of recent analyses in this regards include the recovery of selected drugs 
and herbicides from soil (De Carlo et al., 2015), the analysis of five neonicotinoid insecticides 
and their primary metabolite in cucumbers and soil (Abdel-Ghany et al., 2017), determination 
of dioctyl diethylenetriamine acetate residues in soil (Gong et al., 2018), determination of 
polyoxin b in cucumber and soil (Song et al., 2019), monitoring of 218 pesticide residues in 
clay loam soil (Acosta-Dacal et al., 2021), Simultaneous multi-residue pesticide analysis in 
southern Brazilian soil based (Do Amaral et al., 2022), and multi-residue method for trace 
analysis of pesticides in soils (Résch et al., 2023). 

The synthetic pesticides molecules, play a major role in contaminating environmental 
ecosystems, and affect many aspects of human health developments as well (Haggblade et al., 
2019). They are poisons and can be hazardous when misused in the farm or garden (FAO- 
WHO, 2019; United Nations, 2020). They poisoned almost all components of the human body 
including hairs, skin, eyes, sliver, and clothes (Lai, 2017; Joko et al., 2020). Many studies 
reported different pesticide cases on human and environmental health impact (Khan, 2009; 
Handford et al., 2015; van Wendel de Joode et al., 2016; Larsen et al., 2017; El-Haoud et al., 
2017; Mohammad et al., 2018; Mora et al., 2018; Cappa et al., 2022; Marks-Perreau et al., 
2023). The complex environmental and human health effects of pesticides reported by these 
studies are largely due to the commercialization of the illicit chemicals (European Union, 
2021). The toxicity level and formulation of the energetic ingredients contained in these 
outdated pesticides are produced without test and quantification (United Nations, 2018; 
TRACIT, 2019). This toxicity of pesticide molecules has detrimental effects on human health 
(Lorenz, 2009). Diseases such as brains (neuroblastoma) cancers, soft tissue sarcomas, 
colorectal and testes carcinomas; hormonal imbalances leading to infertility, breast pain, 
menstrual disturbances, adrenal gland exhaustions and early menopause among others, have 
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been linked to high level of exposure to toxic pesticides (Larsen et al., 2017; van Wendel de 
Joode et al., 2020). The common symptoms reported include eyes and skin irritation, body 
shivering, headaches, body aches, skin rashes, poor concentration, nausea, dizziness, impaired 
vision, cramps, panic attacks and in severe cases coma and death could probably occur 
(Mohammad et al., 2018; Joko et al.,2020). The human birth effect, infants and children 
multiple diseases cases are also noted as a result of pesticide use in agriculture (Larsen et al., 
2017; Mora et al., 2018; Widyawati et al., 2020). This has also linked to problems related to 
inability to speak fluently, hormonal imbalances leading to infertility, breast pain, menstrual 
disturbances, adrenal gland exhaustions, early menopause, immune system dysfunction 
leading to immune suppression that cause potentially serious health risks and also subject to 
malnutrition (Xavier et al., 2004; Tsimbiri et al., 2015). These complex health issues are 
considerable (WHO, 2017; Zikankuba et al., 2019), and need to be on a regular check and 
balance assessment particularly in northern Nigeria where many rural farmers are mingling 
with pesticides on daily times (Usman, 2020). 

Techniques involving analysis of pesticides residues in soil 

The human relationship with soil is a long-term friendship and has recorded many 
benefits and detriments as well. The entire soil body represents one of the most complex 
mediums with diverse functional services and interactions of pesticides that take place 
through various processes and compound factors (Bollag et al., 1992; Usman and Kundiri, 
2016). The use of pesticides in soil has infiltrated the existence of many residues, which are 
unfriendly to the soil itself, and biological organisms living within (Pelosi et al., 2021; 
Geissen et al., 2021). The employment of soil-based techniques for the analysis of these 
pesticide residues in soil has become necessary (Leo et al., 2016; Usman, 2020). This is to 
provide open access to short- and long-term worldly risk assessment of eco-toxicological 
effects of pesticides in soil bodies and biological lives including humans and _ the 
environmental ecosystem. In this development, many techniques and ideas have been used to 
achieve this risk assessment. The use of traditional techniques, which employed the 
application of examples of liquid—liquid extraction, Soxhlet extraction, sonication assisted 
extraction among others, are laborious, expensive, time-consuming, and require large amounts 
of organic solvents, which also involve many steps, leading to loss of some analyte quantity, 
consequences with the use of hydrocarbon solvents such as depletion of the ozone layer, and 
producing of carcinogenic waste (Mondal, 2020). 

However, considering the historical catastrophes of pesticides in creating soil 
contamination and soil pollution (Usman et al., 2017; Usman, 2018), the scientific 
advancements, has provided some alternative techniques that are quick, easy, cheap, effective, 
sharp, and safe to replace previous, less efficient extraction methods for pesticide analysis 
(Varela-Martinez et al., 2020). These alternative techniques, include the use of solid phase 
extraction (SPE), solid phase micro-extraction (SPME), supercritical fluid extraction (SFE), 
microwave assisted extraction (MAE), microwave-assisted micellar extraction (MAMEB), 
accelerated solvent extraction (ASE), matrix solid phase dispersion (MSPD) extraction, and 
QuEChERS (quick, easy, cheap, effective, rugged and safe) method (Gonzalez-Curbelo et al., 
2022). These methods have been developed for a purpose of overcoming the shortcomings of 
the traditional techniques, even though some of them (e.g. SFE, ASE and MAE) are 
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instrumental techniques, useful for the desired purpose such as purification (SPE) and 
concentration (SPME) of obtained extracts (Mondal, 2020). 

Step by step procedure and soil sampling 

The general step by step procedures for the removal of pesticide residues from the soil 
is believed to begins with an advanced design of the experiment or choosing the best and 
friendly techniques. This will be followed by sampling, extraction, clean-up and estimation. 
In this regard, soil samples should be taken using pedon idea as introduced by Soil Survey 
Staff (2022) or from growing farm using the idea of grid pattern that should consistently 
distributed all over the intended farm (Mondal, 2020). For practical interest, Mondal (2020) is 
suggesting 3 x 3 grid with nine total sample proportions for smaller fields, 4 x 4 with sixteen 
sample propositions for the medium-sized fields, and 55 with more than sixteen 
proportions for larger grids and very large fields. According to him each sample site stands as 
a reference of one portion of the total sample. This means that soil sample has to be taken at 
each site with two soil plugs of about 15 cm deep and 3 to 5 cmp; whereas the combined two 
plugs at each site will become a one sample portion (Mondal, 2020). 

Examples of some modern techniques of pesticide residues analysis 

1. QUEChERS technique: This technique means quick, easy, cheap, effective, 
rugged, and safe (QUEChERS) approach, and was first introduced in 2008 by Lesueur et al., 
(2008) after publishing their work titled ‘Comparison of four extraction methods for the 
analysis of 24 pesticides in soil samples with gas chromatography-mass spectrometry and 
liquid chromatography-ion trap-mass spectrometry’. Since its known procedure and 
techniques, the QuEChERS has received the attentions of many studies year after year 
(Lesueur et al., 2008; Pinto et al., 2009; Rashid et al., 2010; Caldas et al., 2011; Temir et al., 
2012; Rouvire et al., 2012; Fernandes et al., 2013; Chai et al., 2014; De Carlo et al., 2015; 
Zhou et al., 2016; Abdel-Ghany et al., 2017; Gong et al., 2018; Song et al., 2019; Chen et al., 
2020; Acosta-Dacal et al., 2021; Do Amaral et al., 2022; Résch et al., 2023). Significant 
achievement has been made in the analyses of organochlorine, organophosphorus pesticides, 
pyrethroid pesticides, neonicotinoids], carbamates, triazole, and urea derivatives, with the 
advanced QuEChERS (Chai et al., 2014; Abdel-Ghany et al., 2017; Ma et al., 2020; 
Mohamed et al., 2021; Résch et al., 2023). Many advanced and modified procedures were 
developed using QUEChERS more especially when analysing different pesticides in soil and 
plant (Lehotay et al., 2007; Mohamed et al., 2021; Rajput et al., 2021). 

2. Multiple Parameter Hazard Assessment Scoring System (MPHASS): Multiple 
Parameter Hazard Assessment Scoring System (MPHASS) is a technique that employed the 
ideas of Chemical Scoring System (CSS) introduced by O’ Bryan and Ross (1988) as designed 
by US-EPA office of toxic substances and the Oak Ridge National Laboratory. This technique 
combines multiple parameters such as human, soil, water, plant and animal components of 
environmental ecosystem. The soil parameters include the water in soil, soil materials in 
illuviation zone, soil materials in eluviation zone, soil materials in rhizosphere zone; soil 
biota: bacteria, fungi, ant, termites and nematodes. 

3. A handheld SRS microscope: A handheld SRS microscope is a direct method, and 
can be used for pesticide analysis in both soil and plant (Figure 1). The technique is a modern 
method, and can detect the pesticide residue in situ, with high speed, under ambient light, and 
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without special sample preparations, providing a convenient examination platform with a 
label-free detection of the pesticide residue (Lin et al., 2022). 
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Figure 1. Typical handheld SPS microscope in situ detection of pesticide on crop* 

* - This image is cited from Lin et al. (2022), and was appeared as: Fig. 37.5. In situ detection of 
pesticide residues on crop product by a handheld SRS microscope. (A) Digital picture of the experimental setup. 
(B) SRS image at 1605 cm- 1 of a spinach leaf. (C) SRS image at 1679 cm- 1 of a spinach leaf. (D) MCR 
output spectra. (E, F) SRS concentration maps of chlorophyll and thiabendazole. The image was Reprinted with 
permission from C.-S. Liao, et al., In vivo and in situ spectroscopic imaging by a handheld stimulated Raman 
scattering microscope, ACS Photonics 5 (3) (2018) 947-954. Copyright 2021 American Chemical Society). 


4. Multi-residue method via QUEChERS techniques: This technique employed the 
idea of QUEChERS in a multi-purpose manner for the analysis of pesticide residues in soil. 
Rosch et al. (2023) have accurately quantify 146 currently used pesticides in agricultural soils 
with varying soil properties using a multi-residue method for trace analysis of pesticides in 
soils with special emphasis on rigorous quality control. According to Résch et al. (2023) 
pesticides were extracted using an QUEChERS technique while chemical analysis was carried 
out by liquid chromatography attached to tandem mass spectrometry (triple quadrupole). 
Similarly, quantification was done based on matrix-matched internal standards calibration, 
using 95 isotopically labeled analyte analogues; the validated was conducted by an in-house 
prepared partly aged soil, which contains all target pesticides and by agricultural field soils 
with native pesticide residues (Résch et al., 2023). They considered their multi-residue 
technique as highly sensitive, precise and true for regular monitoring of soils contaminated 
with pesticide residues. 

5. Solid Phase Micro-Extraction (SPME) technique: This is one of the newly 
developed technique widely used for the pesticide residues analysis in soil. The technique is 
set to run concurrently for purification and concentration of the sample extract. According to 
Mondal (2020), the SPME syringe is used as the main part of this method that is visually 
resembling on the chromatographic system contained a | cm long fiber which is located 
within the needle of the syringe. This syringe needle is made of an appropriate polymer 


Bulgarian Journal of Soil Science® 2024 Volume 9. Issue 1. www.bsss.bg 103 


Usman et al., 2024 


deposited on the holder of fused silica. However, Mondal (2020) noted that the micro- 
extraction process is fixed based on the rearrangement of extract between micro-extraction 
fiber and sample matrix which he described as thermal in the case of the gas chromatography 
or by solvent elution in the case of liquid chromatography. 


Conclusions 

This evaluation discussed the concept of pesticide in soil with reference to the general 
concept of pesticide and its major groups, pesticide toxicity (acute and chronic effects), 
processes affecting pesticides in soil, pesticides and its ecological effects (soil, crop, biota and 
human), step by step procedure and soil sampling, techniques involving analysis of pesticides 
residues in soil. An intended study reflecting on pesticides analysis in soil, will be enlighten 
with the information discussed in this paper. As other similar reviews have confirmed, the 
content of this paper, is an appraisal towards better understanding of pesticides and how it 
affects the health condition of soil, plant, and human ecosystem. Analysis of pesticide in soil 
with modern techniques, is a good step towards better sustainable management of best use of 
pesticide in farms and garden. This is also vital towards better understanding of the complex 
effects of pesticides in soil and human environment. Some examples of the ways to be used 
for sustainable analysis of pesticide residues in soil are provided. These examples can be 
served as a guide towards better understanding of pesticide residue analysis in soil and plant. 
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Abstract 


The understanding of soil water retention characteristics, including field capacity, 
plays a pivotal role in various agricultural studies and practices. This knowledge informs the 
development of irrigation and drainage schedules, assessments of soil water storage capacity 
(commonly referred to as plant available water), analyses of solute movement, evaluations of 
plant growth, and assessments of water stress levels. This study the aim of creating a 
regression model or equation to predict field capacity of soils in a typical Alfisol, which 
originated from parent material of basement complex, at the University of Ilorin Teaching 
and Research Farm. A total of forty-five (45) disturbed and forty-five (45) undisturbed soils 
samples were collected along a toposequence (upper, middle and bottom slope) at depths: 0 — 
30 cm, 30 — 60 cm and 60 — 90 cm. Soil properties of the disturbed and undisturbed samples 
were ascertained through standard laboratory techniques and/or calculated utilizing 
established protocols. The measured soil properties include sand, silt, clay, bulk density, total 
porosity, field capacity and organic matter. A significant negative relationship between field 
capacity and bulk density (r = -0.582, P<0.05) was found and a significant positive 
relationship between field capacity and total porosity (r = 0.581, P<0.05) was found. A total 
of four different models were developed for predicting field capacity of the soil. Model four 
of field capacity was identified as the best model with the highest R* adjusted value of 
0.4494. The model explains 45% of variance in the mean square errors of field capacity with 
sand, bulk density, silt and organic matter contributing statistically to the model. 


Keywords: field capacity, regression model, basic soil properties, A/fisols 


Introduction 


Field capacity refers to the maximal water content present in soil after two or three 
days since saturation, with minimal drainage occurring. This parameter is a measure of 
maximum amount of water soil can retain for plant use (Hillel, 1998). Historically perceived 
as a constant soil characteristic, field capacity plays a critical role in assessing water retention 
and supply properties, determining irrigation requirements, and evaluating soil quality (Lal, 
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2006; Shao et al., 2006). Optimal plant growth typically coincides with soil water levels 
maintained close to or at field capacity throughout the plant's growth cycle. Therefore, 
maintaining soil moisture levels near field capacity is essential for effective irrigation 
management. As water resources become increasingly valuable, accurately measuring field 
capacity becomes essential for implementing efficient soil and water management strategies 
(Madhumita and Verma, 2011). 

In the last few years, more researchers have undertaken the creation of models that 
simulate soil processes. These models serve the dual purpose of enhancing comprehension of 
crucial soil dynamics and serving as valuable tools for assessing agricultural and 
environmental challenges. As a result, simulation models have become commonplace in both 
research and management endeavors (Minasny and McBratney, 2002). 

Soil moisture characteristics such as Field capacity (FC) is a fundamental soil 
characteristic essential for understanding plant water availability, infiltration, drainage, and 
solute movement. It plays a critical role in various soil water studies. However, directly 
measuring FC can be expensive, time-consuming, and laborious due to the high variability of 
soil properties across space and time. Therefore, scholars have developed alternate methods 
to estimate FC using more readily obtainable soil properties. These properties often include 
proportions of soil particle sizes (sand, silt, and clay), bulk density, total porosity, and organic 
matter (OM) content. Many of these estimation methods are referred to as pedotransfer 
functions (PTFs) as they translate existent, easily measured data into soil hydraulic data 
(Bouma, 1989). While PTFs are usually developed based on experimental observations, their 
pertinency can be limited to the specific datasets used in their creation (Donatelli et al., 1996; 
Wosten et al., 1999). Additionally, different PTFs can generate significantly diverse 
estimations for the same soil. This variability makes it challenging for users to select the most 
suitable PTF for their peculiar purposes (Acutis and Donatelli, 2003). Furthermore, PTF 
estimations often have a moderate level of accuracy. Ideally, PTF predictions would 
incorporate measures of reliability to enhance user confidence (Schaap et al., 2001). 

Soils of the study area have been classified as Typic Haplustalf. The region exhibits 
a mosaic vegetation pattern, alternating between forest and savanna. Gravelly Alfisols are the 
predominant soil type, reflecting their formation over the basement complex. 

Rawls et al. (1982) presented regression equations for the valuation of the soil water 
content at different water potentials in soils with different particle size distribution. The 
general form of the linear regression equation is 
Q, =atb (sand =) +c (siie2) +d (clay2) +e (organic matter 2) +f (bulk density 2) +9 (Total porosity=.) 


Where a, b, c, d, e, f and g are regression coefficients. 

This study has as an objective to create a regression model or equation for predicting 
field capacity of soil on a typical Alfisols at the University of Ilorin Teaching and Research 
Farm. 
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Materials and Methods 

Study Area 

This study took place along a toposequence at the University of Ilorin Teaching and 
Research Farm, Ilorin, Nigeria. The toposequence included upper, middle, and bottom slope 
positions. The dominant vegetation in the location is an intermix of forest and savanna, and 
the underlying soils formed over basement complex (Olaniyan, 2003). Gravelly Alfisols are 
the prevalent soil type. While the land has been cultivated in the past, it was fallow during the 
sampling period. The study area is situated within the Southern Guinea Savanna zone of 
Nigeria at about 4° 35' East longitude and 9° 29' North latitude. The elevation is 307 meters 
above sea level. The climate is tropical and characterized by bimodal rainfall featuring peaks 
in June and September and a break from mid-July to August. Precipitations rate ranges 
between 1000 mm to 1240 mm per annum. Daily temperatures typically fluctuate between 
20°C and 35°C (Kolo et al., 2012). 

Soil Sampling and Laboratory Analysis 

We collected soil samples at three depths 0-30 cm, 30-60 cm, and 60-90 cm along a 
toposequence. Fifteen small pits were dug, with five pits located at each position on the 
toposequence. This resulted in a totality of 90 samples collected (45 undisturbed and 45 
disturbed). At each depth, undisturbed core samples were collected using metal cylinders 
with a height of 8.3 cm and an interior diameter of 5.5 cm. The soil was secured with a piece 
of calico cloth held in place with a rubber band. The samples were then labeled for 
identification. Disturbed soil samples were also collected at each depth using a soil auger. 
These disturbed samples were placed in labeled polyethylene bags. All samples were carried 
to the laboratory for further analysis of their physical and chemical properties using standard 
methods. 

Soil Samples Preparation 

The disturbed samples underwent a process of air drying, crushing, and sieving 

through a 2 mm mesh before being utilized for analysis. 


Soil Properties Analysis 

Particle Size Analysis 

Particle size analysis was conducted utilizing the hydrometer technique outlined by 
Gee and Or (2002), with the inclusion of sodium hexametaphosphate (calgon) as a 
dispersant. 

Bulk Density 

Bulk density was measured using the core method with oven drying at 105°C to 
constant weight (Blake & Hartge, 1986). The calculation followed the standard mass- 
volume relationship: 


Ms 


i (2) 
Where, p, = bulk density (~2), 


Ms = Dry soil mass (Kg), and 
Vt = Total volume of soil (m?>) 
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Total Porosity 
Total porosity was calculated as: 


eas | aid (3) 


p 


s 
Where,@ = Total porosity CS) 
p; = bulk density C3), 


p, = particle density asummed to be 2650 kg/m? 

Soil Organic Matter (OM) 

Soil organic carbon content was determined by the Walkley-Black wet oxidation 
method (Nelson & Sommers, 1982). Following standard practice, the percent organic carbon 
was converted to estimate total soil OM by multiplying by a factor of 1.724, which reflects 
the assumption that OM is composed of approximately 58% carbon (Brady & Weil, 1999). 

Field Capacity 

The determination of field capacity was done by saturating the core samples with 
water for 24 hours, weighing the samples after three days, and then oven-drying them at 
105°C. 


Weightof saturatedsoil after three days—oven dried weight 


Percent water at field capacity = x 100% 


(4) 


oven driedweight 


Statistical Tool 

Correlation Analysis 

In statistics, when two or more variables change together in a predictable way, they 
are considered correlated. 

A multiple correlation was used to investigate the relationship between measured pair 
of values. The statistical tool is given as: 


ex (Zi) (Zi) (5) 
nm 
Rij —- 
and Z;;= - : (6) 


j 
is the standardized values of the response variable Xi. 


L(Xij - x) 
= j 
= 2S (7) 


is the sample variance of the response variable Xi 


This correlation coefficient is generally called Karl Pearson’s coefficient of 
correlation. Correlation coefficient measures the degree of association (linear relationship) 
between two or more random variables. Correlation coefficient may be positive or negative. 
A significant positive correlation coefficient implies that an increase in one variable is 
accompanied by an increase in the other variable. While significant negative correlation 
coefficient implies an inverse relationship between the two random variables in question. In 
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contrast, a non-significant correlation coefficient indicates zero relationship between the 
variables under study. 

A test statistic suggested by Morrison (1976) shall be used to test for significance of 
this correlation coefficient. This test is given as: 


t=r |= (8) 


i-r? 
where, r is the so called Karl Pearson’s correlation coefficient and n is the sample 
size. It can be shown that when Hy is true, t © tx/2 (nm — 2) (9) 


Model Development 

Field Capacity 

It should be noted that in multiple regression the predictions (independent variables) 
are assumed to be known without any uncertainty in their given values. Our model can be 
found useful if some or all of the assumptions hold on the datasets depending on the purpose 
for which we intend to use it. Because we like to use our models for both estimation and 
predictions we will investigate all the assumptions and correct necessary violations in the fit 
of the models. 

Majorly, we will use diagnostic plots from the regression models to check 
assumptions including normality of residuals and presence of potential outliers using 
"Residuals versus Fitted" charts to show if there is a trend to the residuals and Shapiro-Wilk’s 
test of residuals normality (Shaporo and Wilk, 1965) in the R software for statistical 
computing and graphics (Rcore team, 2013) to confirm residuals normality. Model residuals 
(error terms) are assumed to be normally distributed if the residual plot shows a random 
scatter pattern with a consistent spread of residuals around zero for each of the six predictor 
variables, it suggests a good model fit (Olorede et al., 2013; Olorede and Mudasiru, 2013). 
When the points in a quantile-quantile plot, especially those in the central region, fall closely 
along a diagonal line, it suggests that the data are likely normally distributed. We will use 
standardized residual criterion to check effects of potential outliers in the dataset on the 
models by removing observations with standardized residuals outside the interval -2.2 
according to Barnett and Lewis, 1994; John and Prescott, 1975; and Stefansky, 1972. To 
assess potential outliers in our models, we will calculate the minimum and maximum 
residuals for each model. We will then standardize these residuals by dividing them by the 
corresponding model's root mean square error (RMSE). Ideally, the standardized residuals 
should fall within a specific range. Values outside this range (expected to be between +1 and 
+2) might indicate potential outliers that require further investigation. 

All statistical analyses, model fitting and diagnostics shall be done using the language R 
version 3.0.2 (R Core Team, 2013). 

For multiple linear regression models of this type, statistical hypothesis tests 
regarding the model's parameters can be used to assess the model's effectiveness. We will 
therefore describe and test hypotheses about parameters of the proposed regression model 
collectively as well as individually. This will allow us to know if there is significant 
relationship between the response variables (field capacity) and a subset of the predictors. We 
will screen the built models based on number of significant parameters, maximum amount of 
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proportion of variability about the response by these parameters and satisfaction of model 
assumptions using diagnostic plots as aforementioned. 

Rejection of null hypothesis about the full model (model with all parameters) suggests 
that one or more of the predictor variables significantly impacts the model. Rejection of null 
hypothesis about individual regression parameter (ANOVA table) indicates that the variable 
cannot be deleted from the model. Coefficient of multiple determination (R’) measures 


extent of decrease in the variability of field capacity gotten by using the 6 predictors in the 
model. High R’ alone doesn't guarantee a good regression model. (Myers and Montgomery, 
1995). Adding more predictor variables will generally inflate R°, even if those variables don't 
meaningfully contribute to explaining the data's variance. This can lead to models with high 
R’ that perform poorly on unseen data. To address this issue, some researchers prefer using 
adjusted R° (Myers and Montgomery, 1995). 

Generally, R°agj can actually decrease when adding irrelevant variables to the model. 


A marked difference between R” and R2, jis a high possibility of the inclusion of irrelevant or 


Statistically insignificant terms. We are also aware of the fact that presence of 
multicollinearity in the data set will make estimates of coefficients from the least squares fit 
imprecise and Statistically insignificant (Multivariate Calibrations test by Naes and Martens). 
If the primary aim is simply to predict the dependent variable (Y) using the X variables, 
multicollinearity might not be a critical issue. The model's predictions for Y could still be 
accurate, and the overall R? (or adjusted R’ aaj) will reflect how well the model fits the data for 
Y value prediction. However, multicollinearity becomes a significant concern when seeking 
to understand how individual X variables influence Y. Individual P-values for regression 
coefficients can be misleading. A high P-value might suggest an unimportant variable, even if 
it contributes to the overall prediction. The confidence intervals around the regression 
coefficients become very wide. These intervals may even include zero, making it difficult to 
determine if an increase in an X variable leads to an increase or decrease in Y. Due to the 
wide confidence intervals, removing or adding a data point can significantly change the 
coefficients and even reverse their signs. 

Sometimes, multicollinearity can lead to seemingly paradoxical results. The overall model fit 
(indicated by a low overall P-value) might be good, yet individual X variables lack statistical 
significance in predicting Y. This occurs because highly correlated X variables provide 
redundant information. While neither variable may contribute significantly on its own, 
together they improve the model's fit. Removing them would significantly worsen the overall 
prediction. Therefore, multicollinearity can lead to a well-fitting model for prediction 
purposes, but hinder the interpretation of the individual X variables' effects on Y. 


Test of Hypothesis about Full Regression Model 

Hypothesis 1: [Test of significance encompassing all regression parameters] 

Null hypothesis (H, ): The model fits the data/ The model is adequate 

Alternative hypothesis (H,): Not H, 

Or mathematically, 

Null hypothesis (H,): 6, = 8, =-**=8,;=0 [No single predictor showed 
significant contribution to the model] 
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Alternative hypothesis (Hy): : 8; + 0 for at least one j [At least one predictor 
showed significant contribution to the model] 


Tests statistic: F = “Regression Goal F-test] 


ratio MSerror 
Decision rule: Reject null hypothesis and accept alternative hypothesis at 0.05 
significance level if P — value < 0.05, or else do not reject null hypothesis. 


Hypothesis 2: [Test of significance encompassing individual regression coefficient] 

Null hypothesis (H, ): 6; = 0 [Predictor x; is not significant given that others are 

included in model] 

Alternative hypothesis (H,): : 6; +0 [Predictor x; is significant given that 

others are included in model] 

value — qa 
jie" 


A 
= diagonal element of the covariance matrix corresponding to B; 


Tests statistic: t [Individual t-test] 


Cy 
o* = variance of B; 

Decision rule: Reject null hypothesis and accept alternative hypothesis at 0.05 
significance level if P — value < 0.05, or else do not reject null hypothesis. 


Results 

The results of the descriptive statistics obtained for soil properties is presented in 
Table 1. Results obtained for the Pair-wise correlation of variables measured during field 
capacity experiment is presented in Table 2. Results of parameter estimates for the models 
generated for field capacity is presented in Table 3. Results of the Analysis of Variance 
(ANOVA) for the multiple regression models created for field capacity are presented in Table 
4. 


Table 1. Descriptive Statistics for Soil Properties 


Soil Properties Mean Min Max Std dev C.V 
Sand (g/kg) 636.96 400 885 140.15 22.00 
Silt (g/kg) 198.73 31 270 70.37 35.41 
Clay (g/kg) 166.53 82 350 88.38 53.08 
Bulk density (kg/m?) 1402.36 1065 1698 140.19 10.00 
Total porosity (m*/m*) 0.471 0.359 0.598 0.053 11.25 
OM (g/kg) 13.12 1.38 27.80 7.40 56.40 
Field capacity (%) 44.05 32.26 71.43 11.89 26.99 
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Discussion 

Statistical Analysis 

Descriptive Statistics Analysis 

The results obtained in Table 1 for the values of sand, silt and clay contents, bulk 
density, total porosity, OM and field capacity ranges from 400 to 885, 31 to 270, 82 to 50 
g/kg, 1065 to 1698 kg/m, 0.359 to 0.598 m*/m*, 1.38 to 27.80 g/kg and 32.26 to 71.43%, 
respectively. 

Pair-wise correlation 

Sand had negative correlation with silt (r = -0.870*) and clay (r = -0.909*) (Table 2). 
As sand increases, silt and clay will decrease. Sand was positively correlated with OM (r = 
0.668*) as sand increases, OM will increase. Silt was positively correlated with clay (r = 
0.603*) and inversely correlated with OM (r = -0.512*) this shows that silt affects OM 
negatively and as silt increases, clay will also increase. Clay had negative correlation with 
OM (r = -0.689*) clay affects OM negatively. Bulk density had negative correlation with 
total porosity (r = -0.999*) and field capacity (r = - 0.582*) as bulk density increases total 
porosity and field capacity decreases Total porosity was positively correlated with field 
capacity (r = 0.581%) as total porosity increases, field capacity will increase. 

Model Development 

Four models with p-values < 0.05 were created for field capacity with the 
corresponding R? adjusted values of 0.3282, 0.3849, 0.4288 and 0.4497 for models 1, 2, 3 and 
4 as shown in Table 7. Model 4 of field capacity was selected as the best predictive model 
due to having the highest R’ adjusted value of 0.4497. This implies that the model fits the 
field capacity data with the four predictors included after removal of observations with 
potential outliers and non-significant predictors clay and total Porosity so as to determine the 
model’s prediction accuracy. The model justifies approximately 45% of the variation in mean 
squared errors of field capacity with sand, bulk density, silt and organic matter being 
statistically significant contributors to the model. 

Examining the significance of each predictor variable in Table 4 for model 4 of field 
capacity, the results show that we can reject the null hypothesis for bulk density, silt and OM 
with p-values of 6.722e-06, 0.01404 and 0.04252, respectively. The simple interpretation of 
this is that each variable is statistically significant within the model when considered 
alongside the other included variables. This suggests they contribute meaningfully to the 
model's explanatory power in combination. It does also mean they cannot be excluded. The 
expression of model 4 as given in Table 3 is expressed in equation (10) as: 


Y = 153.787791 — 0.065790 * Sand — 0.113318 * Silt — 0.037695 * Bulk Density 
+ 0.487612 * Organic Matter (10) 


For model 4 of field capacity, the residuals are now normal and there are no potential 
outliers as shown in Figure 1. In fact, the Cook’s distance plot now confirms this by having 
all the Cook’s distances for all the predictors less than 1 as shown in Figure 1. The residuals 
plot and normal quantile-quantile plots presented in Figure | also support this. Normality of 
the residuals is further confirmed by the histogram presented in Figure 2. 
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The model equation generated for model 4 was used in predicting the field capacity 
values shown in Table 5. Table 5 presents observed, predicted and residuals of the prediction. 
The values obtained for observed and predicted field capacity as presented in Table 5 were 
analyzed for extent of relationship and presented a correlation coefficient of of 0.7086. Based 
on the guideline outlined at http://www.westgard.com/lesson42.htm for assessing correlation 
coefficients, it is observed that when r falls within the range of 0.90 and 1.00, 0.70 and 0.89, 
0.50 and 0.69, 0.30 to 0.49, and 0.00 to 0.29 are said to show very high, high, moderate, low, 
and little if any correlation, respectively. It indicate that field capacity predictions of model 4 
of field capacity have high correlation relationship with the observed field capacity values. 

The graph of observed and predicted field capacity values is presented in Figure 3. 
This shows the trend of observed and predicted values of field capacity. 
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Fig. 3. Graph of Observed versus Predicted Field Capacity 
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Table 5. Summary of Model 4 Prediction of Field Capacity 


Observations Observed Y (%) Predicted Y (%) Residuals (Y — Y) 
1 32.26 33.07635 -0.81635 
2 32.26 37.02047 -4,76047 
3 38.46 48.16551 -9.70551 
s 32.26 33.96041 -1.70041 
5 32.26 32.50082 -0.24082 
6 57.69 50.70623 6.983774 
a 71.43 63.37098 8.059022 
8 32.26 38.74145 -6.48145 
9 38.46 49.08697 -10.627 
10 32.26 38.98002 -6.72002 
11 38.46 44.33006 -5.87006 
12 38.46 51.89833 -13.4383 
13 38.46 46.34247 -7.88247 
I 37.31 32.24192 5.068075 
15 37.31 27.7532 9.556797 
16 57.69 49.4147 8.275303 
17 57.69 52.42 5.269999 
18 57.69 49.63758 8.052422 
19 57.69 53.22095 4.469047 

20 57.69 53.8081 3.881896 
21 32.26 37.5815 -5.3215 

22 38.46 48.55325 -10.0933 
23 48.39 33.43684 14.95316 
24 38.46 47.08808 -8.62808 
25 32.26 40.62179 -8.36179 
26 57.69 49.80101 7.888989 
27 57.69 49.68176 8.008239 
28 57.69 45.74357 11.94643 
29 57.69 45.29467 12.39533 
30 38.46 40.567 -2.107 

31 32.26 41.06191 -8.80191 
33 57.69 49.85784 7.832161 
34 48.39 34.08673 14.30327 
35 32.26 35.84598 -3,58598 
36 38.46 46.11329 -7.65329 
37 38.46 43.909 -5.449 

38 38.46 39.81466 -1.35466 
39 38.46 41.38441 -2.92441 
40 38.46 41.76351 -3.30351 
4] 38.46 41.77301 -3.31301 
42 38.46 38.70256 -0.24256 
43 38.46 39.62808 -1.16808 
44 32.26 28.61995 3.640048 


Note: Observations 32 and 45 are not included because they contain potential outliers 
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Conclusion 

In the experiment carried out at the University of Ilorin Teaching and Research Farm 
for developing a model for predicting field capacity of the soil, a total of four distinct models 
were created. Model four of field capacity was chosen as the best model, having the highest 
R’ adjusted value of 0.4494. Sand, silt, bulk density and OM was found a fitting option to 
ascertain soil moisture at field capacity. 
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Abstract 


An investigation was carried out on selected properties of basaltic soils along the 
toposequence of Aya river in Ogoja Cross river state. The soils are deep and occur on a gentle 
and undulating landscape with a gradient of <4%. Three profile pits were dug in each of the 
mapping units identified as crest, mid-slope and valley bottom. Each Pedon had four distinct 
genetic horizons (Ap, AB, Btl and Bt2). The dominant matrix for the Ap horizon was dark 
brown (7.5YR3/4), brown (7.5YR6/4) and yellowish brown (10YR6/8) for pedons I, II and 
II respectively. The particle size distribution data indicates the preponderance of sand 
fraction while the percentage sand mean values are 82.45 %, 81.45 % , 85.78 % and 81.0 % 
for pedons HI, II and I respectively. The soils reaction pH p20 is assessed as strongly acidic 
ranging from 5.25-5.75. The exchangeable Ca**, Mg”* and K* are in low concentrations but 
are within the threshold of critical limits. Organic carbon and total nitrogen and available P 
were generally low (N: 1.7-1.1 g/kg; OC:7.3-5.1 g/kg and P:18.25- 12.60 g/kg). The Organic 
C/N ranged from 3 — 7 and Ca/Mg had mean of 1.40, 1.36 and1.31 in pedons I, Il, and HI 
respectively. The base saturation (%BS) was high to very high (64.70-79.26 %) while 
effective cation exchange capacity was generally low (2.03 cmol/kg, 2.29 cmol/kg, and 1.93 
cmol/kg for pedons I, II, and III respectively). The soil on the toposequence of Aya river 
showed evidence of oxic properties with absence of weatherable minerals in the sub stratum. 
The soil suitability class were moderately suitable with limitations in pH and low soil organic 
matter content. The soils were classified as Isohyperthermic ustalfs which corresponded to 
Arenic Luvisols (pedons 1 and 2) with fluvic characteristics in pedon 3 (WRB, 2006). Soil 
quality of the studied soils is low due to the established unfavourable fertility indicators; 
hence soil conservation and use of integrated soil fertility management- organic sources 
(compost manure), inorganic sources (NPK fertilizers) and crop rotation as well as 
conservational tillage is recommended. 


Keywords: basaltic materials, toposequence, morphology, classification, suitability 
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Introduction 

Soils formed on basaltic parent materials are usually rich in ferromagnesian minerals 
mainly plagioclase feldspar (which are rich in Ca and/or Na- CaAl,Si,O4 or NaA13Si2O0) and 
pyroxene [(Mg, Fe)2SiO,] (Esu et al., 2015) and (Ojo-Atere et al., 2011). Esu et al., (2015) in 
their investigation in the morphological, physicochemical and mineralogical properties of 
soils developed from basalt at not too far away Ikom, they established that 77.7% of the 
mineralogy is composed of kaonilite, 10.75% is Fe bearing minerals (goethite and hematite) 
and they concluded that these soils are dominated by highly resistant minerals, a sign of soil 
in an oxic state of weathering where leaching, eluviation-illuviation, laterization and faunal 
pedoturbation are common features (Esu et al., 2015 and Eshett, 1987). In many tropical 
soils, where excessive rainfall and leaching of the bases results in the acidification of the soils 
through Fe- Al saturation in the exchange site while Nsor and Akamigbo (2015), also 
observed the problem of plinthization in the soils of central cross river state. 

Topography depends on the configuration of the system in the of elevation, slope and 
landscape expositions. Differences in topography can cause wide variations in soils within 
the confines of a single field. Topography determines the local distribution or disposal of the 
precipitation and determines the extent to which water tables influence soil genesis (Foth, 
1990). Topography hastens or delays the work of climatic factors in pedogenesis with respect 
to soil horizon development (Brady and Weil, 1999). Toposequence is used interchangeably 
with catena by many soil scientists, and catena concepts have emanated as slope- soil 
evolutionary process (Esu et al., 2008). The soils differ as a result of erosion, transportation 
and deposition of surficial materials as well as leaching, translocation and deposition of 
chemical and particulate constituents in the soil. These processes result in morphological 
changes associated with soil colour and horizonation as a result of hydrology related to 
topographic position (Juo and Moorman, 1980). Soil colour is one of the important basic 
properties which help in identifying the kinds of soils and recognize the successions of soil 
horizons or layer in the soil profile and it gives the physical structure of the soil some 
intrinsic imprint that are related to topography. It has been used for soil identification and 
qualitative measurements of soil properties and serves as a useful tool in field studies 
(Noshadi et al., 2013). While Wakene (2003), stated soil colour is function of pH, redox 
reaction, and organic matter. A change in soil colour from adjacent soil also indicates a 
difference in the soil mineral origin or soil development, geologic origin and degree of 
weathering of the soil material and leaching or accumulation of chemical compounds such as 
iron, which may greatly influence soil quality (Fisher and Binkley, 2000). Though not all 
soils are strongly correlated with terrain attributed, much less can be linked to them in a 
straightforward functional way (Sobieraj et al., 2002), this is even more complex within 
agricultural field where different types of land use types are practiced. Landscape position 
influences rainfall, drainage and erosion, water velocity on a slope affects deposition of 
materials in suspension, the largest sized particles like sand are the first to drop out of 
suspension, fine clay size particles can be carried further away from the base of the slope 
before they are deposited. Hill slope orientation also affects the microclimate of a place, 
inclined surface facing into the sun tend to be warmer and drier than flatter surface facing 
away from the sun (Atofarati et al, 2012). Soils developed on a toposequence have been 
found to be suitable in the production of crops such as tree crop (Theobroma cacoa and 
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Gmelina aborarea), food crops (yam, cassava and maize) and vegetables (tomato, fluted 
pumpkin and garden egg etc). These were reported by Eshett (1987) for basaltic soils on the 
toposequence at Ikom in Northern Cross river state, Nigeria. Soil in the tropics are extremely 
diverse, they ranged from young fertile to senile soils. The objectives of this study were to 
characterise and classify the basaltic soil along the toposequence of Aya river in Ogoja Cross 
River state and study its suitability classification for selected crops. 


Materials and Methods 

Site Description 

The study area was along the toposequence of Aya River in Ogoja Local Government 
Area of Cross River State. It lies between latitude 6°41.38’N and longitude 8°58.03’E. The 
area has distinct rainy (March to October) and dry seasons. The mean annual rainfall in the 
study area ranges from 1728 - 2200 mm with double maxima in the months of July and 
September. Mean annual temperature varies from 22.4° to 32.5°C while the mean annual 
relative humidity varied from 74 to 79%. The climate of the area is sub-humid tropical, and 
its vegetation is forest-grassland-savanna. The soils are generally derived from the basaltic 
and alluvium of either false-bedded sandstone or the upper coal measure (Asadu, 1990) with 
gentle undulating plains with scattered rock outcrops. The hydrology is governed by Aya 
River (it raises from Obudu plateau) which is a major tributary of Cross River. The soils are 
well drained with few poorly drained soils (loamy sand to sandy clay) which are used for 
paddy rice and sugarcane cultivation. Its vegetation is characterized as forest- grassland- 
savannah (transitional vegetation type). Agriculture is the major socioeconomic activities in 
the area with cottage industries in fishing and agricultural related enterprises. The major farm 
produce is yam (Dioscorea spp), plantain (Musa paradisiaca), cassava (Manihot esculenta), 
rice (Oryza sativa), ground nut, maize, soybean, bambara nut (Vigna subterranean) and 
vegetables (water leaf and fluted pumpkin- Telifairia occidentalis). The area is rich in forest 
resources but with the rate of current deforestation amidst uncontrolled access. Oil palm 
(Elaeis guineensis) and cocoa (Theobroma cacoa) are common but value addition to oil palm 
fruits is also a major cottage industry alongside cassava processing. 


Field Work 

Three mapping units were identified and classified based on physiographic position 
along the toposequence and associated landuse. A stratified random sampling technique was 
adopted in the study. The profile pits were dug at upper, mid-slope and lower part of the 
slope. A total of three profile pits was sunk and described following the guidelines of FAO, 
(2006). Soil samples were taken for routine physical and chemical analysis. 


Laboratory Analysis 

The pH was determined by glass electrode pH meter both in 1:2.5 soil/ liquid 
suspension of water and KCl (Henderson et al., 1993). Particle size analysis was by 
hydrometer method (Gee and Or, 2002). Organic carbon was determined by wet dichromate 
method by Udo et al., (2009). Total nitrogen was determined by micro-Kjeldahl digestion 
method as modified by Udo et al., (2009). Extraction of available phosphorus was done using 
Bray 1 method. Exchangeable cations (K*, Ca, Mg**, and Na‘) were extracted by neutral 
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normal ammonium acetate, K* and Na’ in the extraction were determined by flame 
photometer (Udo et al.,, 2009) while Ca’* and Mg** were by atomic absorption 
spectrophotometer. Cation exchange capacity was by summation method. Exchangeable 
acidity was determined in 1NKCI extracting solution with 0.5N NaOH using phenolphthalein 
indicator by titration method of Mclean as described by Udo et al., (2009). Bulk density was 
measured by core method (Grossman and Reinsch, 2002). Total porosity (Po) was obtained 
from bulk density (Cp) values with assumed particle density (Cs) 2.65 g cm” as follows, 
Porosity (Po) =100 — (tp/€s) <100/1 

Statistical analysis was performed using Genstat statistical package for the analysis of 
variance and treatment mean compared using least significant difference at 5% level of 
probability. 


Table 1. Location Coordinates and Elevation 


Position Profile code Elevation (M.A.S.L) North-East coordinates 

Lower fat 110 6 39°27.9"N, 8 59’23.0°E 
Mid-slope I 108 6 39°27.3"N, 8 59’21.8"E 
Upper (Crest) I 119 6 39°27.3"N, 8 59’20.4”E 


Results and Discussion 

Morphological properties of the studied soils are shown in table 2, the physiographic 
mapping units have different colour matrix range, and all the horizons are deep and well 
drained characterized by loose top soil. The peds are generally arranged in weak to moderate 
granular and sub angular blocky, clear wavy horizon boundaries and soils have fine and 
medium roots and pronounced faunal activities in the epipedons. The soil at the slope appears 
to be shallow with a depth of 92 cm while the depth at both the midslope and valley bottom 
are well above 110 cm, this was in agreement with the findings of Ofem and Esu (2015). The 
soil colour varies from the crest to the valley bottom, the crest is dominated by brownish 
yellow (1OYR6/8) at 0-25 cm and strong brown (7.5YR5/8) at depth of 25-47cm, reddish 
yellow (7.5YR5/8) at 47-76 cm and 5YR7/8 at 76-92 cm. The mid-slope is dominated by 
brownish colour 7.5YR4/4 at 0-17cm, dark yellowish brown (1OYR4/6) at 17-39cm, 
yellowish red (SYR6/6) at 39-100cm. the soils at the valley bottom had dark brown 
(7.5YR3/4) at 0-20cm, strong brown (7.5YR5/8) at the sub surface profiles, the variation in 
colour matrix of the soils is a reflection of the sequence of drainage condition and 
physiographic positions and this conformed with the works of Nuga et al., (2006) and 
Akamigbo and Asadu (1986). The presence of many fine pores, many fine roots and fauna 
activities in Ap horizon of the pedons are indicative of the desirable effect of soil organic 
matter which enhances aggregate stability (Daniel, et al., 2000). Termites and earthworms 
were observed as agents of micro-faunaturbation in the study areas. 
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Table 2. Morphological Properties of Soils Along the Toposequence of Aya River 


Depth 5 Text Horizon ; 

(em) Horizon Colour ve Structure ae Consist. Veg. Root presence 

Crest Pedon 1 
Yellowish brown 

0-25 Ap (10YR6/8) LS Gr SC NSNP Many(vf-f) 
Strong brown 

25-47 AB (7.5YRS5/8) LS Gr SC NSNP Common(vf-f) 

SR 

Reddish yellow 

47-76 Btl (7.5YR6/6) SL Gr SC NSNP Few(f) 
Reddish yellow 

76-92 Bt2 (SYR7/8) SL Gr SC FRSST Very few(F) 

Midslope Pedon 2 

0-17 Ap Brown (7.5YR6/4) LS Gr CW Crumby Many(vf-f) 
Dark yellowish ' 

17-39 AB brown(10YR6/4) LS SAB CW Friable Common(vf-f) 
Yellowish Red Ext 1 a 

: ellowish Re xtremely 

39-74 BT1 (SYR6/6) SCL SAB CW firm Few(f) 
Yellowish Red Extremely 

74-135 BT2 (SYR7/8) SCL SAB CW firm Very few(F) 

Valley Pedon3 
Dark brown : 

0-20 Ap (7.5YR3/4) LS Gr CS Friable Many(vf-f) 
Strong brown : 

20-44 AB (7.5YR5/8) LS Gr CS Friable SR common(vf-c) 
Strong brown Firm and 

44-72 Btl (7.5YR4/6) S Gr CS sii Few(f) 
Strong brown Very firm 

72-110 Bt2 (7.5YR5/6) SL Gr CS and sticky Very Few/(f) 


Colour: YR = reddish yellow, Structure: Gr= granular, SAB = subangular blocky, Consistence: NSNP= non 
sticky non plastic. FRSST = friable slighty sticky, Roots: C= coarse, vf = very fine, vf = very few, f= few, f=fine, 
vegetation: SR= secondary regrowth 
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The physical properties of the studied soils are shown in table 3. The particle size 
distribution data indicates the preponderance of sand in all the pedons in the study area; this 
suggests they have quartzite in its geology. The mean values of sand are 82.45 % (pedon III), 
81.45 % (pedon II), 85.78 % and 81.0 % (pedon I). This is in contrast with the several authors 
who reported the dominance of clay in the textural class. The difference in the textural class 
may be due to intense and prolonged precipitations in the study area. 

The silt content showed no regular pattern in pedons I and [I but without a 
fluctuation the profile in pedon II. Clay content was generally higher than silt content in all of 
the studied pedons. The highest clay content was established in pedon II (17.15%). This trend 
is a common physical feature with soils in this zone and this is due to weathering of silt 
fraction into clay. The silt clay ratio (SCR) in the studied profile showed mean values of 
<1l.which indicates of reserve of weatherable minerals in the substratum as typical of 
kandiustalfs (USDA, 2015). The SCR is an index for gauging extent of weathering in soils 
(Oleghe and Chokor, 2015). 

The textural classes range from loamy sand in the top soil (Ap horizon) and sandy 
loam in the subsoil. Pedon II has loamy sand in (Ap horizon) and sandy clay loam (subsoil), 
pedon III had loamy sand (Ap horizon) and sandy loam (subsoil). The mean values of bulk 
density were 1.24 gicm®, 1.17 g/cm”, L227 g/em® for pedons I, II and III respectively. The 
bulk density decreased with profile depth though values obtained are within the range for 
good rooting crops (Soil Survey Staff, 2006 and Ladon, 1991). 

The porosity values range from 52 % - 56 %. These porosity values indicate good 
soil-air conditions. This meets with the study of Chude et al., (2011); which noted that soil 
physical properties are dominant factors affecting the usage of soils. This is because the 
success or failure of agricultural venture (crop wise) is often hinges on the physical properties 
of the soil and physical properties are more difficult to change than chemical properties 
(Okolo et al., 2015). Several authors (Sanchez, 2019) pointed out that the soil texture is less 
susceptible to change 
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Table 3. Physical Properties of Soils Along the Toposequence of Aya River 


Horizon Depth Sand Clay Silt SCR TC BD Po MC 
(cm) (%) (%) (%) g/cm’ % 
Crest (pedon I) 
Ap 0-22 86.2 11.4 2.4 0.21 LS 1.27 52 8.31 
AB 25-47 85.2 13.4 1.4 0.1 LS 1.42 46 8.48 
Btl 47-76 82.2 14.4 3.4 0.24 SL 1.33 50 9.44 
Bt2g 76-92 76.2 18.4 5.4 0.24 SL 1.05 60 14.9 
Mean 82.45 14.4 3.15 0.21 1.27 52 10.28 
SD+ 4.5 2.94 1.71 0.08 0.08 3.12 3.123 
CV 5.46 20.42 54.29 38.1 12.6 30.35 30.35 
Midslope (pedon II) 
Ap 0-17 89.2 9.4 1.4 0.15 LS 1.19 29 8.03 
AB 17-39 85.2 13.4 1.4 0.1 LS 1.29 51 6.32 
Btl 39-74 78.2 20.4 1.4 0.07 SL 1.11 58 9.16 
Bt2g 74- 73.2 25.4 1.4 0.06 SCL 1.07 60 11.7 
Mean = 81.45 17.15 1.4 0 1.17 56 8.8 
SD+ 7.14 7A4 0 0.1 0.1 2.26 2.26 
CV 8.77 41.63 0 4.45-7.55 48.55 25.68 25.68 
Valley bottom (pedon III) 
Ap 0-20 86.8 7.3 5.9 0.81 LS 1.28 52 8.9 
AB 20-44 89.8 8.3 1.9 0.23 LS 1.16 56 20.1 
Btl 44-72 91.8 7.9 0.3 0.04 S 1.25 53 6.62 
Bt2g 72- T4.7 13.1 12.2 0.93 SL 1.28 52 92 
Mean = 85.78 9.15 5.06 0.5 1.24 53 10.39 
SD+ 7.66 2.67 23 0.43 0.06 6.6 6.6 
CV 8.93 29.18 104.74 86 4.84 63.52 63.52 
LSD(5%) 10.92 30.2 39.16 16.5 0.01 17.6 59 
SCR= silt/ clay ratio, TC= textural class, BD g/cm’, Po= total porosity, MC = moisture content 
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The chemical properties of the studied soil are as presented in table 4. The pH had 
mean values of 5.50, 5.46 and 5.64 for pedons I, II and III respectively, the soil reaction 
generally ranges from strongly acidic in the topsoil to moderately acidic in the subsoil 
(acidity decreased with profile depth). The acidity may be due to leaching caused by 
precipitations and urea fertilizer application (many of the farmers are found to apply urea 
fertilizers in crop production). This indicates that the soils are strongly acidic and shows that 
there are significant exchangeable Al°** and H* hence the entire soil will have net negative 
charges on soil colloids and this will engender cation adsorption on the soil exchange 
complex (Kolay, 2002). 

According to Osodeke, (2017) such acidic cause fixing of nutrients such as 
phosphorous, potassium, and sulphur hence making them unavailable for plant uptake. The 
exchangeable Ca’*, Mg”* and K* are in low concentrations but are within the threshold of 
critical limits for deficiencies and likely in response to fertilization (Adisa et al., 2016 and 
Haby et al., 1990). The mean value for exchangeable Ca”* was 0.86 cmol/kg, 0.78cmol/kg 
and 0.82 cmol/kg in pedons I, II and III respectively. The concentration of exch. Mg** iS 
close to that of exch. Ca’* with mean values of 0.61 cmol/kg, 0.57 cmol/kg and 0.63 cmol/kg 
in pedons I, II and III respectively. Exchangeable Mg** below 0.28 cmol/kg would manifest 
as deficiency in maize (FMANR, 1990). 

Organic carbon and total nitrogen and available P were generally rated low to 
moderate (N: 1.7-1.1 g/kg; OC: 7.3-5.1 g/kg and P: 18.25- 12.60 g/kg). The obtained value of 
TN, OC and available P obtained were higher than those reported by Ofem and Esu (2015). 
These variations may be due to climatic factors and agronomy practices. 

The soil would require higher amount of N and P fertilizers to maintain their adequate 
supply. The application of gypsum is required to increase the soil pH, release the fixed 
phosphorous and boost Ca and Mg in the soil (Osodeke, 2017, Esu et al., 2015). The mean 
values of exch. K* were assessed as very low 0.01 cmol/kg in all of the studied pedons. This 
is very critical for crop development as K* is macro element in plant nutrition. Therefore, K* 
in the profile is inadequate for crop utilization (Imadojemu et al, 2017). The org. C/N ratio 
and Ca/Mg ratio were range from very low to low. 

The org. C/N is important for determining mineralization and immobilization of 
nitrogen, the C/N ranged from 3.15 — 6.64 which is below 25 favours mineralization set by 
Paul and Clark (1989). The established Ca/Mg ratio of <5.1 indicates low activity clay 
(Brady and Weil, 2014). The Ca/Mg ratio had mean values of 1.40, 1.36 and 1.31 in pedons I, 
Il, and III respectively. The percentage base saturation (%BS) was high to very high (64.70- 
79.26 %), a similar result was reported by Ofem and Esu (2015). 

The effective cation exchange capacity was assessed as low, the highest mean value 
was obtained in pedon II (2.29 cmol/kg). This means that the soil requires integrated soil 
fertility management (organic composting), inorganic (N.P.K fertilizers) and conservation 
agriculture (crop rotation, sub soil tilling etc). 
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Correlation Matrix 

The correlation matrixes are presented in table 5. There was a negative but highly 
significant correlation between sand and clay (-0.85**) while %BS and ECEC was also 
negative but highly significant correlation (-0.97**). The positive and highly significant 
correlations are the silt and clay, clay and ECEC, TN and OC. 


Table 5. Correlation matrix of selected physical and chemical Properties of the soils 
along the Toposequence of Aya River 


Sand Silt Clay SCR pH OC T™N AvP TEB ECEC BS 
Sand 1 
Silt -0.48 1 
Clay -0.85** — -0.05 1 
SCR -0.22 0.927** -0.29 1 
pH(water) -0.35 -0.13 -0.32 -0.25 1 
Oc 0.08 -0.28 0.08 -0.14  -0.06 1 
TN -0.14 0.56 -0.18 0.56 -0.27 -0.8** — 1 
AvP 0.19 -0.44 0.05 -0.42 = 0.01 -0.31 0.03 1 
TEB -0.09 0.58* -0.24 0.49 0.14 -0.51 0.67* -042 1 
ECEC -0.88** 0.31 0.81** 0.13 -0.65* = -0.23 0.33 -0.02 0.01 
BS 0.86** -0.26 -0.82** -0.09 —0.64* 0.13 -0.19 -0.04 0.19 -0.97** 1 


* Correlation is significant at the 0.05 level, ** Correlation is significant at the 0.01 


Soil Classification 

The soils of the study area were classified using USDA Soil Taxonomy system 
(2015), and correlated with World Reference Base (2006). In classifying these soils, the 
criteria considered were nature of the epipedon, the types of diagnostic master horizon, the 
clay content, the organic matter content, the percentage base saturation, the presence, or 
absence of concretions (plinthitic features), the soil moisture, temperature regimes, diagnostic 
horizons and soils colour. The soils tends to have more bright hues (reddish brown) due to 
good aeration, these pedons was classified as in the order as Alfisols because of the presence 
of Bt horizon, high percentage base saturation and the absence of lamellae, suborder as 
utsalfs owing to the Ustic moisture regime, great group as Argillic typic ustalfs as the most 
subsoil horizons were argillic along with high OC contents as well as the high %BS while the 
preponderance of sand showed psammentic (USDA) correlates to Arenic (WRB). The 
subgroup was classified Eutric arenic Haplustalfs. The family was classified coarse sandy 
clay loam, Isohyperthermic, ustalfs and classified as Ogoja sandy clay loam series. This 
corresponds to Arenic Luvisols (pedons | and 2) with fluvic characteristics in pedon 3 
(WRB, 2006). 


Conclusion 

Basaltic soils are rich in intermediate plagioclase feldspar which means they have 
reserve of weatherable minerals. The physical and chemical limitations of tropical Alfisols 
can be ameliorated by applying both conservational and good tillage practices. The soil 
texture was dominated by sand (85.78 %) followed by clay (17.15%), define texture as loamy 
sand. Soil reaction ranges from strongly acidic in the topsoil to moderately acidic in the 
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subsoil. The soil pH differ among the different pedons on the toposequence with mean values 
of 5.50 (Pedon I), 5.46 (Pedon II) and 5.64 (Pedon III). The soils had low BS (%). The total 
N was also low in concentrations. Available P was in medium concentrations and exch. K* 
and the Na’ was assessed as very low concentrations and do not pose a threat to soil 
salinization. Most prominent problems to these soils are high temperature and precipitation. 
The soil quality and fertility indicators are generally low for these soils; hence soil 
conservation and use of integrated soil fertility management- organic sources (compost 
manure), liming, inorganic sources (NPK fertilizers) and crop rotation as well as 
conservational tillage is recommended. The soils were classified as moderately suitable with 
limitations in pH and low soil organic matter content. 
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